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Schlioicn  photograph  of  the  laminar  separation  and  turbulent 
reattachiiu'nt  ori  the  hemisphere  nose  body.  Flow  is  from  left 
to  right.  Maximiim  height  of  ihe  separated  region  is  approximately 
0.01  inches.  Rej,  = 2.6  X 1 0‘ . Body  diameter  is  two  inches. 
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ABSTRACT 


The  influence  upon  the  basic  viscous  flow  about  two  axisym- 
metric  bodies  of  (i)  freestream  turbulence  level  and  (ii)  the  injt^ction 
of  small  amounts  of  a drag -reducing  polymer  (Polyox  WSR  301)  into  | 

the  test  model  boundary  layer  was  investigated  by  the  schlieren  flow  j 

visualization  technique.  The  changes  in  the  type  and  occurrence  of 

I 

cavitation  inception  caused  by  the  subsequent  modifications  in  the  j 

viscous  flow  were  studied.  A nuclei  counter  using  the  holographic  j 

I 

technique  was  built  to  monitor  freestream  nuclei  populations  and  a I 

few  preliminary  tests  investigating  the  consequences  of  different  pop-  j 

ulations  on  cavitation  inception  were  carried  out. 

Both  test  models  were  observed  to  have  a laminar  separation 
over  their  respective  test  Reynolds  number  ranges.  The  separation  on 

1 

one  test  model  was  found  to  be  insensitive  to  freestream  turbulence  j 

levels  of  up  to  3.  75  percent.  The  second  model  was  found  to  be  very  j 

susceptible  having  its  critical  velocity  reduced  from  30  feet  per  second 
at  a 0.04  percent  turbulence  level  to  10  feet  per  second  at  a 3.75  per- 

i 

cent  turbulence  level.  Cavitation  tests  on  both  models  at  the  lowest  I 

I 

turbulence  level  showed  the  value  of  the  incipient  cavitation  number  i 

and  the  type  of  cavitation  were  controlled  by  the  presence  of  the  laminar  i 

separation.  Cavitation  tests  on  the  second  model  at  0.65  percent  tur- 
bulence level  showed  no  change  in  the  inception  index,  but  the  appear- 
ance of  the  developed  cavitation  was  altered. 

The  presence  of  Polyox  in  the  boundary  layer  resulted  in  a 


cavitation  suppression  comparable  to  that  foun<J  by  other  i/jv<-sli^alors . 
The  elimination  of  the  normally  occurring  laminar  separation  on  thesi? 
bodies  by  a polymer -induced  instability  in  the  laminar  boundary  layi:r 
was  fovuid  to  be  responsible  for  the  suppression  of  inception. 

Freestream  nuclei  populations  at  test  conditions  were  measur- 
ed and  it  was  found  that  if  there  were  many  freestream  gas  bubbles  the 
normally  present  laminar  separation  was  elminated  and  travelling  bub- 
ble type  cavitation  occurred  — the  value  of  the  inception  index  then  de- 
pended upon  the  nuclei  population.  In  cases  where  the  laminar  separa- 
tion was  present  it  was  found  that  the  value  of  the  inception  index  was 
insensitive  to  the  freestream  nuclei  populations. 
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Fig.  1 
Fig.  2 

Fig.  3 
Fig.  4 

Fig.  5 

Fig.  6 
Fig.  7 

Fig.  8 
Fig.  9 


LIST  OF  FIGURE  CAPTIONS 

Cavitation  inception  number  versus  velocity  for  a series  of 
six  hemisphere  nose  bodies.  (Ref.  128) 

Results  of  a comparative  cavitation  inception  test  on  a modi- 
fied ellipsoidal  headform  sponsored  by  the  International 
Towing  Tank  Conference.  Taken  from  F'ig.  8 of  Ref.  4 
Photographs  of  incipient  cavitation  on  the  ITTC  headform  in 
various  facilities.  Taken  from  Ref.  5. 

Diagrams  of  the  test  bodies  showing  the  geometric  details. 

(a)  hemisphere  nose  body 

(b)  NSRDC  body 

Schlieren  photographs  of  the  laminar  separation  and  cavitation 

inception  on  the  hemisphere  nose  body.  The  photographs  illus - 

trate  the  development  of  cavitation  from  "steady  bubble"  type 

to  "band  type"  cavitation.  Inception  can  clearly  be  seen  to 

start  in  the  reattachment  zone  of  the  separation  . Flow  is  from 

5 

right  to  left.  Re,,  = 3.3  x 10  , o = 0.66  for  each  photograph. 
Schematic  diagram  of  the  LTWT. 

Photograph  of  nozzle  and  stagnation  section  of  the  LTWT. 
Panels  in  the  lower  right  of  the  picture  are  the  tunnel  power 
and  velocity  controls.  On  top  of  the  stagnation  section  is  the 
expansion  pot  and  the  piping  to  the  vacuum  pump  and  drain. 
Schematic  diagram  of  the  details  and  positions  of  the  turbulence 
manipulators  in  the  Low  Turbulence  Water  Tunnel. 

The  second  honeycomb  located  in  LTWT  stagnation  section. 

(a)  View  of  honeycomb  in  its  supporting  frame.  Honeycomb 
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Fig.  11 


Fig.  12 


Fig.  13 


Fig.  14 


Fig.  15 
Fig.  16 
Fig.  17 


Fig.  18 


ia  4'  X 4'  X 2"  deep.  (b)  Close-up  view  oi  I bexanorial 
cells . 

One  of  the  three  turbulence  tlampinn  screens  loi  alcd  in  lln 
stagnation  section  of  the  I/FWl  . (a)  View  of  sc:reen  on  ils 

frame.  Screen  is  4'  y 4'.  (b)  Close-up  view  of  screen. 

Screen  is  made  of  0.0075”  diameter  stainless  steel  wire  and 
has  22  mesh/lineal  inch. 

Photograph  of  turbulence  damping  screen  being  installed  in 
the  stagnation  section  of  the  LTWT. 

Photograph  of  turbulence  generating  grid  being  installe'^  in 
the  test  section  entrance  of  the  LTWT. 

Photograph  of  the  Disa  hot -film  probe  used  to  measure  tur- 
bulence levels  in  the  LTWT. 

Photograph  of  the  hot-film  probe  mounted  in  the  LTWT. 

Flow  is  from  right  to  left. 

Summary  of  the  turbulence  level  measurements  in  the  LTWT. 
Schematic  diagram  of  the  schlieren  flow  vi sualization  sy stem. 
Photograph  of  the  schlieren  light  sources.  The  silver  box  is 
a Hg  arc  lamp  used  for  steady  observation.  The  removable 
mirror  reflects  the  light  along  the  optical  axis  of  the 
schlieren  system.  The  black  cylinder  in  the  right  of  the 
photograph  is  the  spark  gap  light  source. 

The  camera  side  of  the  schlieren  system.  The  lens  simul- 
taneously focuses  the  collimated  beam  on  the  knife  edge  and 
the  image  of  the  model  on  the  film.  The  camera  box  is  about 
SOinches  long  giving  an  overall  magnification  of  about  8 times . 


Fig.  19 
Fig.  20 

Fig.  21 
Fig.  22 

^ig.  23 
-g.  24 

Fig.  25 

Fig.  26 
Fig.  27 

Fig.  28 
Fig.  29 
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Light  source  side  of  schlieren  system  mounted  at  HSW'I  . 
Camera  side  of  schlieren  system  mounted  at  IISW'J'.  System 
is  suspended  from  ceiling  to  isolate  it  from  tunnel  vibration. 
Example  of  quality  of  .schlieren  photography  obtained  in 
LI  WT  which  has  flat  glass  test  section  windows,  h'low  is 
from  left  to  right  on  hen ) isphe  re  nose  body  . Hep  - 2 . 6 ^ 1 f)  . 
Example  of  quality  of  schlieren  photography  obtained  m 
HSWT  which  has  curved  plexiglass  test  section  windows. 

Note  loss  of  detail  in  reattachment  zone  as  compared  to 
Fig.  21.  Flow  is  from  left  to  right  on  a hemisphere  nose 
body.  Rep  = 4.  2 x 10^. 

Schematic  diagram  of  the  injection  system. 

Photograph  of  the  injection  system  at  the  LTWT  showing  the 
holding  tank  and  the  flowmeter. 

Schematic  drawing  showing  the  injector  mounted  inside  the 
hemisphere  nose  body. 

Photograph  of  injector  installed  in  hemisphere  nose  body. 
Photograph  of  dismantled  injector  showing  porous  foam  and 
sintered  disc  used  to  help  smooth  the  injection  flow. 
Schlieren  photograph  of  an  example  of  stable  injection  of 
water  at  stagnation  point  of  hemisphere  nose  body. 

Uj  /U^=  0.  16,  U^=  12.  2 feet  per  second.  Freestream  flow 

is  from  right  to  left. 

Schlieren  photograph  of  an  example  of  unstable  injection  of 
water  at  stagnation  point  of  hemisphere  nose  body. 

Uj /Ug,=  0.  29,  = 1 2 . 2 feet  per  second.  Freestream  flow 
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is  from  ri)>ht  to  left. 

Fig.  30  Plot  of  the  regions  of  stable  and  unstable  injection  with 

water  on  the  hemisphere  nose  body  for  two  injection  liil>c- 
diamete  r s . 

Fig.  31  Schematic  drawing  of  holocamera:  A)  dielectric  mirror 

B)  iris  C)  spectrophotometer  cell  D)  ruby-flash  lamp  as- 
sembly E)  iris  F)  dielectric  mirror  G)  beamsplitter 
H)  neutral  density  filter  I)  1 . 6 cm  focal  length  lens  J) 
pinhole  K)  43.2  cm  focal  length  collimating  lens  L)  front 
surface  mirror  M)  pin  diode  N)  filmholder. 

Fig.  32  Schematic  drawing  of  the  reconstruction  system. 

Fig.  33  Photographs  of  the  components  of  the  nuclei  counting  system. 

(a)  holocamera  set-up  at  the  LTWT.  (b)  the  reconstruction 
system.  On  the  monitor  screen  can  be  seen  the  reconstruct- 
ed images  of  some  lOOu  diameter  polystyrene  spheres. 

Fig.  34  Schematic  drawings  showing  the  definitions  of  the  laminar 
separation  dimensions. 

Fig.  35  Effect  of  freestream  turbulence  level  upon  laminar  separa- 
tion on  hemisphere  nose  body.  Flow  is  from  right  to  left. 

Rco  = 2.6  X 10^. 

(a)  u'/U  = 0.05%  (b)u'/U=1.2% 

(c)u'/U=2.3%  (d)u'/U  = 3,6% 

Fig.  36  Effect  of  freestream  turbulence  level  upon  the  laminar  sep- 
aration on  the  NSRDC  body.  Flow  is  from  right  to  left. 

5 

Reo^  1.6  X 10  . 

(a)  u'/U  = 0.05%  (b)  u'/U  = 0.65%  (c)  u'/U  = 1.  1% 
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Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 


(d)  u'/U  = 2.3%  (e)  u7U  = 3.6% 

37  Effect  of  freestream  (urbulence  level  on  the  laminar  .separa 

tion  on  the  NSRDC  body.  Flow  is  from  right  to  left, 

5 ^ 

Rep  = 2.  5 X 10  . 

(a)  u'/U  = 0.05%  (b)  u7u  = 0.05%  (e)  u aj  1 1% 

(d)  u'/U  = 2.3%  (e)  u'/U  U6% 

38  Effect  of  freestream  turbulence  level  on  the  length  of  the 
separated  region  on  the  hemisphere  nose  body. 

39  Effect  of  freestream  turbulence  level  on  the  length  of  the 
laminar  separation  on  the  NSRDC  body. 

40  Effect  of  freestream  turbulence  level  on  the  height  of  the 
separated  region  on  the  hemisphere  nose  body. 

41  Effect  of  freestream  turbulence  level  on  the  height  of  the 
separated  region  on  the  NSRDC  body, 

42  Effect  of  freestream  turbulence  level  on  the  position  of  sep- 
aration for  both  the  hemisphere  nose  body  and  the  NSRDC 
body . 

43  Comparison  of  calculated  and  measured  values  of  the  criti- 
cal Reynolds  number  for  the  hemisphere  nose  and  NSRDC 
bodies . 

44  Effect  of  injection  of  water  on  the  laminar  separation  on  the 
hemisphere  nose  body.  Flow  is  from  right  to  left. 

Re„=  3. 9 X 10^. 

(a)  Q (injection  rate)  :=  0.0  ml/sec  (b)  Q = 3 . 57  ml/sec 
(c)Q  = 8.35ml/sec  (d)Q  = 13.21ml/sec 
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Fiji.  Effect  of  injection  of  water  on  tin-  latninar  separation  on  tin- 

NSR  DC  body.  Flow  is  from  rij>ht  to  left.  C > 1 

(a)  Q (injection  rate)  ; 0.  0 ml  /sec  (l>)  Q 1 . HO  ml  /s<-c 

(c )  Q - . 60  m 1 / sec  (d)  Q - 6 . <>()  ni  1 / se  c (e  ) Q 0 . HO  m 1 / sec 

Fig.  46  Effect  of  injection  uf  J^olyox  (WSR  101)  on  the  laminar  s«-par 
ation  on  the  hemisphere  nose  body.  Flow  is  from  right  to 

5 

left.  Rep=4.2xl0.  Concent  ration  - 600  wpprn  . 

(a)  Q (injection  rate)  =0.0,  G = 0. 0 

(b)  Q = 0.  1 ml/sec,  G = 0.47  y lO'^ 

(c)  Q = 0.3  ml/sec,  G = 1.40  y 

(d)  Q = 0.  5 ml/sec,  G = 2.34yl0’^ 

(e)  Q = 2.0  ml/sec.  G = 9.35  y lO'^ 

Fig.  47  Effect  of  injection  of  Polyox  (WSR  301)  on  the  laminar  separ- 
ation on  the  hemisphere  nose  body.  Flow  is  from  right  to 

5 

left.  Rep  = 3.9  X 10  . Concentration  = 100  wppm . 

(a)  Q (injection  rate)  = 0.  0 ml/sec,  G = 0.  0 

(b)  Q = 0.  53  ml/sec,  G = 0.  5 y lO'^ 

(c)  Q = 1.10  ml/sec,  G = 1.  1 y 10‘^ 

(d)  Q = 1 . 75  ml/sec,  G=1.7yl0‘^ 

(e)  Q = 3. 00  ml/sec,  G = 2.  9 y lO'^ 

Fig.  48  Effect  of  injection  of  Polyox  (WSR  301)  on  the  laminar  separ- 
ation on  the  NSRDC  body.  Flow  is  from  right  to  left. 

5 

Reo  = 1.6  X 10  . Concentration  = 500  wppm. 

(a)  Q (injection  rate)  = 0.0  ml/sec,  G = 0.  0 

(b)  Q = 0.  1 ml/sec,  G = 0.  5 x 10  ^ 

(c)  Q = 0.3  ml/sec,  G = 1 . 5 x 10  ^ 
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Fig.  49 

Fig,  50 
Fig.  51 
Fig.  52 
Fig.  53 
Fig.  54 
Fig.  55 
Fig.  56 

Fig.  57 
Fig.  58 


(d)Q  - 0.5  ml/scc,  (<  2.5  y 10 

Fffpcts  of  injection  of  Polyox  (WSK  iOl)  on  (In-  laminar  sop 

aration  on  the  NSR  DC  body,  r'U)w  is  from  right  to  left. 

5 

Re,)  = 3.4  y 10  Concentration  • 500  wp|>m. 

(a)  Q (injection  rate)  -•  0.0  nil /sec,  C 0.0 

(b)  Q =0.1  ml /sec,  G = 0.36  y 10*^' 

(c)  Q = 0.3  ml /sec,  G = 1.08  v 10  *■’ 

Effect  of  injection  of  Polyox  (WSR  301)  on  the  length  of  the 
laminar  separation  on  the  hemisphere  nose  body. 

Effect  of  injection  of  Polyox  (WSR  301)  on  the  length  of  the 
laminar  separation  on  the  NSR  DC  body. 

Effect  of  injection  of  Polyox  (WSR  301)  on  the  height  of  the 
laminar  separation  on  the  hemisphere  nose  body. 

Effect  of  the  injection  of  Polyox  (WSR  301)  on  the  height  of 
the  laminar  separation  on  the  NSRDC  body. 

Effect  of  injection  of  Polyox  (WSR  301)  on  the  position  of  separ- 
ation on  both  the  hemisphere  nose  and  the  NSRDC  bodies. 

The  estimated  additional  boundary  layer  disturbance  amplifi- 
cation caused  by  the  presence  of  Polyox  (WSR  301). 
Photographs  of  the  physical  appearance  of  cavitation  on  the 
NSRDC  body  at  two  turbulence  levels.  Flow  is  from  right  to 
left.  Rep  = 3.4  y 10^. 

(a)  u' /U  = 0.  05%,  0=0.44  (b)  u ' /U  = 0.  65%,  o = 0.  36 

(c)  u'/U  = 0.  65%,  0 = 0.36  (d)  u'/U  = 0.65%,  o = 0.35 

The  effect  of  freestream  turbulence  level  upon  the  cavitation 
inception  number  on  the  NSRDC  body. 

Photographs  of  the  physical  appearance  of  cavitation  on  the 
hemisphere  nose  body  with  polymer  injection.  Flow  is  from 
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Fig.  59 


Fig.  60 
Fig.  61 
Fig.  62 

Fig.  63 


left  to  right.  Rcj,  = 6. 7 x 10^  , o 0.  59 
Concentration  = 500  wppni . 

(a)  Q (injection  rate)  - 0.0  ml/sor,  C 0.0 

(b)  Q = 0.  51  ml/sec,  G = 1.8«yl0‘^' 

(c)  Q = 1.  12  ml/sec.  G = 4.  14  x lo'^ 

(d)  Q = 1.42  ml/sec,  G-  5.24x10"^ 

Photographs  of  the  physical  appearance  of  cavitation  on  the 
NSRDC  body  with  polymer  injection.  Flow  is  from  right  to 
left.  Rep  = 3.4  X 10^. 

(a)  Q (injection  rate)  = 0.0  ml/sec.  , G = 0.  0,  0 = 0.  44 

(b)  (i  = 2.40  ml/sec.  , c (concentration)  = 20  wwpm, 

G = 0. 34  X 10'^,  o =0.45 

(c)  Q = 2.40  ml /sec.  , c = 20  wppm,  G = 0.34  x 10  ^ , 

O = 0. 34 

(d)  Q = 2.  00  ml/aec.  , c = 500  wppm,  G = 7,n  xlO  * . 
o = 0.  34 

The  effect  of  injection  of  Polyox  (WSR  301)  upon  cavitation 
inception  number  on  the  hemisphere  nose  body. 

The  effect  of  injection  of  Polyox  (WSR  301)  on  the  cavitation 
inception  number  on  the  NSRDC  body. 

The  effect  of  injection  of  Polyox  (WSR  301)  on  the  cavitation 
inception  number  on  the  hemisphere  nose  body  as  a function 
of  injection  rate  for  two  concentrations.  Rep  = 7.  5 y 10^. 

The  same  data  as  presented  in  Fig.  62  but  now  plotted  versus 
the  parameter  G to  show  the  correlation  of  the  cavitation  sup 
pression  with  the  amount  of  polymer  injected  into  the  boundarv 
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layer. 

Fig.  64  Comparison  of  the  cavitation  suppression  results  of  thi;  pre- 
sent study  with  those  of  several  other  investigations. 

Fig.  65  Comparison  of  nuclei  number  density  distribution  functions 
derived  from  the  experimental  results  of  several  investiga- 
tions . 

Fig.  66  Schlieren  photographs  of  the  effect  of  increasing  numbers  of 
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Chapter  I 
INTRODUCTION 

I.  1 Background 

The  capacity  to  withstand  tensile  forces  is  a property  normally 
associated  with  solids  only.  Liquids,  however,  are  also  capable  of 
sustaining  tensile  stresses  and  theoretical  estimates  for  a pure,  homo- 
geneous liquid  predict  ultimate  strengths  of  thousands  of  atmospneres  [l]  . 
Yet,  the  most  careful  experiments  yield  strengths  of  only  a few  hundred 
atmospheres  [l].  To  explain  this  discrepancy  it  is  postulated  that 
weak  spots  or  "nuclei"  exist  in  the  liquid  samples.  The  study  of  the 
"breaking"  of  the  liquid  (the  formation  of  a hole  or  cavity)  at  a nucleus 
site  and  the  subsequent  growth  and  eventual  collapse  of  the  cavity  is 
called  "cavitation."  If  the  resulting  cavity  is  filled  predominantly  with 
gases  formerly  dissolved  in  the  liquid,  the  process  is  called  "gaseous" 
cavitation. 

Liquids  are  usually  "broken"  either  by  heating  at  constant  pres- 
sure (boiling)  or  by  reducing  the  pressure  at  constant  temperature.  It 
is  the  latter  case  which  is  of  concern  to  hydraulic  engineers.  Hydrau- 
lic machines  (such  as  propellers,  pumps,  and  turbines),  valves,  hydro- 
foils, and  any  hydraulic  device  in  which  the  dynamically  reduced  pres- 
sure falls  to  a "low  enough"  value  are  susceptible  to  cavitation.  The 
presence  of  cavitation  (unless  specifically  designed  for)  causes  a loss 

* Numbers  in  brackets  refer  to  References  at  the  end  of  the  paper. 
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of  performance,  erosion  damaf^e  and  noise.  Thus  it  becomes  desir- 
able to  be  able  to  predict  the  set  of  operating  conditions  for  a parti  c.dar 
device  which  form  the  boundary  between  cavitalini;  and  nun  - lav  j t at  i nn 
regions . 

The  impracticality  of  observing  a full  scale  device  in  operation 
has  led  to  the  development  of  facilities  for  model  testing.  The  bound- 
ary between  the  cavitating  and  non- cavitating  regions  can  then  be  de- 
termined for  the  model.  To  transfer  these  model  results  to  a proto- 
type operating  at  dynam.ically  similar  conditions  requires  knowledge  of 
appropriate  similarity  parameters  and  scaling  laws.  A great  deal  of 
theoretical  and  experimental  work  has  been  done  to  determine  what 
these  significant  parameters  and  their  scaling  laws  are.  However,  at 
the  present  time  the  process  of  cavitation  onset  (or  inception)  and  the 
factors  controlling  the  process  are  not  fully  understood  [2,3]. 

The  main  parameter  used  to  describe  a cavitating  condition  is 
the  cavitation  number,  o.  It  is  defined  as: 


a 


P -P^ 

00  V 

1 IT  2 

2 P^oo 


where  p and  are  the  reference  pressure  and  velocity  respectively ; 

0 and  p are  respectively  the  liquid  density  and  vapor  pressure  at  its 

bulk  temperature.  In  water  tunnel  testing,  p^  and  are  the  test 

section  static  pressure  and  velocity  measured  at  some  distance  from 
the  body.  Standard  procedure  at  most  water  tunnel  facilities  during 
cavitation  tests  is  to  fix  the  water  velocity  and  slowly  reduce  the 
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pressure  until  cavitation  occurs.  The  value  of  o at  which  c.ivitation 
first  occurs  is  called  the  "incipient"  cavitation  number  and  i .s  demoted 
a..  By  raising  the  tunnel  pressure  the  cavitation  can  be  made  to  dis- 
appear. The  value  of  O at  which  it  first  disappears  is  called  the 
"desinent"  cavitation  number  and  is  denoted  In  almost  all  cases 

a hystersis  exists  and  O. . Now'  if  a was  the  only  significant  pa- 

rameter involved  in  the  inception  process,  then  a geometrically  sind- 
lar  model  operating  at  a dynamically  similar  condition  should  have  the 
same  value  of  a at  inception.  The  results  presented  in  Fig.  1 show 
this  is  not  the  case  and  it  can  be  seen  in  this  figure  that  the  value  of 
j.  depends  on  body  diameter  and  freestream  velocity.  Thus  n is  not 
the  only  pertinent  parameter  and  a complete  set  of  governing  parame- 
ters is  yet  to  be  found  [z].  The  dependence  of  a on  such  variables 
as  body  size,  freestream  velocity,  and  liquid  properties  is  called 
"scale  effect.  " 

Research  on  the  mechanism  of  cavitation  inception  and  scale  ef- 
fects is  conducted  in  many  facilities  worldwide.  Each  of  these  facil- 
ities has  its  characteristic  environment  of  the  test  section.  For  ex- 
ample, some  water  tunnels  have  resorbers  which  put  back  into  solu- 
tion the  gases  brought  out  by  the  cavitation.  These  facilities  usually 
have  very  few  freestream  bubbles  present  in  the  test  section,  while 
water  tunnels  without  resorbers  typically  have  a very  bubbly  test  sec- 
tion environment.  To  compare  cavitation  inception  in  these  various 
facilities  the  International  Towing  Tank  Committee  (ITTC)  organized 
a comparative  test  [4].  Each  facility,  using  its  own  standard  proce- 
dure, measured  cavitation  inception  on  a model  of  prescribed  geome- 
try. Some  results  of  these  tests  are  presented  in  Fig.  Z.  Neithe r the  quan- 
titative nor  the  qualitative  agreement  between  the  facilities  is  acceptable  . 
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Associated  with  the  large  variance  in  the  values  of  o.  is  an 

1 

equally  disconcerting  variation  of  the  physical  appearance  of  the  cavi- 
tation. Figure  3 shows  a number  of  photographs  of  cavitation  incep- 
tion on  the  ITTC  standard  headform  body  taken  in  different  facilities  [5], 
Photographs  4,5,9  show  "band"  type  cavitation  inception.  It  usually 
occurs  suddenly  without  any  warning,  precurser  bubbles  forming  an 
attached  "sheet”  cavity.  Photograph  6 is  an  example  of  "spot"  type 

cavitation  inception.  It  starts  at  a fixed  spot  near  Cq  and  has  a 

^min 

region  of  V-shaped  attached  cavitation  behind  it  [3].  Photographs  1,3, 
7,8  show  "traveling  bubble"  type  inception.  These  three  types  of  in- 
cipient cavitation  are  the  ones  that  typically  occur  on  smooth  two-dimen- 
sional bodies  and  can  occur  singly  or  in  any  combination  [3].  It  should 
be  expected  that  the  cavitation  indices  of  these  various  forms  will  dif- 
fer greatly  and  the  type  that  does  occur  depiends  upon  the  liquid  envi- 
ronment. It  is  then  obvious  that  it  is  necessary  to  document  the  liquid 
environment  during  cavitation  inception  experiments  — a fact  recog- 
nized but  not  extensively  practiced  at  the  present  time. 

Scale  effects  and  liquid  environment  effects  by  themselves  are 
the  source  of  formidable  experimental  problems.  Yet,  recently  an- 
other consideration  has  been  added  by  the  introduction  of  small  amounts 
of  drag- reducing  polymers  into  the  liquid.  This  has  been  motivated  by 
the  discovery  that  certain  naturally  occurring  organisms  (such  as  al- 
gae) secrete  drag- reducing  chemicals  into  the  environment  [6]  and 
also  by  the  intentional  introduction  of  polymers  into  the  fluid  to  reduce 
viscous  drag  on  a vehicle.  Thus  the  propulsive  system  of  a vehicle 
(usually  a propeller,  but  in  the  case  of  a jet  boat  a pump)  may  operate 
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in  a drag- reducing  solution.  Tests  to  deterinine  the  cavitation  char- 
acteristics of  the  propulsive  systeni  in  a d rag  - reducing  er  vi  ronnwnt 
must  then  be  carried  out.  However  before  discussing  the  preceding 
topics  in  detail  it  is  a good  idea  to  digress  slightly  and  consider  some 
physical  aspects  of  cavitation  inception. 


I.  2 Cavitation  Inception  — .Some  Physical  Considerations 


Previously  it  was  mentioned  that  cavitation  started  at  "weak 
spots"  or  nuclei  in  the  liquid.  And  since,  in  essentially  all  cases  of 
engineering  concern,  ca'.itation  would  not  occur  if  it  were  not  for  the 
presence  of  cavitation  nuclei,  it  is  worthwhile  giving  them  a great  deal 
of  consideration.  There  are  tnree  questions  of  interest  about  nuclei: 
ncimely,  what  is  their  physical  nature,  their  origin,  and  their  role  in 
the  inception  process.  First,  what  are  the  nuclei?  Since  it  is  an  ob- 
vious candidate  for  the  physical  m.odel  of  a nucleus,  the  free  gas  bub- 
ble has  received  the  most  attention.  However,  Epstein  and  Plesset 
[7]  showed  that  theory  predicts  free  gas  bubbles  will  disappear  from 
an  under  saturated  liquid  either  by  dissolving  or  by  rising  to  the  high 
points  of  the  system.  Yet,  in  the  most  careful  experiments  where  a 
water  sample  has  been  carefully  treated  to  remove  any  free  gas  bub- 
bles, the  sample  still  "breaks"  or  cavitates  at  tensions  much  less 
than  theoretically  predicted  [l].  To  explain  the  evident  persistence  of 
the  nuclei,  several  additional  models  have  been  proposed. 

Since  these  models  have  been  discussed  extensively  in  the  lit- 
erature, the  comments  here  will  be  limited  to  a brief  description  of 
the  model  and  its  current  rating  among  cavitation  investigators.  Fox 
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and  Herzfeld  [h]  su^geslod  that  an  t)ri>anio  skin  protects  the  gas  hutj- 
ble  by  preventing  diffusion.  Peas<’  and  Blinks  [ | discardt  d the  ga.s 

bubble  completely  and  postulated  that  the  nucleus  consi.st.s  of  a solid, 
hydrophobic  particle.  Both  these  rnod»-ls  have  bet  ri  considered  by 
Plesset  [ij  and  Holl[lO]  anci  for  various  leaso  is  they  iiave  givi  n un- 
favorable opinions  as  to  their  sig.’iiticance . The  niodel  presently  con- 
sidered with  the  most  favor  is  that  proposed  by  Harvey  et  al  111]. 
Harvey  suggested  tliat  undissolved  uas  could  exist  in  pockets  or  pores 
of  microscopic,  hydrophobic  crevices  in  container  walls  or  in  freely 
suspended  solid  particles.  These  gas  pockets  would  then  be  the  nuclei. 

I Experiments  by  Keller  [l2,  13]  and  Knapp  [H]  have  provided  some 

indirect  evidence  supporting  this  model.  However,  the  absolute  evi- 
dence as  to  what  the  nuclei  actually  are  that  could  be  provided  by  ili- 
rect  observation  is  unavailable.  This  is  so  because  of  the  extreme' 

I difficulties  encountered  in  trying  to  observe  the  very  sniall  nuclei 

(usually  less  than  100  microns  in  diameter)  as  they  are  moving  at  high 

i 

i velocities  (normally  greater  than  20  feet  per  second).  The  only  pres- 

t ently  available  observations  of  cavitation  nuclei  under  test  conditions 

. are  those  made  by  pulsed  laser  holography  [ 1 5,  16,  17]  . In  these  in- 

vestigations it  was  found  that  both  gas  bubbles  and  solid  particles 
were  present.  However,  the  resolution  was  insufficient  to  determine 
if  the  solid  particles  had  any  pockets  of  trapped  gas.  Thus,  at  the 
present  time  a complete  description  of  the  cavitation  nucleus  is  still 
not  available. 

The  second  question  that  we  asked  about  nuclei  was:  what  are 
their  origins?  This  can  be  answered  easily  (if  somewhat  superficially) 
for  experimental  facilities.  There  are  two  sources  of  cavitation  nuclei 
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in  a water  tunnel  — the  incominjii  flow  and  the  Lest  tjody  itself.  Ihese 
are  referred  to  as  "freestrearn"  and  "surfar-e"  nuclei  respe<  tively  . 

The  third  question  was:  what  roles  do  these  surface  and  treestreani 

nuclei  play  in  the  inception  process?  First  consider  surtact'  nuclei 
and  then  freestrearn  nuclei. 

I.  3 Surface  Nuclei  and  Inception 

Assuming  that  surface  characteristics  control  surface  nuclei 
behavior,  several  investigations  have  been  conducted  on  the  influences 
of  surface  nuclei  in  cavitation  inception  by  controlling  the  test  body 
surface  quality.  In  these  investigations  roughness  effects  were  not  a 
consideration.  Acosta  and  Hamaguchi  [l8]  dipped  the  nose  of  the 
ITTC  test  body  into  a mixture  of  silicon  oil  and  solvent.  Since  silicon 
oils  have  the  ability  to  dissolve  large  amounts  of  air,  the  effect  of  the 
dipping  was  to  create  a layer  of  potential  nuclei  on  the  surface  of  the 
model.  Ensuing  inception  tests  showed  the  incipient  cavitation  indices 
to  be  systemically  higher  than  for  the  pure  water  case.  But  as  the 
coating  aged  and  became  deaerated  by  the  evolution  of  cavitation  bub- 
bles, the  values  of  the  cavitation  index  fell  below  those  for  pure  water 
and  the  appearance  of  the  cavitation  changed.  This  latter  result  was 
never  resolved  by  the  investigators. 

Ho  11  and  Treaster  [ 19],  Holl  [ZO],  Reed  [Zl],  and  Gupta  [ZZ] 
studied  cavitation  inception  on  hemisphere  nose  bodies  having  surface 
finishes  of  smooth  stainless  steel,  sandblasted  stainless  steel,  waxed 
("simonized")  stainless  steel,  teflon,  nylon,  and  glass.  They  found 
that  whereas  teflon  could  be  made  to  cavitate  easily,  glass  could  be 
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made  to  cavitate  only  with  difficulty  and  the  siinonized  stainless  sli  » l 
was  harder  to  cavitate  than  the  smooth  stainless  steel  suiface.  'I  hcsi 
observations  led  Holl  [Z]  to  saf.’gest  that  surlacc-  porosily  is  the  con- 
trolling factor. 

V'an  der  Meulen  [23,24]  also  studied  inception  and  desinence  on 
stainless  steel  and  teflon  hemisphere  nose  bodies.  He  found,  like  Holl 
[20],  that  the  teflon  body  was  easier  to  cavitate  than  the  stainless  steel 
one  and  that  the  cavitation  appearance  was  different  for  the  different 
materials.  He  suggested  tliat  whereas  surface  nuclei  were  responsible 
for  cavitation  on  the  teflon  body,  freestream  nuclei  were  responsible 
for  cavitation  on  the  stainless  steel  model. 

Further  tests  of  this  type  were  carried  out  by  Peterson  [25]  on 
a modified  ellipsoidal( NSRDC  )body( see  Fig.  4)  which  was  first  cleaned  with 
acetone  and  then  air  dried.  When  tested,  it  produced  a cavitation  index 
of  0.780  at  20  feet  per  second.  The  headforrn  was  then  soaked  in  ace- 
tone for  65  hours,  washed  in  distilled  water,  quickly  pressurized  to 
19,000  psi  and  installed  in  the  filled  water  tunnel.  All  steps  were  car- 
ried out  so  that  the  model  was  never  exposed  to  the  atmosphere.  The 
average  incipient  cavitation  index  after  this  treatment  was  0.625  at 
20  feet  per  second  and  was  found  to  be  the  same  five  days  later  after 
continuous  submergence.  The  headforrn  was  then  air  dried  and  retested. 
The  value  of  the  cavitation  index  then  returned  to  its  initial  value  of  0 . 780 . 

As  a summary  of  the  above  review,  it  can  be  concluded  that  not 
only  the  surface  characteristics  but  also  the  pressure -time  history 
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expe  rienced  by  the  surfacearc  important  in  dote  r mining  the  occur  re  nee  of 
cavitation.  These  observations  , though,  are  for  cases  where  surface  nuclei 


control  inception.  What  about  the  cases  where  freestream  nuclei  con- 
trol inception^ 

1.4  Freestream  Nuclei  and  Inception 

Before  discussing  the  role  of  freestream  nuclei  in  cavitation,  let 
us  briefly  return  to  the  nature  of  these  nuclei.  In  Section  1.4  it  was 
pointed  out  that  a free  gas  bubble  should  be  eliminated  frori'  considera- 
tion as  a nucleus  since  it  will  either  dissolve  or  rise  to  the  high  points 
of  the  system  and  be  removed.  However,  this  statement  assumes  that 
the  gas-water  system  is  in  equilibrium.  This  is  hardly  the  case  in  a 
water  tunnel  under  test  conditions  where  free  gas  bubbles  are  being  gen- 
erated throughout  the  circuit  — especially  at  the  tunnel  pump  f 28]  . Al- 
so recall  that  the  holographic  observations  of  freestream  nuclei  men- 
tioned earlier  showed  both  free  gas  bubbles  and  solid  particles  to  be 
present.  Therefore,  the  free  gas  bubble  will  be  brought  back  into  con- 
tention as  a model  for  the  cavitation  nucleus.  And  as  we  shall  soon  see, 
until  very  recently  the  free  gas  bubble  was  the  major  model  used  as  the 
cavitation  nucleus  in  experimental  and  theoretical  investigations. 
Therefore  before  reviewing  these  investigations  it  is  helpful  to  briefly 
discuss  some  aspects  of  the  growth  of  a free  gas  bubble. 

An  analysis  of  the  static  stability  of  a stationary  spherical  gas 
bubble  containing  the  liquid  vapor  and  a small  amount  of  permanent  gas 
shows  (see  for  example  Knapp  et  al  [26]i  that  for  an  initial  bubble  radi- 
us Rq  there  exists  a critical  ambient  pressure  p^  below  which  the 
bubble  cannot  be  in  static  equilibrium.  If  the  ambient  pressure  falls 
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below  this  value,  the  bubble  will  grow  dynaniical  l>  umi«T  the  conditions 
of  constant  vapor  pressure  and  gas  content.  In  this  case  the  bul)ljle  is 
said  to  grow  by  "vaporous”  cavitation.  As  the  initial  size  of  the  bubble 
increases,  the  critical  pressure  p^  approaches  the  fluid  vapor  pres- 
sure asymptotically'.  Thus  in  order  for  vaporous  cavitation  to  occur, 
the  local  ambient  pressure  must  be  below  the  vapor  pressure.  And  if 
the  initial  bubble  is  very  small,  it  may  actually  be  necessary  to  impose 
a tension  on  the  liquid  before  the  bubble  becomes  unstable. 

As  can  be  seen  from,  the  i*xperimental  results  of  Fig.  2,  cavita- 
tion also  occurs  when  the  liquid  pressure  is  above  the  vapor  pressure. 
This  can't  be  vaporous  cavitation  and  it  has  been  suggested  by  Holl  [z] 
that  gaseous  diffusion  is  responsible  for  the  bubble  growth.  "Gaseous" 
cavitation,  as  this  process  is  called,  can  occur  anytime  the  ambient 
liquid  is  supersaturated  with  respect  to  the  partial  pressure  of  the  gas 
in  the  bubble.  Thus  bubbly  growth  by  diffusion  is  possible  even  when 
the  ambient  pressure  is  above  the  vapor  pressure.  Vaporous  cavitation 
is  a much  faster  phenomenon  compared  to  gaseous  cavitation.  Typical 
times  for  vaporous  growth  are  of  the  order  of  milliseconds  [ 27  j,  where- 
as gaseous  growth  may  be  the  order  of  a few  seconds  [7],  However,  to 
differentiate  between  vaporous  and  gaseous  cavitation  in  a test  situation 
is  often  difficult. 

Now  let  us  begin  the  discussion  of  freestream  nuclei  and  their 
role  in  cavitation  inception.  Ripken  and  Killen  [28]  studied  the  influence 
of  freestream  gas  bubbles  on  cavitation  inception.  Using  an  acoustic 
technique  to  measure  free  gas  volume  they  tried  to  relate  the  free  gas 
volume  to  the  occurrence  of  inception.  They  found  that  many  of  the  tests 
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did  not  correlate  well  with  free  gas  content  and  suggeb‘e<l  that  a meas- 
urement of  the  gas  bubble  size  and  number  dist  ril>ution  would  be  more 
meaningful.  They  farther  concluded  that  the  principal  biiblde  source  is 
the  pump  and  that  the  resulting  bubble  numbers  anti  sizes  at  the  test  sec- 
tion would  be  determined  by  the  time-pressure  characteristics  ol  the 
flow  circuit,  i.e.  each  facility  will  produce  its  own  characteristic  bub- 
ble spectrum.  Killen  and  Ripken  [Z9]  continued  work  on  the  acoustic 
free  gas  meter  to  give  it  the  capability  of  determining  bubble  sizes  and 
concentrations.  Schiebe  [30],  using  this  upgraded  version,  tried  to  es- 
timate the  number  of  cavitation  events  that  should  occur  for  a given 
body,  tunnel  conditions  and  measured  nuclei  distribution.  He  found  that 
the  predicted  number  of  events  compared  poorly  with  the  measured  num- 
ber of  events  and  claimed  that  insufficient  accuracy  of  the  acoustic 
method  was  the  probable  cause. 

The  problem  of  accurately  determining  nuclei  concentrations  and 
size  distributions  has  hampered  the  study  of  the  role  of  freestream  nu- 
clei in  the  inception  process.  Morgan  [80]  has  reviewed  the  existing 
nuclei  counting  techniques  and  has  proposed  certain  criteria  which  an 
ideal  counter  should  satisfy.  However,  most  of  the  present  techniques 
fail  to  meet  these  criteria  and  it  has  only  been  with  the  development  of 
laser  technology  that  more  promising  methods  have  become  available. 
One  of  these  new  techniques  has  been  used  by  Keller  [ IZ,  13]  to  investi- 
gate the  role  of  freestream  nuclei  in  cavitation  inception. 

Keller[12,  13,31, 3Z]has  investigated  the  effects  of  nuclei 
on  cavitation  inception  by  the  direct  measurement  of  nuclei  distributions 
^ using  an  optical  technique.  Since  Keller's  nuclei  measuring  system 
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was  a major  advance  in  this  type  of  work,  it  is  worthy  of  a brief  des- 
cription before  discussing  his  results.  Basically  the  sy  stem  dete  r mi ne s 
the  size  of  a nucleus  by  measuring  the  intensity  of  light  scattered  Ijy  the 
nucleus  as  it  passes  through  a laser  beam.  In  Keller's  configuration, 
a control  volume  was  defined  by  the  intersection  of  a small  diameter 
laser  beam  and  the  focal  volume  of  the  receiving  optics.  When  a nucleus 
passed  through  this  control  volume,  it  would  scatter  some  of  the  laser 
light  into  the  receiving  system.  In  the  receiving  system  the  intensity  of 
the  scattered  light  would  be  measured  and  from  a previous  calibration 
the  nucleus  size  could  be  determined.  Keller  connected  the  output  of  the 
receiving  system  into  a recorder  where  it  was  stored  until  the  data  col- 
lection was  finished.  Then  the  recorder  would  provide  a nuclei-size 
histog  ram. 

Using  this  counting  method,  Keller  [12]  measured  nuclei  popu- 
lations in  four  types  of  water  (tap  water,  degassed  tap  water,  filtered 
water,  degassed-f iltered  water)  and  then  tested  a half  body  (formed  by 
the  potential  flow  solution  of  a point  source  in  a uniform  flow)  in  each 
type.  The  results  of  these  tests  showed  that  by  changing  from  highly 
nucleated  tap  water  to  a low  nuclei  population  in  filtered-degassed  water 
the  cavitation  number  at  inception  could  be  reduced  by  as  much  as  a fac - 
tor  of  ten!  The  kind  of  cavitation  was  always  of  the  traveling  bubble 
type. 

Unfortunately,  Keller's  nuclei  counter  could  not  distinguish  be- 
tween free  gas  bubbles  and  particulates.  However,  some  observations 
of  the  effects  of  various  water  treatment  (degassing,  filtering,  aging)  on 
the  nuclei  populations  provided  some  useful  conclusions.  For  example. 
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Keller  found  that  whereas  an  "alteration  of  the  total  ^as  content  aflei  ted 
the  measured  nucleus  spectrum  only  moderately,  many  nuclei  wer»* 
withdrawn  from  the  water  by  filtering"  [13],  He  also  found,  that  if  the 
water  was  slowly  circulated  for  an  hour  before  the  cavitation  test,  the 
cavitation  index  could  be  reduced  by  a factor  of  two  in  tap  water  whereas 
it  would  be  reduced  by  a factor  of  ten  in  filtered  water.  Any  longer 
waiting  time  before  the  cavitation  test  showed  no  further  suppression  of 
cavitation  or  reduction  in  the  nuclei  population.  We  infer  from  tfiese 
results  that  although  free  gas  bubbles  and  particulates  can  both  act  as 
nucleation  sites,  the  particulates  are  more  important  (at  least  in  Keller's 
facility).  Keller  apparently  agrees  with  this  inference,  since  he  for- 
warded an  inception  theory  based  on  the  Harvey  "pore"  model  which 
will  be  discussed  in  a later  section. 

Keller's  work  showed  very  dramatically  that  the  freestream 
nuclei  population  is  a major  factor  influencing  cavitation  inception 
Thus,  any  similarity  law  cannot  be  useful  if  it  does  not  take  proper  ac- 
count of  the  nucleus  spectrums  of  the  prototype  water  and  the  teat  water 
However,  although  this  fact  is  generally  appreciated  in  the  field,  the 
lack  of  full  understanding  of  how  these  nuclei  initiate  cavitation  restricts 
one  to  passing  only  qualitative  judgements  as  to  the  quantitative  effect 
of  these  various  nuclei  populations  on  the  actual  cavitation  events. 

In  the  previous  two  sections  we  have  discussed  experimental 
investigations  of  the  role  of  the  cavitation  nucleus  in  the  inception 
process.  Theoretical  investigations  of  this  problem  have  also  been 
made.  And  although  the  role  of  cavitation  nuclei  is  still  not  fully 
understood,  it  is  worth  briefly  mentioning  these  theories  to  gain  insight 
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into  how  the  thinking  about  cavitation  inception  has  devc-loped. 

I.  5 Cavitation  Inception  — Theories 

From  the  earlier  analysis  of  the  stat^  stability  of  a spiierical 
gas  bubble,  it  can  be  shown  that  for  greater  than  250  microns,  it 

would  be  expected  that  they  would  grow  vaporously  as  soon  as  the  pres- 
sure in  the  neighborhood  of  the  body  fell  to  the  vapor  pressure.  And 
since  the  lowest  pressure  available  is  the  minimum  pressure  on  the 
body  it  has  been  customary  to  set  the  cavitation  inception  index  equal 
to  the  negative  of  the  minimum  pressure  coefficient.  A quick  glance 
at  Figs.  1 and  2 show  this  to  be  a poor  correlation. 

An  early  attempt  to  more  rigorously  model  the  cavitation  incep- 
tion process  theoretically  was  made  by  Parkin  in  1952  [33].  Parkin 
assumed  that  cavitation  was  initiated  from  small  solid  particles  with 
stabilized  pockets  of  air  and/or  vapor.  Parkin  recognized  that  a spec- 
trum of  nuclei  sizes  probably  existed;  but  having  no  experimental  or 
theoretical  evidence  as  to  what  the  spectrum  would  be,  he  assumed  a 
constant  size.  These  nuclei  were  then  assumed  to  travel  along  the 
body  surface  at  the  local  freestream  velocity  and  the  nuclei  growth  was 
estimated  using  the  Rayleigh-Plesset  [34]  equation  neglecting  viscous 
and  thermal  terms.  Inception  was  said  to  have  occurred  when  the  nu- 
clei had  grown  to  a 1 millimeter  radius.  His  predicted  values  of  the 
inception  index  showed  qualitative,  but  poor  quantitative  agreement 
with  experiment  and  Parkin  attributed  the  lack  of  agreement  to  the  ne- 


glect of  the  influence  of  the  boundary  layer. 

To  investigate  what  the  effects  of  the  boundary  layer  could  be. 
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Parkin  and  Kermeen  [35]  carried  out  further  experinienta  on  a hi  toi- 
sphere  nose  body.  In  what  is  now  considered  a i:lassic  experiment 
tliey  showed  that  the  boundary  layer  has  a significant  influence  on  cavi- 
tation and  from  their  observations  they  concluded  the  following  events 
take  place  at  inception.  Inception  would  first  occur  as  a ring  of  mac- 
scopic  bubbles  around  the  body.  These  macroscopic  bubbles  are  sus- 
tained by  very  small  microscopic  bubbles  which  grow  in  the  boundary 
layer  upstream  of  the  macroscopic  cavitation.  These  microscopic 
bubbles  would  grow  in  time  at  a fixed  position  until  the  bubble  radius 
reached  a value  equal  to  one  half  of  the  boundary  layer  displacement 
thickness.  Then  the  microscopic  bubble  would  be  entrained  by  the 
flow  and  travel  downstream  to  be  distorted  into  irregular  patches 
which  sustained  the  macroscopic  cavitation.  A reduction  in  the  tunnel 
pressure  from  this  point  would  increase  the  rate  of  microscopic  bubble 
evolution  and  cause  more  profuse  macroscopic  cavitation  until  a clear 
attached  cavity  was  formed.  These  interpretations  made  by  Parkin 
and  Kermeen  of  their  observations  are  very  important  because  they 
were  incorporated  into  many  of  the  subsequent  theories  of  cavitation 
as  we  shall  see. 

Roll  and  Kornhauser  [36]  carried  out  an  analysis  similar  to 
Parkin's  [33]  but  included  the  effect  of  initial  conditions,  stream  versus 
surface  nuclei  and  thermodynamic  effects  using  the  complete  Rayleigh- 
Plesset  equation.  Boundary  layer  effects  were  taken  into  account  by 
assuming  that  the  bubble  had  a velocity  equal  to  that  in  the  boundary 
layer  at  one  nuclei  radius  away  from  the  body  surface.  Again,  incep- 
tion was  said  to  occur  when  the  bubble  diameter  reached  an  arbitrary 
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diameter  — in  this  case  250  microns.  Comparison  of  theory  with  ex- 
perimental results  for  a hemisphere  nose  body  showed  qualitative 
agreement  for  diameter  and  size  effects.  But  tlieory  predicted  tiie  in- 
cipient cavitation  number  should  dei  rease  with  teni[)e  nature , when  as 
experiment  showed  it  increased.  The  discrepancy  was  attributed  to 
convective  heat  transfer  effects  which  were  neglected  in  the  cali;ula- 
tions. 

Another  theory  was  proposed  by  van  der  Walle  [37]  who,  moti- 
vated by  Parkin  and  Kermeen's  observations,  presented  a theory  which 
is  concerned  mainly  with,  the  stabilization  of  bubbles  in  the  boundary 
layer  under  the  action  of  drag,  pressure  gradient,  and  adhesion  forces 
and  the  resulting  importance  of  gaseous  cavitation.  Taking  the  above 
factors  into  account  and  doing  calculations  similar  to  Parkin,  he  ob- 
tained cavitation  inception  values  for  bodies  with  different  pressure 
gradients.  Here  too,  comparison  with  experimental  results  for  a hemi- 
sphere nose  body  showed  only  qualitative  agreement. 

At  about  the  same  time  van  der  Walle  presented  his  theory, 
Oshima  [38]  attempted  to  formulate  scaling  laws  based  on  similarity 
ideas.  He  assumed  that  bubbles  undergo  dynamically  similar  growth  as 
specified  by  the  Rayleigh-Plesset  equation  and  that  this  growth  occurs 
at  geometrically  similar  locations  for  geometrically  similar  bodies. 

He  also  assumed  that  bubbles  move  with  the  liquid  at  a distance  from 
the  wall  equal  to  a constant  percentage  of  the  boundary  layer  thickness 
and  inception  was  assumed  to  take  place  when  the  bubble  diameter 
equalled  the  local  boundary  layer  displacement  thickness.  Oshima's 
results  agreed  reasonably  well  with  experiments  for  hemisphere  nose 
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bodies  with  small  diameters  but  departed  from  the  ocpe rimeiital  values 
as  the  diameter  increased. 

The  above  theories  an-  basi'd  on  inte rprt-talions  of  Parkin  and 
Kermeen's  observations  and  thus  boundary  layiT  considerations  play  a 
main  part.  Johnson  and  Hsich  [19]  introduced  a theory  in  which  bound- 
ary layer  considerations  are  neglected.  They  assumed  that  freestream 
gas  bubbles  convected  into  the  minimum  pressure  region  of  the  body  are 
responsible  for  cavitation  inception.  They  then  calculated  bubble  tra- 
jectories and  examined  the  bubble  stability  along  the  trajectory  for  dif- 
ferent combinations  of  body  size,  freestream  velocity  and  initial  nuclei 
sizes.  The  results  of  the  calculations  showed  a "screening"  effect 
whereby  the  pressure  gradient  force  deflected  the  large  bubbles  away 
from  the  minimum  pressure  region.  This  "screening"  effect  was  found 
to  be  body  size  and  tunnel  velocity  dependent.  Thus  for  a given  body 
size,  tunnel  velocity  and  tunnel  pressure^a  range  of  "cavitatable"  bub- 
bles was  found.  Below  this  range  the  bubbles  were  too  small  to  become 
unstable  whereas  the  bubbles  above  this  range  would  be  deflected  away. 
As  the  tunnel  pressure  increased,  this  range  would  decrease  in  size 
until  there  would  be  only  one  unstable  bubble  size.  Therefore,  accord- 
ing to  this  theory,  for  a particular  velocity  and  body  a unique  entrained 
bubble  size  should  exist  whose  stability  defines  inception. 

Schiebe  [30]  used  this  theory  in  conjunction  with  some  actual  nu- 
clei distributions  obtained  by  using  the  acoustic  method.  He  estimated 
the  number  of  cavitation  events  that  would  occur  for  a given  body  size, 
velocity  and  cavitation  number.  However,  comparison  with  the  actual 
number  of  events  as  determined  by  a hydrophone  located  in  the  body 


showed  only  qualitative  af>recmenl. 

Most  recently,  Keller  [13]  has  sugi^estird  a theory  liased  cmi  th«- 
Harvey  nucleus.  Keller  estimated  the  growth  of  a pocket  of  gas  trappc-d 
in  a crevice  of  a hydrophobic  particle  as  it  traveled  along  the  body  sur- 
face with  the  local  freestream  velocity.  This  was  accomplished  by  con- 
sidering the  growth  of  a volume  equivalent  spherical  gas  bubble  and 
using  the  Rayleigh-Plesset  equation.  Inception  was  said  to  occur  when 
the  radius  of  the  phase  boundary  arching  outside  the  crevice  equalled 
the  crevice  radius.  He  then  compared  his  critical  pori-  size  with  his 
measured  nuclei  spectrum  and  found  it  was  approximately  equal  to  the 
.argest  measured  nuclei.  Assuming  then  that  the  largest  nuclei  cavi- 
tates first,  Keller  said  he  had  good  agreement  between  his  predicted 
inception  and  experiment.  However,  the  concept  did  not  work  well  for 
the  blunter  bodies  tested,  and  Keller  attributed  this  to  viscous  effects 
inside  the  crevice  walls  which  were  not  included  in  his  analysis. 

I.  6 Summary  of  Cavitation  Inception  Theories 

In  general,  the  prediction  theories  reviewed  above  calculate 
the  growth  of  a spherical  bubble  from  an  initial  radius  to  a pre- 

scribed final  radius  as  it  experiences  a pressure -time  history  deter- 
mined by  the  tunnel  velocity  and  test  model.  The  bubble  growth  is  es- 
timated by  the  Rayleigh-Plesset  equation  and  it  is  assumed  that  the  be- 
havior of  a single  particle  characterizes  the  problem.  Attempts  to  in- 
clude viscous  effects  are  based  on  Parkin  and  Kermeen's  [35]  inter- 
pretation of  their  observations  made  by  high-speed  photography.  How- 
ever, the  theories  used  this  interpretation  without  any  knowledge  of  the 
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actual  viscous  flow  over  the  body.  For  examplt',  on  the  hemisphere 
nose  body,  Holl  and  Kornhauser  [3b]  assumed  a laminar  to  turbulenl 
boundary  layer  transition  occurred  at  the  minimum  pressure  poinl; 
Oshima  also  assumed  a laminar  to  turbulent  transition  oceurreel.  In 
the  next  section  it  will  be  shown  that  not  only  is  this  assumption  wront’ 
but  also  that  Parkin  and  Kermeen's  interpretation  of  their  observations 
was  incornect.  Thus,  it  is  not  too  surprising  that  the  quantitative  re- 
sults of  these  theories  do  not  agree  well  with  experiment. 

I.  7 Viscous  Effects  and  Cavitation  Inception 

Even  though  Parkin  and  Kermeen  urged  that  further  work  on 
boundary  layer  effects  be  done,  it  was  not  until  very  recently  that  a 
careful  investigation  of  viscous  effects  was  undertaken.  Arakeri 
and  Acosta  [40,4l]  used  the  schlieren  flow  visualization  technique  to 
observe  the  fully  wetted  flow  and  the  inception  of  cavitation  on  two 
axisymmetric  bodies  (a  hemisphere  nose  body  and  an  ITTC  body). 
Arakeri  observed  that  Holl  and  Kornhauser's  [36]  and  Oshima's  [38] 
assumption  of  a laminar  to  turbulent  transition  was  wrong  and  that  in- 
stead a laminar  separation  existed  on  the  two  bodies  at  all  test 
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Reynolds  numbers  (3.8  x 10  to  9.  1 X 10  for  the  hemisphere  nose 
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body  and  3 . 7 X 10  to  8.8x10  for  the  ITTC  body).  Arakeri  also 
determined  that  Parkin  and  Kermeen's  interpretation  of  their  obser- 
vations was  incorrect  and  that  the  actual  sequence  of  events  in  the  in- 
ception process  on  the  hemisphere  nose  body  is  as  follows:  Inception 
begins  in  the  reattachment  zone  of  the  laminar  separation  and  appears 
as  a ring  of  microscopic  and  macroscopic  bubbles  around  the  body. 
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Some  of  the  microscopic  bubbles  were  entrained  in  the  reverse  flcjw 
in  the  separated  region  and  carried  to  the  separation  point  where  they 
would  then  be  entrained  by  the  flow  external  to  the  separated  region 
and  carried  downstream  within  the  free  shear  layer  back  to  the  posi- 
tion of  the  macroscopic  cavitation.  A few  of  the  microscopic  bubbles 
were  observed  to  grow  at  a fixed  point  just  downstream  of  the  separa- 
tion. When  these  bubbles  grew  to  the  height  of  the  separated  region, 
they  were  entrained  into  the  main  flow  and  carried  downstream  — as 
was  observed  but  misinterpreted  by  Parkin  and  Kermeen.  If  the  tun- 
nel pressure  was  then  lowered  slightly,  the  cavitation  would  become 
more  profuse  and  a clear,  attached  cavity  was  formed.  This  sequence 
of  events  from  the  macroscopic  bubble  cavitation  to  "band"  type  cavita- 
tion and  its  relation  to  the  laminar  separation  are  clearly  demonstrated 
in  the  schlieren  photographs  of  Fig.  5. 

The  authors  speculated  that  large  pressure  fluctuations  in  the  re- 
attachment were  responsible  for  cavitation  inception  occurring  at  that 
position.  Arakeri  [42]  then  measured  the  pressure  fluctuations  in  this 
region  and  obtained  values  of  peak  fluctuations  of  25  percent  of  the  dy- 
namic head  and  thus  supports  their  speculation.  Further,  the  authors 
found  a good  correlation  between  the  cavitation  inception  index  and  the 

negative  of  the  pressure  coefficient  at  the  position  of  separation  (-  c ). 

^s 

However,  for  the  hemisphere  nose  body  it  was  found  that  a.  > - c 

^ ^s 

whereas  for  the  ITTC  body  o.  < - c . The  authors  had  no  explanation 

^ ^s 

for  this  difference  in  behavior,  but  speculated  that  the  height  of  the 
separated  region  might  be  responsible  (the  height  of  the  separated  re- 
gion on  the  hemisphere  nose  body  is  two  to  three  times  greater  than  on 
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the  ITTC  body). 

Arakeri  and  Acosta[44jcontinued  their  investigation  of  viscous  < f- 
fects  on  bodies  with  a laminar  separation  by  next  consid(?ring  the  effei  1 of 
removing  the  separation.  They  investigated  this  topic  fur  three  axisym  - 
metric  bodies:  a hemisphere  nose , an  I'J'TC , and  a modified  i 1 lipsoidai 
body  provided  by  Peterson.  The  laminar  separation  was  eliminated  Ijy 
stimulating  the  boundary  layer  so  that  an  early  transition  would  occur. 

This  was  done  artificially  on  the  hemisphere  nose  and  ITTC  bodies  by  the 
addition  of  a trip  ring , whereas  a natural  transition  on  the  modified  ellip- 
soidal body  was  obtained  by  increasing  the  freestream  velocity.  The  effects 
of  the  trips  on  the  hemisphere  and  ITTC  bodies  was  then  to  lower  the  c ritical 
Reynolds  number  to  a value  easily  obtainable  in  the  HS  WT.  The  critical 
Reynolds  niomber  is  that  value  of  the  Reynolds  number  at  which  the  positions 
of  laminar  separation  and  transition  coincide . 

In  the  subsequent  cavitation  tests  on  the  tripped  hemisphere  nose 
body,  the  following  observations  were  made.  At  Reynolds  numbers  less 
than  the  critical  value,  band  type  inception  took  place  indicating  the 
presence  of  a laminar  separation.  As  the  Reynolds  number  approached 
the  critical  value  from  below,  the  cavity  surface  changed  from  smooth 
to  very  "rough"  and  the  cavity  separation  changed  from  a well  defined 
line  to  a very  "wavy"  line.  If  the  Reynolds  number  was  increased  fur- 
ther until  it  was  near  the  critical  value,  the  cavitation  became  very 
"patchy"  and  intermittent.  Finally,  at  supercritical  values  of  the 
Reynolds  number  they  found  that  no  attached  cavitation  could  be  sus- 
tained on  the  model'.  Similarly  dramatic  cavitation  results  were  ob- 
served on  the  ITTC  and  the  modified  ellipsoidal  bodies.  These  obser- 
vations plus  the  original  ones  [40]  show  that, on  bodies  normally 
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having  a laminar  separation,  the  condition  of  the  separated  region 
(i.e.  extent,  steadiness)  is  a major  factor  in  determining  when  and 
what  type  of  cavitation  occu^.  This  important  conclusion  will  be  used 
extensively  in  the  discussion  of  the  results  of  the  present  experiments. 

I.  8 Freestream  Turbulence  Level  and  Cavitation  Inception 

Since  Arakeri  and  Acosta  found,  on  the  axisymmetric  bodies  they 
tested,  that  the  type  and  occurrence  of  cavitation  is  controlled  by  the 
condition  of  the  fully  wetted  flow,  it  should  be  expected  that  any  factor 
which  would  influence  the  fully  wetted  flow  would  also  influence  the 
cavitation.  One  such  factor  is  the  freestream  turbulence  level.  It  is 
well  known  [47,48]  that  the  laminar  to  turbulent  boundary  layer  tran- 
sition is  affected  by  not  only  the  intensity  but  also  by  the  power  spec- 
trum of  freestrecim  disturbances.  It  is  suspected  then  that  some  of 
the  observed  differences  between  facilities  may  be  caused  by  different 
disturbcinces  present  in  the  facility  test  section.  Unfortunately,  meas- 
urement of  turbulence  levels  in  water  is  difficult  and,  unlike  wind  tun- 
nels, there  has  been  no  great  demand  to  determine  the  freestream  tur- 
bulence levels  in  water  tunnels.  Tabulated  on  page  23  is  the  free- 
stream turbulence  level  for  the  few  tunnels  that  have  this  information 
available  [49].  Unfortunately  scale  and  spectrum  information  was 
not  available  with  these  figures. 

Whether  this  variation  in  turbulence  level  will  influence  the  fully 
wetted  flow  sufficiently  to  cause  the  inconsistencies  between  facilities 
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previously  noted  is  unknown  and  will  require  an  experimental  investi- 
gation. 


Polymer  Additives  and  Cavitation  Inception 


In  the  preceding  sections  the  effects  upon  cavitation  inception  of 
cavitation  nuclei  and  the  fully  wetted  flow  have  been  discussed.  Re- 
cently a new  factor  has  been  introduced  — drag  reducing  polymers. 
The  addition  of  a few  parts  per  million  by  weight  of  drag- reducing 
polymers  to  the  liquid  enviro.nment  has  been  found  to  suppress  the  oc- 
currence of  cavitation  in  most  situations.  Ellis  et  al  [50]  carried  out 


experiments  in  a blowdown  tunnel  to  explore  the  influence  of  polymer 

additives  on  cavitation  inception  for  two  stainless  steel  hemisphere 

nose  bodies  with  ^ inch  and  inch  diameters.  The  polymers 

tested  include  Polyox  (WSR301),  algae  (Porphy ridium  aerugineum), 

and  guar  gum  (J2FP)  at  various  concentrations.  Over  the  Reynolds 

4 5 

number  range  tested  (7  x 10  to  3.  1 x 10  ) the  incipient  cavitation 
number  was  reduced  by  as  much  as  50  percent  by  the  addition  of  poly- 
mers. Polyox  was  observed  to  be  the  most  effective  of  the  additives 
in  that  for  a given  ccncentratio.n  it  produced  the  greatest  suppression. 

Similar  experiments  have  been  carried  out  at  the  Garfield 
Thomas  Water  Tunnel  (GTWT)  as  part  of  a long  term  investigation  of 
polymer  effects  in  naval  hydrodynamics  [ 51,  52,  53,  54  ].  The  results 
of  their  tests  on  hemisphere  nose  bodies  indicate:  (1)  the  presence  of 

Polyox  (WSR301)  inhibits  cavitation  inception — maximum  inhibition  oc 
curring  when  the  polymer  is  fresh.  (2)  there  is  a definite  trend  of  re- 
duced polymer  effectiveness  as  model  size  increases  (tests  on  an  8.0 
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inch  diameter  hemisphere  body  did  not  show  a sinnificanl  cfli'ct  due  lu 
polymers).  (3)  there  seems  to  be  a correlation  between  reduced  dr.ie, 
on  a standard  body  and  cavitation  inception.  Work  on  con/ined  jets 
has  also  been  started  at  GTWT  but  no  results  are  available  yet. 

Van  der  Meulen  [ 55,  56,  57,  58,  59]  has  also  investigated  the  in- 
fluence of  polymers  on  cavitation  for  two  axisymmetric  bodies:  a 

hemisphere  nose  body  and  a "Schiebe"  body.  The  Schiebe  body  is  one 
of  a series  developed  by  placing  a distributed  source  disc  in  a uniform 
flow  and  these  bodies  are  characterized  by  having  no  laminar  separa- 
tion [88],  Van  der  Meulen  injected  polymer  into  the  body  boundary 

layer  at  the  stagnation  point.  He  found  that  polymer  suppressed  cav- 
itation on  a stainless  steel  hemisphere  nose  body  but  had  little  to  no 
effect  on  a teflon  hemisphere  nose  body.  For  the  Schiebe  body  no 
change  in  the  incipient  cavitation  index  was  measured. 

To  investigate  why  this  different  behavior  should  occur, 
van  der  Meulen  and  Oosterveld  [57]  developed  a holographic  technique 
to  simultaneously  visualize  the  viscous  flow  and  cavitation  inception. 
Experiments  with  the  hemisphere  nose  body  with  500  wppm  Polyox 
showed  that  the  polymer  eliminated  the  laminar  separation  by  appar- 
ently causing  an  early  "transition"  to  occur.  He  did  not  however  make 
any  observations  on  the  teflon  hemisphere  nose  body  but  instead  sug- 
gested that  surface  nuclei  were  responsible  for  the  different  behavior 
of  the  teflon  hemisphere  body.  Observations  on  the  "Schiebe"  body 
substantiated  that  no  laminar  separation  existed  but  instead  a 
boundary  layer  with  an  attached  laminar  to  turbulent  transition  was 
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present.  The  effect  of  the  polymer  injection  on  this  body  was  to  move 
the  transition  point  upstream  on  the  body  but  did  not  produce  any  ni»-as- 
urable  difference  in  the  inception  index. 

Other  investigations  into  the  influence  of  drag- reducing  polymers 
on  cavitation  have  been  carried  out  for  different  configur ations . Hoyt 
[60,  6 l]  has  studied  jet  cavitation  with  polymer  additives.  He  found  the 
addition  of  a few  parts  per  million  Polyox  to  the  fluid  reduced  the  in- 
cipient cavitation  number  by  a factor  of  two.  Walters  [62]  has  found 
similar  cavitation  suppression  on  discs  and  White  [63]  noted  that  bub- 
ble production  by  a small  propeller  was  reduced  by  the  presence  of 
25  wppm  Polyox.  Substantial  reductions  in  the  cavitation  number  in 
flows  of  polymer  through  a venturi  have  been  found  by  Vyazmenskii  [64]. 

Associated  with  the  suppression  of  cavitation  inception  in  these 
experiments  is  a change  in  appearance  of  the  developed  cavitation. 

Ellis  [50]  noted  that  the  developed  cavity  on  a hemisphere  became  stri- 
ated. Hoyt  [65]  found  this  in  an  earlier  experiment  and  also  noted  that 
the  cavity  collapse  was  less  violent  in  the  polymer  solution.  Brennen 
[66]  found  very  dramatic  changes  in  the  cavity  appearance  for  cavi- 
tating  spheres  and  cylinders  with  the  addition  of  different  polymers. 

To  explain  these  observed  changes  in  cavitation  appearance  and 
the  suppression  effect  two  possible  polymer-induced  changes  in  the 
fluid  behavior  have  been  suggested.  One  possibility  is  that  the  presence 
of  the  polymer  alters  the  dynamics  of  bubble  growth  and  collapse. 
Fogler  and  Goddard  [67]  analyzed  the  collapse  of  a spherical  cavity  in 
a viscoelastic  liquid  and  found  that  if  the  relaxation  time  of  the  liquid 
is  much  greater  than  the  Rayleigh  collapse  time,  the  elasticity  can 
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induce  a prolonged  vibrational  motion  of  the  cavity  boundary.  Ellis 
and  Ting  [68,69],  however,  found  no  change  in  growth  or  collapse 
rates  of  spark  generated  bubbles.  Another  possible  mechanism  is  that 
the  polymer  modifies  the  viscous  flow  around  the  body  and  thus  alters 
cavitation.  The  visual  observations  of  van  der  Meulen  [59]  and  the  re- 
sults of  Arakeri  and  Acosta  [41,42,44]  suggest  that  this  second  expla- 
nation is  the  much  more  likely. 

I.  10  Definition  of  Cavitation  Inception 

In  the  introduction  of  this  chapter  it  was  simply  stated  that  cavi- 
tation inception  occurred  when  cavities  or  bubbles  became  visible  on  or 
near  the  body  surface.  This  definition  is  not  always  easily  applied  in 
practice.  It  has  been  the  usual  procedure  to  observe  the  test  model 
under  stroboscopic  light  and  estimate  inception  by  observing  so  many 
cavitation  events  per  unit  time.  Recently  a number  of  investigators 
have  applied  the  same  idea  of  counting  cavitation  events,  but  they  have 
developed  equipment  to  do  the  counting.  At  the  present  time  there  are 
two  techniques  for  counting  cavitation  events,  one  acoustical  and  one 
optical.  The  group  at  St.  Anthony  Falls  Hydraulic  Laboratory  have  de- 
veloped the  acoustical  technique  [30].  A pressure  transducer  is  lo- 
cated inside  the  test  body.  When  a cavitation  occurrence  takes  place 
on  the  surface  of  the  body,  a spike  occurs  in  the  output  voltage.  Cavi- 
tation inception  occurs  when  the  occurrence  rate  reaches  a certain  ar- 
bitrarily defined  level. 

Ellis  [so]  and  Keller  [ 12]  have  used  optical  techniques  to  deter- 


mine cavitation  inception.  Using  a Smw  He-Ne  laser  Ellis  passed  a 
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0.  1 inch  diameter  beam  through  the  test  section.  The  beam  was  ad- 
justed to  just  graze  the  body  at  the  point  where  incipient  cavitation 
had  been  observed  to  occur.  A photocell  was  located  at  approximately 
90  degrees  to  the  laser  beam.  At  incipient  cavitation,  bubbh  s in  the 
beam  would  scatter  light  which  was  picked  up  by  the  photocell  and  the 
signal  was  recorded  on  an  oscillograph.  The  sudden  discontinuity  on 
the  oscillograph  record  was  taken  to  be  inception. 

Keller  focused  the  output  from  a strong  lamp  onto  a point  on  ei- 
ther the  top  or  the  bottom  of  the  body  where  inception  was  expected  to 
occur.  The  light  was  collected  on  the  other  side  of  the  test  section  and 
directed  into  a photocell.  Whenever  a bubble  would  pass  through  the 
focus  point  it  would  cause  scattering  and  thus  a drop  in  the  light  inten- 
sity received  at  the  photocell.  Cavitation  inception  was  said  to  occur 
when  five  events  would  occur  in  a tenth  of  a second.  It  should  be  noted 
that  the  optical  technique  determines  inception  by  the  events  at  one  par- 
ticular spot  on  the  body,  whereas  the  acoustic  method  gathers  informa- 
tion from  the  entire  body  surface. 

The  acoustic  method  is  more  useful  in  cases  where  noise  is  an 
important  design  factor  or  the  geometry  makes  use  of  the  optical  tech- 
nique difficult.  The  optical  method  has  the  advantage  that  it  is  not  dis- 
turbed by  cavitation  events  in  other  parts  of  the  water  tunnel.  The  ad- 
vantages of  these  techniques  is  that  they  provide  a method  for  consist- 
ent, repeatable  measurement  of  inception.  However,  these  techniques 
only  apply  to  cavitation  which  is  of  the  traveling  bubble  type  and  are 
not  applicable  to  cases  where  band  type  or  other  attached  forms  of  cav- 


itation occur. 
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I.  1 1 Scope  of  Present  Work 

The  liquid  environment  and  the  viscous  flow  are  known  to  he  ma- 
jor factors  in  the  inception  process.  One  purpiose  of  this  present  work 
is  to  investigate  two  aspects  of  the  environment  on  inception:  free- 

stream  turbulence  level  and  freestream  nuclei.  To  carry  out  this  in- 
vestigation the  Low  Turbulence  Water  Tunnel  (LTWT)  at  Caltech  was 
modified  so  that  the  axisymmetric  test  models  could  first  of  all  be  cavita- 
ted  and  secondly  exposed  to  a series  of  freestream  turbulence  levels. 

The  schlieren  flow  visualization  technique  was  used  to  observe  the  ef- 
fect of  the  turbulence  level  on  the  fully  wetted  flow.  Of  the  two  bodies 
tested,  one  was  found  to  be  completely  insensitive  to  the  turbulence 
level  while  the  other  was  found  to  be  very  sensitive.  Some  cavitation 
tests  were  carried  out  at  the  higher  turbulence  levels. 

To  monitor  the  freestream  nuclei  population  a pulsed  ruby  laser 
was  obtained  and  a hologram  reconstruction  system  built.  The  purpose 
of  measuring  freestream  populations  is  two-fold:  (1)  to  gather  infor- 
mation about  nuclei  distributions  in  different  facilities  and  (2)  to  in- 
vestigate the  effect  of  changing  nuclei  distributions  on  the  type  and  oc- 
currence of  cavitation.  It  was  found  in  the  present  investigation  that 
the  nuclei  distribution  obtained  in  the  LTWT  agree  well  with  those 
measured  in  other  facilities  and  that  the  presence  of  many  large  free- 
stream nuclei  can  alter  the  viscous  flow  past  a body  and  cause  travel- 
ing bubble  type  cavitation  to  occur. 

As  for  the  investigation  of  the  polymer  effect,  the  purpose  of  the 
present  investigation  was  to  determine  whether  or  not  alteration  of  the 
viscous  flow  was  responsible  for  the  observed  cavitation  suppression. 
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Polyox  (WSR301)  was  injected  into  the  boundary  layer  at  the  stagna- 
tion point  on  two  axisymmetric  bodies:  a hemisphere  nose  body  and 

the  NSRDC  body.  The  schlieren  flow  visualization  system  was  again 

« 

used  to  observe  simultaneously  the  fully  wetted  ilow  and  cavitation  in- 
ception. It  was  found  that  the  presence  of  a small  amount  of  drag- 
reducing  polymer  can  cause  an  instability  to  develop  which  results  in 
aui  early  "transition"  to  turbulence  wi'th  its  subsequent  removal  of  the 
laminar  separation.  This  work,  which  began  long  before  van  der  Meulen's 
results  were  available  [ 59  J , confirms  van  der  Meulen's  findings  at  ’ 

least  on  the  hemisphere  nose  body. 

Essentially  then,  the  present  work  is  an  investigation  into  the 
combined  effects  of  freestream  nuclei  and  the  fully  wetted  flow  on  cav- 
itation inception,  how  they  interact  with  one  another,  and  the  influence 
of  factors  which  will  change  either  the  viscous  flow  or  the  nuclei  pop- 
ulation. 


bin 
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Chapter  II 


EXPERIMENTAL  EQUIPMENT 


II.  1 Low  Turbulence  Water  Tunnel 

Two  water  tunnels  were  used  in  the  present  investigation  — 
the  Low  Turbulence  Water  Tunnel  (LTWT)  and  the  High  Speed  Water 
Tunnel  (HSWT).  All  the  preliminary  work  and  the  bulk  of  the  experi- 
ments were  carried  out  in  the  LTWT;  the  HSWT  only  being  used  to  ob- 
tain high  Reynolds  number  data  for  the  polymer  tests  on  the  hemisphere 
nose  body. 

Although  the  LTWT  has  been  around  for  a number  of  years 
(completed  in  1951),  it  has  just  recently  become  a "cavitation"  facility. 
Since  it  is  "new"  in  this  respect  and  since  it  has  a number  of  unique 
capabilities  and  advantages  as  compared  to  other  facilities,  it  is  worth- 
while describing  the  LTWT  circuit  in  detail.  The  LTWT  was  completed 
in  the  fall  of  1951  and  was  built  as  part  of  a program  to  study  the  dynamics 
of  particulates  in  fluid  suspension  [70,  7l].  In  particular,  it  was  intend- 
ed to  study  transportation  and  diffusion  of  particles  in  the  atmosphere 
and  thus  it  was  designed  to  have  a low  freestream  turbulence  level  in 
the  test  section  — and  thus  its  name.  A schematic  drawing  of  the  tun- 
nel circuit  is  shown  in  Fig.  6. 

The  working  section  of  the  tunnel  is  8-j  feet  long  and  expands 
from  a 12  inch  square  cross-section  at  its  entrance  to  a 12  inch  by  14 
inch  cross-section  at  its  exit  — the  constant  dimension  being  the  width. 
The  test  section  is  expanded  in  this  fashion  to  compensate  for  boundary 
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layer  growth  On  the  test  section  walls  and  maintain  a constant  core  ve- 
locity. The  expansion  is  too  much  however,  and  the  velocity  decreases 
by  about  4 percent  over  the  length  of  the  test  section.  The  top  and  bot- 
tom test  section  windows  are  plexiglass  and  the  side  windows  are 
glass . 

When  the  flow  leaves  the  test  section,  it  enters  a diffuser  which 
exp2uids  the  tunnel  cross-section  to  24  inches  square;  the  included  angle 
between  opposite  sides  of  the  diffuser  being  3 degrees  and  13  minutes. 

At  the  end  of  the  diffuser  is  a vaned  90  degree  elbow  which  directs  the 
flow  into  a vertical  section.  This  vertical  section  makes  a transition 
from  the  24  inch  square  section  to  a 32  inch  diameter  section  at  the 
next  elbow.  This  elbow,  which  is  also  vaned,  is  just  upstream  of  the 
pump  and  carries  the  outer  bearing  for  the  pump  shaft.  The  tunnel 
pump  is  a mixed  flow  type  and  is  driven  by  a 30  horsepower  D.  C.  elec- 
tric motor  through  a V-belt  drive.  This  combination  of  pump  and  mo- 
tor provides  a maximum  test  section  velocity  of  about  25  feet  per  sec- 
ond. 

The  iower  "leg"  of  the  tunnel  expands  from  the  3 2 inch  diameter 
cross-section  at  the  pump  exit  to  a 48  inch  square  cross-section  at  the 
next  elbow.  This  elbow  has  turning  vanes  made  of  ^ inch  steel  plate  with 
a chord  length  of  1 1 inches  and  are  spaced  3-^  inches  apart.  After  a 
short  vertical  section,  there  is  another  elbow  which  is  "heavily  " vaned  — 
the  vanes  being  only  one  inch  apart.  The  leading  edges  of  these  turning 
vanes  are  rounded  off  to  approximately  a circular  section  and  the  trail- 
ing edges  are  tapered  down  to  a inch  thickness.  These  particular 

7 

vanes  have  a chord  length  of  2-g-  inches. 
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The  section  of  the  tunnel  between  this  elbow  and  the  tunnel  noz- 
zle is  called  the  "stagnation"  section.  Two  honeycombs  and  three  tu> - 
bulence  damping  screens  are  installed  in  this  section.  The  purpose  of 
the  screens  and  honeycomb  is  to  eliminate  any  large  scale  disturbances 
in  the  oncoming  flow  and  produce  a low  turbulence  level  in  the  working 
section.  Finally,  as  can  be  seen  in  Fig.  7,  the  tunnel  cross  section  js 
reduced  from  the  48  inch  square  to  the  12  inch  square  at  the  test  section 
entrance  by  a nozzle  with  a 16:1  contraction  ratio. 

As  already  noted,  honeycombs  and  screens  are  used  in  the  stag- 
nation section  to  control  the  freestream  (.urbulence  level.  The  config- 
uration of  these  turbulence  manipulators  is  shown  in  Fig.  8.  The  up- 
stream honeycomb  is  located  immediately  downstream  of  the  elbow  and 
has  equilateral  triangular  cells  with  a 1 inch  altitude  and  is  7 inches 
deep.  Approximately  9 inches  downstream  from  the  first  honeycomb  is 
another  honeycomb.  This  one  (see  Fig.  9)  is  2 inches  deep  and  has  ^ 
inch  hexagonal  cells.  Utilizing  the  results  of  Loehrke  and  Nagib  [72], 
a screen  is  located  immediately  downstream  of  this  honeycomb.  Two 
more  screens  follow  the  first  screen  — the  spacing  between  screens  be- 
ing about  1-^  inches.  All  three  screens  have  22  meshes  per  lineal  inch 
and  are  woven  wi  th  0 . 007 5 inch  diamete r s tainle s s steel  wi re  (see  Fig. 

10).  These  screens  replace  the  original  ones  which  were  made  of  0.  018 
inch  diameter  brass  wire  and  had  18  meshes  per  lineal  inch.  These 
screens  had  an  open  areaofonly  48.3  pe  rcent(vs.  69.7  pe  rcent  for  the 
new  screen)  and  were  replaced  for  this  reason.  As  is  shown  in  Fig.  11, 
the  second  honeycomb  and  the  three  screens  are  removable. 

The  two  honeycombs  and  three  screens  produce  a low  "refe rence" 
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turbulence  level  in  the  test  section.  Turbulence  ^eneratinn  can 

be  placed  at  the  inlet  to  the  test  section  to  produce  the  desired  tt;st 
section  turbulence  levels  (see  Fig.  12).  Four  turbulence  generating 
grids  were  used  in  the  experiments  thus  providing  five  turbulence  lev- 
els for  study.  The  five  configurations  of  turbulence  manipulators  are: 

1.  no  turbulence  generating  grids 

2.  grid  #1  — ^ inch  diameter  bars  with  2 inch  mesh. 

3.  grid  #2  — inch  diameter  bars  with  1 inch  mesh. 

4.  grid  #3  — three  1 inch  diameter  horizontal  bars  on  3 inch  cente  rs . 

3 

5.  grid  #4  — 0.  025  inch  diameter  fishing  line  with—  inch  mesh. 

Grids  1 , 2 and  4 are  located  at  the  entrance  to  the  test  section.  The  distance 
from  the  grid  location  to  the  test  model  position  is  approximately  48  inches . 
Grid  #3  is  located  in  the  stagnation  section  immediately  after  the  final 
screen.  Grid#3  has  this  particular  configuration  and  location  because 
(after  much  trial  and  error)it  was  found  to  produce  a turbulence  level  which 
is  close  to  the  levels  measured  in  a number  of  other  facilities— see  Table  1.  1 . 

The  turbulence  level  in  the  test  section  at  each  configuration  of 
turbulence  manipulators  was  measured  as  a function  of  velocity  by  a 
DISA  constant  temperature  anemometer.  A DISA  wedge-shaped  hot 
film  probe  (see  Fig.  13)  was  firmly  mounted  on  the  tunnel  center  line 
(see  Fig.  14)  at  the  model  location.  The  results  of  these  measurements 
have  been  summarized  and  are  presented  in  Fig.  15.  Details  of  the  tur- 
bulence measuring  equipment  and  calculations  are  given  in  Appendix  1. 

Since  the  tunnel  was  originally  designed  and  built  without  pres- 
sure control,  a system  to  reduce  the  tunnel  pressure  so  that  the  models 
could  be  cavitated  had  to  be  added.  A Nash  5 hp  vacuum  pump  was 
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connected  by  2 inch  PV C piping  to  two  5 gallon  tanks  mounted  at  the 
high  points  at  each  end  of  the  tunnel.  The  tank  at  the  upstream  end  of 
the  tunnel  and  the  two  inch  line  to  the  vacuum  pump  can  be  seen  in  Fig. 7. 
Pressure  control  is  obtained  by  a bypass  valve  to  atmosphere  connected 
to  the  2 inch  vacuum  line.  With  this  system  the  lowest  value  of  the  cav- 
itation number  obtainable  is  approximately  0.30. 

The  tunnel  velocity  is  determined  from  the  pressure  drop  across 
the  nozzle  which  is  measured  by  a Statham  differential  pressure  trans- 
ducer and  the  tunnel  test  section  static  pressure  is  measured  by  a 
Statham  absolute  pressure  transducer  positioned  at  the  entrance  to  the  test 
section.  Each  transducer  output  is  measured  by  a Ballantine  digital  voltme- 
ter (two  models  are  used,  83 05  A and  83 10  A ) . The  two  voltmeter  displays 
can  be  frozen  simultaneously  by  a remote  switch.  Each  transduce r was 
calibrated  using  a Hg-water  manometer. 

To  inhibit  corrosion  in  the  LTWT  a passivator  (sodium  chromate ) 
has  been  added  to  the  water.  The  concentration  of  this  passivator  is 
kept  between  0.04  to  0.  2 percent  and  the  pH  is  adjusted  to  between  7.  5 
and  9.5  by  the  addition  of  KOH  [73].  Addition  of  the  inhibitor  does  not 
affect  any  properties  of  the  water,  but  it  does  give  it  a yellow  color. 

Since  the  details  of  the  HSWT  are  available  in  many  other  places 
(see,  for  example,  Knapp  et  al  [26]),  no  description  of  the  tunnel  will 
be  made  here. 

II.  2 Flow  Visualization 

To  observe  the  fully  wetted  flow  on  the  test  bodies,  the  schlieren 


flow  visualization  technique  was  used.  The  schlieren  system  used  in 
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the  present  experiments  is  shown  schematically  in  Fig.  16  and  is  essen- 
tially the  same  as  that  developed  by  Arakeri  [40].  However,  a number 
of  improvements  were  made  on  Arakeri' s equipment  to  increase  the  sys- 
tem resolution.  First,  the  spark-gap  light  source  was  replaced.  1 he 
new  light  source,  which  can  be  seen  in  Figs.  17  and  19,  has  four  O.OS 
microfarad,  10  kilovolt  capacitors  connected  in  parallel  across  two 
stainless  steel  electrodes  and  is  very  similar  to  the  unit  described  by 
Kova'sznay  [74].  It  provides  a 0,040  inch  diameter  light  source  with  a 
2 microsecond  duration  (versus  a 0.  14  inch  by  0,  20  inch  source  with  a 
5 microsecond  duration  for  the  previous  light  source).  Also,  if  required, 
the  new  unit  can  provide  two  flashes  0,  5 milliseconds  apart.  Between 
the  two  electrodes  is  a 0.  1 inch  thick  fiberglass  disc  with  a 0.012  inch 
diameter  hole  in  the  center.  The  purpose  of  t»iis  disc  is  to  confine  the 
spark  and  produce  a high  intensity  source.  The  discs  are  cut  from  cir- 
cuit board  blanks  because  it  was  found  that  it  was  the  only  readily  avail- 
able material  that  could  withstand  the  high  temperatures  created  by  the 
discharge.  For  the  same  reason  the  trigger  or  tickier  electrode  is 
made  from  tungsten  rod.  Even  with  these  precautions  the  disc  would 
have  to  be  replaced  and  the  tickler  reground  after  about  100  discharges 
to  maintain  the  light  intensity  at  the  film  plane  when  <i  magnification  of 
8 or  greater  was  used.  Maximum  magnification  obtainable  is  about  11, 
since  at  higher  magnifications  there  is  insufficient  light. 

The  other  change  to  the  original  system  to  improve  the  resolu- 
tion was  the  replacement  of  the  lenses.  The  original  lenses  were  "sim- 
ple" biconvex  ones  with  a 17  inch  focal  length.  The  replacement  lenses 
are  high  <|uality  "complex,"  corrected  lenses  which  were  originally 
used  in  aerial  cameras.  The  details  of  these  lenses  are  given  in  Fig.  16  and 
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can  be  seen  in  the  photographs  of  the  schlieren  system  in  Fi^s.  17 
through  20.  As  in  the  original  system  [40],  the  prerequisite  dtusiiy 
gradient  was  produced  by  heating  the  body  slightly  with  electric  heat- 
ers. The  subsequent  schlieren  effect  was  recorded  on  a 120  film  for- 
mat. This  allowed  the  use  of  roll  film  and  greatly  facilitated  film  han- 
dling. 

To  obtain  acceptable  schlieren  photographs  in  the  HSWT  it  was 
necessary  to  add  correction  lenses  to  the  side  teat  section  windows  to 
compensate  for  the  window  curvature.  Also,  the  windows  did  not  have 
optically  finished  surfaces  so  it  was  necessary  to  put  a thin  film  of  glyc- 
erine between  the  correction  lens  aind  the  tunnel  window  to  reduce  the 
resulting  "orange  peel"  effect.  No  such  problems  were  encountered  in 
the  LTWT  since  it  has  flat,  good  quality  glass  test  section  windows  and 
excellent  quality  schlieren  photographs  were  obtained.  An  example  of  a 
schlieren  photograph  obtained  in  the  LTWT  is  shown  in  Fig.  21  aind  can 
be  compared  with  a photograph  obtained  in  the  HSWT  shown  in  Fig.  22. 

One  final  note  about  the  schlieren  system  that  should  be  made  is 
about  the  supporting  structure.  As  can  be  seen  in  Fig.  20,  the  schlieren 
system  at  the  HSWT  was  supported  from  the  ceiling  to  isolate  it  from 
tunnel  vibrations.  This  proved  unnecessary  at  the  LTWT  and  as  can  be 
seen  in  Fig.  18  the  schlieren  system  was  mounted  directly  on  a table. 

II.  3 Injection  System 

Two  methods  for  introducing  the  polymer  into  the  flow  have  been 
used  in  the  past.  The  first  is  to  create  a polymer  "ocean"  by  adding 
polymer  additive  to  the  tunnel  water;  the  second  is  to  inject  the  polymer 
into  the  boundary  layer  from  orifices  in  the  body  surface.  Creating  a 


UUU«.AU3^ 


-38- 

polyiner  ocean  has  the  advantaf>e  that  it  eliminates  concern  over  possi- 
ble disturbances  caused  by  the  injection  — liut  that  is  about  the  only 
vantage.  Since  the  polymer  degrades  under  mechanical  action  and  with 
time  [75],  it  is  necessary  to  continuously  monitor  the  state  of  the  poly- 
mer. However,  recalling  Arndt  et  al's  [54]  experience,  obtaining  a 
valid  correlation  between  monitor  results  and  cavitation  inhibition  is  not 
a trivial  task.  Further  it  is  a major  undertaking  to  change  the  polymer 
concentration  since  this  involves  draining  the  tunnel  completely  and 
flushing  it  at  least  once.  Ellis  [50]  eased  this  particular  difficulty  by 
using  a blowdown  tunnel  which  required  much  less  water  than  a closed 
circuit  tunnel.  Another  difficulty  with  the  polymer  ocean  approach  is 
that  anomalous  readings  are  obtained  from  pressure  taps  as  noted  by 
Ellis  [50]  in  his  experiments  cind  by  Astarita  and  Nicodemo  [76], 
Metzner  and  Astarita  [77],  and  Smith  et  al  [78]  in  other  investigations 
of  polymer  effects. 

The  injection  method,  on  the  other  hand,  allows  exact  knowledge 
of  the  state  of  the  polymer,  the  injection  concentration  can  easily  be 
changed,  and  there  is  no  effect  on  tunnel  instrumentation.  However, 
the  injection  must  be  done  carefully  to  prevent  any  permanent  flow  dis- 
turbance and  also  the  tunnel  water  must  be  monitored  for  polymer  build 
up.  However,  from  van  der  Meulen's  experience  [55]  this  does  not 
seem  to  be  a problem  since  the  polymer  is  rapidly  degraded  by  mechan 
ical  action  in  the  circuit. 

Jin  Wu  [79]  has  studied  the  effectiveness  of  various  injection 
configurations  and  found  that  a porous  injector  gave  the  greatest  initial 
mixing.  However,  the  initial  estimates  of  the  polymer  injection  rates 
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required  (based  on  van  der  Meulen's  experiment)  were  too  ureat  to  pass 
through  a small  porous  injector  — a porous  injector  woul<l  have  to  he 
kept  small  to  avoid  roughness  effects.  It  was  then  derided  t«'  <^ollow 
van  der  Meulen's  example  and  inject  the  polymer  into  the  boundary  layer 

through  a hole  at  the  stagnation  point. 

As  can  be  seen  in  Fig.  23,  the  injection  system  consists  basically 
of  a supply  tank,  a filter,  a flowmeter  and  an  injector.  The  supply  tank, 
which  has  a capacity  of  about  75  liters,  stores  the  mixed  polymer  solu- 
tion until  it  is  ready  for  use.  A regulated  air  pressure  of  80  psi  is  ap- 
plied to  the  top  of  the  supply  tank  to  force  the  polymer  solution  into  the 
injector.  The  polymer  is  removed  from  the  tank  through  a fitting 
about  six  inches  from  the  tank  bottom.  As  is  shown  in  Fig.  24,  the  pfdy 
mer  solution  is  then  immediately  filtered  by  a 5 micron  fibre  filter.  A 
— inch  diameter  Poly-Flo  tube  connects  the  filter  through  a metering 
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valve  to  the  rotameter.  Finally,  a ^ inch  O.D.  by  inch  I.  D. 

Tygon  tube  connects  the  output  of  the  flowmeter  to  the  injector.  The 
rotameter  was  calibrated  for  each  concentration  used  and  interestingly 
it  was  found  that  the  calibration  curves  did  not  depart  much  from  each 
other  or  the  water  curve. 

The  first  injector  used  simply  consisted  of  a 0.  085  inch  I.  D. 
copper  tube  which  ran  the  length  of  the  test  body  and  emerged  at  the 
stagnation  point.  However,  initial  flow  visualization  tests  of  the  injec- 
tion at  the  stagnation  point,  at  flowrates  which  at  the  time  were  consid- 
ered necessary,  found  the  flow  there  to  be  very  unsteady.  The  hole  size 
was  then  increased  to  0.  167  inch  diameter  to  make  the  injection  diame- 
ter to  body  diameter  ratio  the  same  as  van  der  Meulen's.  At  the  same 
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time,  a rather  more  complex  model  injection  configuration  was  made 
first  to  minimize  any  disturbances  in  the  injection  system  which  could 
get  into  the  boundary  layer  and  second  to  standardize  the  injection 
configuration  inside  the  body  so  it  could  easily  be  fitti'd  to  other  models. 

The  injector,  shown  assernhled  inside  the  hemisphere  nose  body, 
is  presented  schematically  in  Fig,  25.  Figs  20  and  27  show  photographs 
of  the  partially  assembled  injector.  The  injector  consists  of  first  a set- 
tling section  4- inch  in  dianiet<«r  and  1—  inches  long.  This  section  was 
packed  with  porous  plastic  foam  held  in  place  by  a sintered  brass  di.sc. 
The  purpose  of  this  section  is  to  disperse  the  jet  entering  the  injector 
and  provide  a smooth  flow  into  the  9:1  contraction  which  follows.  After 
the  smooth  contraction  there  is  a tube  with  a length  to  diameter  ratio  of 
22.  The  end  of  this  tube  exits  at  the  stagnation  point. 

The  schlteren  technique  was  used  to  visualize  the  injected  flow 
and  the  maximum  ratio  of  injected  velocity  to  freestream  velocity  for 
stable  injection  was  determined  . Figs.  28  and  29  show  examples  of 
stable  and  unstable  injectant  flow  respectively  and  the  results  of  these 
tests  are  shown  plotted  in  Fig.  30.  When  experiments  were  finally 
started,  it  was  found  that  the  required  injectant  rates  were  substantially 
less  than  initially  thought.  In  fact  schlieren  photos  at  the  lower  injec- 
tion rates  showed  no  trace  of  an  injectant  standoff  like  those  shown  in 
Figs.  28  and  29.  A similar  check  on  the  injection  stability  was  carried 
out  for  the  second  body  studied  with  results  similar  to  those  found  on 
the  hemisphere  nose  body.  Thus  for  all  test  injection  rates  on  both 
bodies,  no  disturbing  effects  came  from  the  injection  process. 
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Degradation  of  the  polymer  as  it  passed  throuj;h  the  ifijection  sysleni  | 

I 

was  found  to  be  minimal. 

i 

II.  4 Models  i 

Two  axisymmetric  models  were  tested  in  this  program,  d he 
first  is  a brass  hemisphere  nose  body  2 inches  in  diameter  made  spe- 
cifically for  these  tests.  The  details  of  this  body  are  shov/n  in  Figs.  4 
and  25.  The  second  model  is  a modified  ellipsoidal  or  NSR  DC  body  and 
again  details  are  shown  in  Fig.  4.  This  body  was  provided  by  Peterson 
from  NSRDC  and  is  made  of  copper.  It  was  modified  so  that  the  inject- 
or and  the  electric  heaters  could  be  installed. 

The  models  were  supported  in  the  LTWT  by  a two-bladed  sting 
such  that  the  nose  of  the  model  was  about  six  body  diameters  upstream 
of  the  sting.  Misalignment  from  tunnel  center-line  in  both  the  LTWT 
and  the  HSWT  was  measured  to  be  about  0.  2 degrees.  In  the  HSWT  the 
models  were  supported  by  a three-bladed  sting  with  the  nose  of  the  mod- 
el again  about  six  body  diameters  upstream  of  the  sting. 

II.  5 Nuclei  Counter 

As  mentioned  in  Chapter  I,  the  lack  of  a technique  for  the  accu- 
rate appraisal  of  nuclei  character,  concentrations  and  size  distributions 
has  hampered  the  study  of  the  role  of  freestream  nuclei  in  the  inception  5 

process.  Morgan  [80],  in  his  review  of  nuclei  measuring  methods, 
enumerated  the  characteristics  that  such  a nuclei  counter  should  possess: 
it  should  measure  both  number  and  size  distributions  accurately  and  de- 
termine the  type  of  nuclei  (solid  or  bubble)  in  the  test  section  under  test 
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conditions,  i.e.  it  should  not  he  necessary  t<j  ren  c)ve  a sampK-  Irojn 
the  test  section  to  count  it  Also,  the  counter  should  operate  (Oiitir.- 
uously,  measure  over  a volume  compatible  with  body  size,  anti  be  di- 
rect i.e.  the  measurement  sl;ou!d  be  independent  of  theoretic^i!  calcu- 
lations involving  many  as s uni ptioi .s  or  extensive  calibrations.  I'inally, 
it  should  do  all  this  without  disturbing  the  flow. 

Until  recently  previous  investigators  have,  at  best,  only  been 
able  to  estimate  nuclei  distributions  by  measuring  the  amplitude  atten- 
uation of  an  acoustic  pulse  as  il  propagates  across  a water  tunnel  test 
section  [29].  And,  in  general,  most  investigators  have  been  forced  to 
rely  on  measurements  of  total  gas  content  as  an  indicator  of  nuclei  con- 
centrations. The  first  method  relies  on  theoretical  calculations  for 
interpretation  of  the  results  and  the  second  method  satisfies  none  of  the 
criteria  of  the  "perfect"  counter.  However,  with  the  development  of 
laser  technology  two  new  techniques  have  been  applied  to  this  problem. 
Both  methods  offer  large  improvements  over  the  previous  methods,  but 
they  still  do  not  satisfy  all  the  requirements  for  the  perfect  counter. 

The  first  technique  that  will  be  discussed  is  the  light  scattering 
technique  employed  by  Keller  [ 12]  which  is  described  in  Chapter  1 . 

The  advantage  of  this  teclmique  is  that  it  provides  continuous  and  ini- 
mediate  results.  The  disadvantages  being  that  it  uses  single  particle 
scattering  and  thus  has  a small  sample  volume.  Second,  the  results 
are  dependent  upon  a calibration  which  has  unknowns  in  it  associated 
with  the  nucleus  index  of  refraction.  Third,  because  of  the  Gaussian 
intensity  distribution  in  the  laser  beam  and  the  small  sample  size,  a 
particle  passing  through  the  edge  of  the  volume  will  appear  to  have  a 


-4i- 


different  size  than  it  it  j)assed  tiirouj^h  the  center  ot  the  volume. 

Fourth,  the  type  of  nuclei  is  not  th^te  rm  ined  by  this  techniqu**. 

The  second  ri’ethod  is  that  of  h'raunhofer  holography,  dhis 
technique  was  first  developed  for  and  aoplied  to  the  problem  of  sizing 
particles  in  a dynamic  aeroso’  [hll  Petersi/n  1 Ibj  and  Feklberg  and 
Shlemenson  [l7]  have  taken  this  techi.iqv.u-  and  v pplie-1  it  to  tl.e  problem 
of  cavitation  nuclei  with  goofi  success.  The  methocl  involves  a two- step 
image  forming  process.  In  the  first  step,  a sampli*  volume  is  illumi- 
nated by  a collim.ated  bean  o!  eo-nerent,  q,jasi-monochromatic  light. 
The  interference  b-'twoen  the  coherent  background  and  the  particle- 
diffracted  radiation  is  photographically  recorded  in  the  far  field  of  the 
individual  particles.  This  photographic  record  is  called  a Fraunhofer 
hologrcim  because  it  records  the  far-field  or  h'raunhofer  diffraction 
pattern.  In  the  second  step,  the  developed  hologram  is  illuminated  by 
another  collimated  beam  of  coherent,  quasi-monochromatic  light  pro- 
ducing a three-dimensional  image  of  the  original  volume.  The  recon- 
structed images  of  the  particles  < an  then  be  sized  and  counU-d  <at  the 
investigator's  leisure.  If  a more  detailed  account  of  the  theory  is  de- 
sired, Thompson  [82],  or  Parrent,  Jr.  [83]  provide  good  starting 
points . 

The  advantage  of  this  technique  is  that  it  is  an  absolute  and  di- 
rect method  for  determining  nuclei  distributions.  And  by  recording  a 
large  volume  it  provides  an  accurate  estimate  of  nuclei  distributions. 
The  disadveintages  are  that  this  technique  is  not  continuous,  does  not 
provide  immediate  results  and  the  reduction  of  data  from  the  recon- 
structed holograms  is  very  tedious. 
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II.  5.  1 Holocanie  ra 

As  menliunefl  in  the  previuus  .section,  Itie  lirst  .step  oi  itu  proc- 
ess involves  the  making  of  a tioloprarn.  I his  is  lione  hy  tlie  os<‘  o(  a 
holocamera  which  is  shovn  in  (letaii  I'l  I i . !•  i'  cotisisi--  ijasit  aliy 

of  three  components:  a light  '/Oiii  , a bear>  expander  and  a recording 
film.  Since  the  nuclei  are  small  (usually  <'  iOOu  in  diam.ete rl  and 
moving  with  velocities  of  up  to  8rn  /'sec,  a cohf»  rent  light  source  capable  of 
a very  short  duration,  high  energy  output  is  required  to  "freeze"  the 
particle  motion.  This  tv  pe  of  cutout  is  readily  available  from  a 
"Q-switched"  ruby  laser. 

The  ruby  laser  used  in  these  experiments  has  a ruby  rod  approx- 
imately 10  cm  long  X 0.  6 cm  in  diameter  which  is  "pumped"  by  a heli- 
cal xenon  flashlamp.  The  ruby  rod- flashlamp  assembly  is  located  be- 
tween two  flat,  dielectric  mirrors  which  are  spaced  1.5m  apart  and 
form  an  "optical  cavity."  Inside  the  cavity  near  each  mirror  is  an  iris 
which  is  adjusted  to  reduce  the  number  of  modes  present  in  the  laser 
output.  Also  inside  the  optical  cavity  is  a spectrophotometer  cell  con- 
taining a saturable  dye  (cry ptocy amne  in  acetone).  By  gradually  in- 
creasing the  concentration  of  cryptocyanine  the  output  of  the  laser  can 
be  reduced  to  a single  20  nanosecond  pulse  (called  Q- switching) . 

The  ruby  laser  output  first  passes  through  a beamsplitter  where 
a small  fraction  of  the  energy  is  reflected  into  a pin  diode.  The  diode 
output  is  monitored  on  a CRT  display  and  gives  the  output  intensity  as 
a function  of  time.  It  is  used  primarily  as  a check  on  Q-switching. 

The  main  portion  of  the  beam  then  enters  the  beam-expander. 
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The  beam-expander  consi.--ts  uf  a I ' irn  tucal  leiif'th  lens  which 
expands  the  beam,  a pinhole  located  in  tii-  local  plane  ol  the  l.heni 

lens  and  a 4^.ilcm  focal  length  lens  whii.h  vidamates  the  exfiandin^; 
beam.  The  purpose  of  the  pinhole  is  to  tdin  li  na*  <■  i rr*i'ulariiie.s  in  the 
beam  caused  by  the  l.bcm  len--.  Thi’  collito  t a Ixarn,  which  is  h.'lScrn 
in  diameter,  is  then  turned  ihrouyh  an  ani^h,  of  40  di^rees  Ijy  a front 
surface  mirror.  The  beam  then  passes  thront;h  the  water  tunnel  test 
section  and  finally  onto  the  recording  film.  The  holograms  are  recorded 
on  70mm  AGFA  CiEVAFR  ' lbF75  lollfilm.  This  is  a high  resolution 
film  (Z800  lines /mm)  with  maximum  sensitivity  at  the  wavelength  of 
the  laser  (6943  A).  The  filmholder  is  a 70mm  camera  back  clamped  to 
a tripod  and  held  against  the  test  section  window.  Film  exposure  is 
controlled  by  neutral  density  filters  placed  between  the  beamsplitter 
and  be  am -expander. 

II.  5.  2 Reconstructor 

The  second  step  of  the  process  involves  the  reconstruction  of 
the  original  sample  volume.  A .Spectra  Physics  5mw  He-Ne  laser  is 
used  as  the  illumination  source  and  as  can  be  seen  from  Figs.  31  and 
32,  the  beam  is  expanded  in  the  same  fashion  as  in  the  holocamera. 

The  slide  mounted  hologram  is  transported  on  an  X-Y-Z  vernier  car- 
riage along  the  axis  of  the  collimated  beam.  Reconstructed  nuclei  im- 
ages are  imaged  through  a lOX  Wild  microscope  objective  onto  a vide- 
con  eind  displayed  on  a closed  circuit  television  monitor.  The  total 
magnification  from  actual  nuclei  to  the  monitor  is  160X.  The  particles 
are  sized  directly  from  the  monitor  by  using  a reticle  on  the  screen. 
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Fig.  33  shows  a photograph  of  Ihi-  reconstruction  system. 

To  obtain  good  estimates  of  nuclei  concentrations  over  tlie  r 
corded  size  range  (l0-500u  diameter)  two  different  size  sample  vol 
umes  from  the  reconstructed  space  are  an.iiy-ed.  For  th(  more  numer- 
ous small  particles  « 50, j diame'er)  a d.Ocm^  volume  is  .-xamined, 
while  for  the  more  scarce  large  parti,  les  (-  50u  diameter)  a 3o.9cm^ 
volume  is  studied.  Typically.  200  particles  per  hologram  are  ana- 
lyzed and  this  requires  about  one  man-day.  As  was  also  found  by- 
Peterson  [15],  the  limit  of  r-e solvability  between  a dirt  particle  and  a 
gas  bubble  is  20-25a.  Although  the  smallest  nuclei  that  can  be  detected 
in  the  present  system  are  5u  diameter,  the  background  noise  in  the 
reconstruction  system  prevents  any  reliable  estimate  of  nuclei  concen- 
trations for  particles  less  than  lOu  diameter. 
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Chaprer  111 


EXPERIMENTAL  PROCEDURES  AND  REDUCTION  OF  DA  J A 


m.  1 p reparation  oi  Polymer  S(.>!ut-ion a 

The  polymer  used  in  all  the  expe  r T.-ienta  is  Polyox  (WSR  301) 
and  was  obtained  from  Union  Carbide.  Before  mixing,  the  polymer  was 
first  sifted  through  a U S.  Series  Equivalent  Mesh  No.  45  wire  screen 
which  has  a mesh  opening  of  0.  01  >8  inches.  This  was  done  to  remove 
any  large  aggregates  in  t:it  pvjwder  because  it  was  found  that  they  usu- 
ally did  not  dissolve.  The  polymer  was  then  mixed  according  to  the  fol- 
lowing procedure: 

1.  The  polymer  was  put  into  suspension 
in  a small  beaker  of  alcohol. 

Z.  The  suspension  was  then  slowly  poured 
into  a beaker  of  water  which  was  being 
gently  stirred.  This  forms  a highly  con- 
centrated solution. 

3.  This  concentrated  solution  was  left  to 
dissolve  overnight. 

4.  At  least  one  hour  before  the  test,  the 
concentrated  solution  was  added  to 
enough  water  to  rnake  a 100  liter  batch. 

5.  Just  before  the  test  was  to  start  the 
solution  was  siphoned  from  the  mixing 
tank  into  the  injection  system  holding 
tank. 

The  proportions  for  each  concentration  are  given  in  the  table  below: 
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TABLE  111  1 

Preparation  of  Poly 

mer  Scdutions 

Cone  'n 

Polyox 

Alcohol 

V'oluriie  of  I '.or.iJ 

(wppm) 

(gni) 

(cc  ) 

(liter; 

20 

Z 

Z ‘j 

* 

50 

5 

■) 

4 

100 

10 

12:> 

4 

500 

50 

500 

It, 

Approximately  Z liters  '.vcre  taken  from  each  batch  and  tested 
in  a turbulent  flow  rheometc”  to  moi.itor  the  consistency  of  the  nuxing 
technique.  The  turbulent  flow  rheometer  [84]  consists  basically  ol  a 
short  length  of  tubing  through  which  a given  amount  ol  polymer  solution 
is  forced.  The  time  required  to  do  this  is  measured  and  compared  to 
the  times  of  other  batches  of  the  same  concentration 

III . 2 Preparation  of  Tunnel  Water 

Since  the  LTWT  has  no  resorber,  it  was  necessary  to  strongly 
deaerate  the  tunnel  water  before  any  cavitation  tests  could  be  done  If 
this  were  not  done,  tremendous  numbers  of  freestream  bubbles  could 
be  seen  in  the  test  section.  The  water  was  deaerated  by  pulling  a vac- 
uum on  the  tunnel  with  the  Nash  pump  and  at  the  same  time  slowly  re- 
circulating the  water  By  this  method  the  total  air  content  could  be  re- 
duced from  a saturated  value  of  approximately  14  ppm  to  7 ppm  (the 
lowest  value  obtained)  in  about  16  hours  of  continuous  operation.  Air 
content  level  in  the  HSWT  was  reduced  by  lowering  the  tunnel  pressure 
until  the  test  body  would  supercavitate  at  a tunnel  velocity  of  30  fps. 
Every  15  minutes  the  tunnel  would  be  stopped  and  the  collected  air 
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would  be  bled  from  the  tunnel  high  points.  By  this  method  the  air  con- 
tent was  reduced  to  between  9-10  ppm  and  it  remained  constant  for  thi- 
duration  of  the  tests.  The  air  content  was  measured  tiy  a Delta  Scien- 
tific 1010  dissolved  oxygen  meter  which  was  occasionally  ctu;cke<l  with 
a van  Slyke  apparatus.  The  ai  i content  was  '-liecked  belore  and  after 
each  day's  testing. 

In  the  LTWT,  in  which  the  bulk  of  the  tests  we  re  done  and  which 
has  a small  capacity  (5000  gallons),  eventual  contamination  of  the  water 
by  the  polymer  was  a worry.  Van  der  Meulen  [55]  used  a turbulent 
flow  rheometer  similar  to  the  one  used  in  these  experiments  to  moni- 
tor his  tunnel  water  and  his  rheometer  results  indicated  that  the  poly- 
mer had  no  effect  on  the  tunnel  water.  Van  der  Meulen  assumed  this 
was  due  to  rapid  mechanical  degradation  of  the  polymer  as  it  circulated 
through  the  tunnel.  However,  in  light  of  Arndt  et  al's  [54]  finding  that 
a test  body  experienced  a drag  reduction  even  though  a rheometer  test 
of  the  solution  indicated  no  drag  reduction,  this  method  is  suspect.  Un- 
fortunately, this  finding  did  not  come  to  my  attention  until  after  the 
completion  of  the  experiments.  But  cavitation  inception  tests  with  no 
injection  were  performed  each  day  of  the  testing  and  these  values 
showed  no  downward  trend.  Also,  schlieren  photographs  of  the  lami- 
nar separation  showed  no  change  over  the  duration  of  the  tests.  Thus 
the  contamination  of  the  tunnel  water  had  no  significant  effect.  Poly- 
mer contamination  in  the  HSWT  was  not  a problem  since  very  little 
polymer  was  injected  into  the  circuit  (tests  were  held  over  a one  week 
period  vs  a several  month  period  for  the  LTWT)  and  the  HSWT  has  a 
much  larger  capacity  (70,000  vs  5000  gallons). 


Ill,  3 Photography  and  Holography 


Single  exposure  schlieren  photographs  were  taken  ot  the  lent 
model  in  silhouette  at  magnifications  from  4 to  11.  The  low  intensity 
light  available  at  the  film  plane  at  magnifications  of  8 or  greater  re- 
quired the  use  of  high  speed  film  and  special  developing.  Royal-X-Pan 
film  with  an  ASA  rating  of  1250  was  used  and  the  fUm  was  "pushed”  in 
the  development  by  using  Acufine  developer.  The  holography  film 
AGFA  10E75  had  no  special  developing  procedures  except  for  a rather 
short  development  time  (3  minu'es). 

III.  4 Test  Procedure  — LTWT 

In  a typical  cavitation  test  involving  polymer  injection,  the  tun- 
nel velocity  and  the  injection  rate  were  first  adjusted  to  the  desired 
values.  The  tunnel  pressure  was  then  lowered  until  inception  took 
place.  Inception  was  determined  visually  under  stroboscopic  light  At 
the  instant  of  inception  two  switches  were  closed.  When  these  switches 
close  the  schlieren  spark  source  and  the  ruby  laser  are  fired  and  the 
DVM  readouts  are  frozen  — all  simultaneously.  Thus  at  inception,  a 
schlieren  photograph  is  taken,  a hologram  made  and  the  tunnel  velocity 
and  pressure  recorded.  The  test  had  to  take  less  than  forty  seconds 
because  by  thac  time  bubbles  generated  at  the  pump  would  reach  the 
test  section  and  dramatically  change  the  conditions  in  the  working 
section.  After  each  test  the  tunnel  pressure  was  raised  and  the  tunnel 
allowed  to  circulate  for  five  minutes  before  the  next  test  began.  This 
recess  between  each  test  was  required  to  let  the  ruby  laser  cool  down 
and  also  let  the  air  go  back  into  solution  or  rise  to  the  tunnel  high 
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points  and  be  removed. 

III.  5 Test  Procedure  — HSWT 

In  the  HSWT  the  same  general  procedure  as  used  in  the  LTWT 
was  followed.  Toggle  valves  in  the  manometer  lines  were  used  to 
"freeze”  the  velocity  and  pressure  readings  at  inception.  After  the 
readings  were  recorded,  the  pressure  was  increased  until  the  cavita- 
tion disappeared  (desinence)  and  the  velocity  and  pressure  readings 
were  again  frozen  and  recorded  Since  no  holograms  were  made  and 
since  the  HSWT  has  a resorber,  it  was  not  necessary  to  wait  five  min- 
utes between  each  test. 

III.  6 Calculation  of  Tunnel  Velocity  and  the  Cavitation  Number 

As  there  is  a boundary  layer  growth  along  the  walls  of  the  work- 
ing section  of  the  tunnel,  it  is  necessary  to  calibrate  the  tunnel  both 
bare  and  with  the  model  installed  in  order  to  determine  the  static  pres- 
sure and  the  velocity  at  the  location  of  the  model.  Such  calibrations 
have  been  done  by  Acosta  and  Hamaguchi  [18]  for  the  HSWT.  Working 
formulas  for  the  tunnel  speed  and  cavitation  number  are  given  in  the 
same  reference  and  were  used  in  the  reduction  of  the  raw  data.  Such 
information  was  not  available  for  the  LTWT  and  the  calibrations  had  to 
be  done.  The  results  of  the  calibrations  and  the  working  formulas  for 
working  section  velocity  and  pressure  are  given  in  Appendix  II. 


n. 
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Chapter  IV 


PRESENTATION  AND  DISCUSSION  OF  THE  FULLY  WETTED  RESULTS 


As  mentioned  previously,  the  influence  of  freestrearn  turbu- 
lence level  and  the  presence  of  a drag  - reducing  polymer  in  the  bound- 
ary layer  upon  the  fully  wetted  flow  about  two  axisymmetric  bodies 
was  studied  by  schlieren  flow  visualization.  The  effects  of  these  fac- 
tors upon  the  viscous  flow  were  quantified  by  measuring  the  position, 
maximum  height  and  the  length  of  the  laminar  separation  which  exist- 
ed on  both  bodies  over  their  respective  test  Reynolds  number  ranges. 
The  length  of  the  separated  bubble  (L)  was  taken  to  be  the  distance 
from  the  position  of  separation  to  the  point  on  the  separated  streamline 
where  an  oscillation  first  begins.  The  height  of  the  separated  bubble  (H) 
was  defined  as  the  maximxim  bubble  height.  These  quantities  are  de- 
fined in  Fig.  34  and  were  measured  directly  from  the  negative  of  the 
schlieren  photograph  with  the  aid  of  a scale  (or  reference)  negative. 

The  definition  of  the  position  of  transition  on  the  separated  streamline 
is  the  same  as  that  for  the  "end"  of  the  separation  bubble,  i.e.  the  end 
of  the  separated  region  and  the  position  of  transition  coincide.  The 
above  quantities  are  illustrated  in  the  schlieren  photograph  of  Fig.  21. 

IV.  I Freestrearn  Turbulence  Level  - Presentation  of  Results 

A sequence  of  photographs  showing  the  laminar  separation  on 
the  hemisphere  nose  body  for  a constant  freestrearn  velocity  but  grad- 
ually increasing  turbulence  level  is  presented  in  Fig.  35.  Similar  se- 
quences for  the  NSRDC  body  are  presented  in  Figs.  36  and  37.  As  can 
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be  seen  in  Fig.  35,  the  laminar  separation  on  the  hemisphere  nose 

body  remained  unchanged  as  the  freestream  turbulence  level  was  in 

creased.  In  Figs.  36  and  37,  however,  it  can  be  seen  that  the  viscous 

flow  on  the  NSRDC  body  was  drastically  influenced  by  the  increasing 

freestream  disturbance  level.  The  freestrc;ani  velocity  in  Fig.  36  is 

5 

about  9 . 7 feet  per  second  1.6  X 10  ) and  in  Fig.  37  is  approxi- 
mately 15  feet  per  second  = 2.  5 x 10^).  The  results  of  measure- 

ments of  the  laminar  separation  quantities  (position,  length,  maximum 
height)  are  presented  in  Figs.  38  through  42.  These  quantities  have 
been  non-dimensionalized  by  the  body  diameter  and  plotted  against  body 
Reynolds  number  with  freestream  turbulence  level  as  a parameter. 
Each  plotted  value  is  the  average  of  at  least  three  and  usually  four 
observations  and  the  maximum  scatter  is  approximately  ±25  percent 
of  the  mean. 

It  is  readily  apparent  from  the  schlieren  photographs  that  the 
hemisphere  nose  body  is  insensitive  to  these  freestream  disturbances 
while,  by  comparison,  the  NSRDC  body  is  very  susceptible.  The  lack 
of  response  of  the  flow  about  the  hemisphere  nose  body  to  these  im- 
posed disturbances  is  surprising  in  view  of  the  strong  effect  on  the 
NSRDC  body  which  is  of  a similar  type.  To  see  if  this  behavior  would 
be  expected,  the  recent  literature  on  the  influence  of  pressure  grad- 
ient and  freestream  turbulence  level  upon  boundary  layer  transition 
was  consulted. 

IV.  2 Transition  Prediction 

There  are  two  basic  approaches  to  the  prediction  of  laminar  to 
turbulent  boundary  layer  transition.  In  the  first  approach  theory  and 
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experiment  have  been  combined  to  provide  an  empirical  relation  for 
transition  prediction  [for  example  see  Ref.  [45]].  The  second  ap- 
proach is  theoretical  and  involves  the  determination  of  the  linear  re- 
sponse of  a laminar  boundary  layer  to  small  usually  two  dimensional 
disturbances,  i.e.  the  solution  of  the  Orr-'^ommerfeld  equation.  The 
advantage  of  the  first  approach  is  that  it  always  gives  a result,  how- 
ever the  validity  of  that  result  is  not  always  known.  The  second  ap- 
proach attempts  to  predict  transition  from  first  principles,  however 
the  lack  of  a full  understanding  of  the  mechanism  of  transition  prevents 
meaningful  calculations  except  in  cases  of  simple  freestream  condi- 
tions and  body  geometries. 

The  available  experimental  data  and  empirical  correlations 

pertaining  to  the  combined  influence  of  freestream  turbulence  level 

and  pressure  gradient  upon  transition  have  been  reviewed  by  Hall  and 

Gibbings  [48].  They  have  summarized  this  information  on  a graph  of 

the  transition  Reynolds  number  (RCq)  versus  the  Polhausen  parameter 

(X  ) for  different  freestream  turbulence  intensities  One  correlation 

9 

which  is  typical  of  the  data  presented  in  Hall  and  Gibbing's  figure  is 
that  of  van  Driest  and  Blumer  [85].  Using  the  correlation  proposed  by 
van  Driest  and  Blumer,  the  critical  Reynolds  numbers  (^le^)  for  the 
present  bodies  were  calculated  for  various  turbulence  levels. 

Another  correlation  mentioned  by  Hall  and  Gibbings,  which  will 
be  used  extensively  in  this  chapter,  is  that  of  Jaffe  et  al  [45].  Basically 
this  technique  involves  the  calculation  of  the  spatial  growth  of  small 
disturbances  in  the  boundary  layer  along  the  surface  of  the  body. 
Transition  is  said  to  occur  when  the  disturbance  has  been  amplified  bj- 
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a factor  of  e . This  method,  however,  does  not  take  into  at  count 
freestream  disturbances  and  thus  the  results  only  apply  to  a /.ero  turbu- 
lence level  flow.  The  details  of  the  application  of  this  method  to  the 
present  bodies  is  given  in  Appendix  III. 

IV . 3 Results  of  Transition  Calculations  and  Comparison  with  Present 
Experimental  Results 

Results  of  transition  calculations  using  the  correlations  of  Jaffe 
et  al  and  van  Driest  and  Blumer  are  presented  in  Fig.  43  along  with 
the  present  experimental  results  and  those  of  Arakeri  and  Acosta  [41  ] 
In  this  figure  the  value  of  the  critical  Reynolds  number  has  been  plotted 
vs.  freestream  turbulence  level  for  both  the  NSRDC  and  the  hemisphere 
nose  body.  The  distinction  between  Eqs.  10  and  1 1 of  van  Driest  and 
Blumer  is  that  Eq.  10  was  developed  assuming  a Polhausen  velocity 
distribution  while  in  Eq.  11  Falkner-Skan  profiles  were  used. 

For  the  NSRDC  body  the  correlation  of  van  Driest  and  Blumer 
agrees  well  with  the  present  observed  results  for  disturbance  levels 
less  than  one  percent.  Above  this  disturbance  level  the  agreement  is 
very  poor.  For  the  hemisphere  nose  body  the  correlation  predicts  criti- 
cal Reynolds  numbers  almost  identical  to  those  of  the  NSRDC  body  and 
thus  varies  substantially  from  the  observed  results.  The  Jaffe  et  al 
correlation  prediction  of  the  critical  Reynolds  number  for  the  NSRDC 
body  also  agrees  well  with  experiment.  However,  for  the  hemisphere 
nose  body  this  correlation  predicts  a critical  Reynolds  number  of 
5 X 10^  a factor  of  10  greater  than  that  predicted  by  the  van  Driest  and 
Blumer  correlation.  Since  Arakeri  and  Acosta  [41  ] found  the  laminar 


r 


-56- 

separation  to  still  be  present  at  a Reynolds  number  of  1 x 10^,  it  is  be- 
lieved the  higher  value  is  more  near  correct.  The  reason  for  the  dif 
ference  in  the  critical  Reynolds  number  predicted  by  the  two  correla- 
tions for  the  hemisphere  nose  body  is  believed  to  be  due  to  the  local 
nature  of  the  van  Driest  and  Blumer  correlation.  This  correlation 
does  not  take  into  account  the  stability  or  previous  history  of  the  bound- 
ary layer  as  the  Jaffe  et  al  correlation  does.  In  view  of  the  last  state- 
ment it  is  believed  that  the  good  agreement  between  the  van  Driest  and 
Blumer  prediction  of  the  cri.  cal  Reynolds  number  and  the  experimen- 
tally observed  value  on  the  NSRDC  body  is  fortuitous. 

IV . 4 Boundary  Layer  — Freestream  Turbulence  Interaction 

Neither  of  the  correlations  predicted  the  lack  of  response  of  the 
hemisphere  nose  body  to  the  increasing  freestream  turbulence  level. 
The  van  Driest  and  Blumer  correlation  failed  because  it  is  inadequate 
and  the  Jaffe  et  al  correlation  does  not  include  the  effects  of  freestream 
turbulence  level.  As  pointed  out  in  the  most  recent  review  of  boundary 
layer  stability  and  transition  (Reshotko  [86],  1976),  one  of  the  basic 
problems  in  transition  prediction  is  lack  of  full  understanding  of  the 
mechanism  by  which  freestream  disturbances  interact  with  the  bound- 
ary layer.  Further  Reshotko  emphasizes  that  a knowledge  of  the  wave 
number  and  orientation  spectra  of  the  disturbance  environment  is  of 
great  importance.  Wells  and  Spangler  [47],  for  example,  have  shown 
that  transition  on  a flat  plate  is  sensitive  both  to  the  nature  of  the  dis- 
turbance and  its  frequency  content  as  well  as  the  intensity.  Unfortu- 
nately this  type  of  information  was  not  able  to  be  obtained  in  the 
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present  experiments.  The  best  that  can  be  done  is  to  make  an  esti- 
mate of  the  ratio  of  the  enerj>y  available  at  the  respective  critical 
boundary  layer  frequencies  for  each  body.  Using  the  approximati;  cal 
culation  technique  of  Jaffe  et  al  [45]  estimates  of  the  critical  boundary 
layer  frequencies  were  made  and  are  presented  in  Table  IV.  1.  Mak- 
ing use  of  Tsuji's  [87]  measurements  of  turbulence  spectra  behind 
grids  it  was  estimated  from  his  power  spectra  that  there  is  approxi- 
mately 60  times  as  much  energy  available  at  the  critical  frequency  of 
the  NSRDC  body  than  at  the  critical  frequency  of  the  hemisphere  nose 
body  at  a Reynolds  number  of  2.  5 x 10^  . 

Another  factor  which  perhaps  should  be  taken  into  consideration 
is  the  distance  available  for  the  growth  of  the  disturbances  on  the  body . 
On  the  NSRDC  body,  the  distcuice  from  the  point  of  neutral  stability  to 
the  position  of  separation  is  0.42  diameters,  whereas  for  the  hemi- 
sphere nose  body  it  is  0.07  diameters'.  Thus  the  distance  available  for 
growth  on  the  NSRDC  body  is  six  times  that  on  the  hemisphere  nose  body 
since  both  bodies  have  the  same  diameter.  Using  values  of  phase  ve- 
locities calculated  by  W'  .zan*  for  aSchiebe  body  [88]  andassuming  the 
values  for  the  NSRDC  and  the  hemisphere  nose  bodies  are  about  the 
same,  the  growth  rates  of  disturbances  were  calculated.  It  was  esti- 
mated that  a 5 percent  growth  in  amplitude  per  wavelength  for  the 
NSRDC  body  and  an  18  percent  growth  in  amplitude  per  wavelength  for 

the  hemisphere  nose  body  would  be  required  to  produce  transition  at 

5 

the  separation  point  at  a Reynolds  number  of  2,  5 x 10  . Whereas  at 
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each  body's  estimated  critical  Ileynolds  number  with  no  freeslrearn 
disturbances,  only  a 1-2  percent  growth  in  amplitude  per  wavokii^^ih 
is  required.  The  above  calculations,  however,  must  h<!  viewed  scep- 
tically since  the  small  disturbance  parallel  flow  theory  is  undoubtedly 
invalid  at  these  amplification  rates.  However,  the  results  do  indicate 
that  the  freestream  energy  available  and  the  time  available  for  the 
growth  of  disturbances  is  much  less  for  the  hemisphere  nose  than  the 
NSRDC  nose  and  these  reasons  probably  account  for  the  insensitivity 
of  the  hemisphere  nose  body  to  the  freestream  turbulence  levels  in  the 
present  experiments.  Yet,  the  above  discussion  is  still  highly  spec- 
ulative in  view  of  lack  of  understanding  about  the  interaction  of  free- 
stream disturbances  with  the  boundary  layer. 

In  this  section  the  influence  of  freestream  turbulence  level  upon 
the  fully  wetted  flow  has  been  discussed.  In  the  next  section  the  influ- 
ence of  small  quantities  of  polymer  in  the  boundary  layer  upon  the  ful- 
ly wetted  flow  will  be  discussed. 

IV.  5 Influence  of  Polymer  Injection  on  the  Fully  Wetted  Flow 
Presentation  of  Results 

Before  conducting  any  polymer  experiments,  the  effect  of  the 
injection  itself  upon  the  flow  was  tested.  The  sequences  of  schlieren 
photographs  presented  in  Figs.  44  and  45  show  the  effects  of  water  in- 
jection upon  the  laminar  separation  on  the  hemisphere  nose  and  NSRDC 
bodies.  It  is  readily  observed  in  these  photographs  that  the  injection 
process  has  no  effect  on  the  viscous  flow  at  the  injection  rates  used  in 
the  polymer  tests. 


Sequences  of  schlieren  pholo^^raphs  showinn  the  effects  of  Krad- 
ually  increasing  injection  rate  on  the  laminar  separation  cjn  the  henii 
sphere  nose  body  for  two  different  concentrations  are  given  in  Figs. 

46  cind  47.  Similar  sequences  are  given  in  Figs.  48  and  4’)  for  the 
NSRDC  body  at  the  same  concentration  but  two  different  velocities. 

In  each  sequence  the  transition  point  is  seen  to  move  upstream  as  poly- 
mer is  added.  The  forward  movement  of  the  transition  point  on  the 
free  shear  layer  causes  a reduction  in  size  of  the  laminar  separation. 
When  the  injection  rate  is  such  that  the  transition  and  separation  points 
coincide,  the  laminar  separation  is  eliminated  and  replaced  by  a turbu- 
lent boundary  layer.  If  the  injection  rate  is  increased  beyond  this  lev- 
el, the  position  of  transition  moves  further  upstream  as  can  be  seen  in 
photographs  47(e)  and  49(c).  However,  when  this  happened  the  position 
of  transition  could  not  be  detected  because  of  the  limited  resolution  of 
the  present  schlieren  system.  Summaries  of  the  separation  dimen- 
sions are  presented  in  Figs.  49  through  54.  The  separation  dimen- 
sions have  been  nondimensionalized  by  the  body  diameter  and  plotted 
against  the  body  Reynolds  numbers  for  various  values  of  the  param- 
eter "G".  This  parameter  is  the  ratio  of  the  mass  flowrate  of  poly- 
mer in  the  boundary  layer  to  the  boundary  layer  mass  flowrate  calcu- 
lated at  the  position  of  separation,  i.e.  : 


G = 


cQ 

U 6* 
oo  s 


where  c is  the  polymer  concentration  in  wppm,  Q is  the  polymer 
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solution  injection  rate,  is  the  freestream  velocity,  and  6^  is  the 

boundary  layer  displacement  thickness  at  the  position  of  separation. 
Why  this  parameter  was  used  will  be  explained  in  the  next  chapter. 

IV . 6 Comparison  of  Present  Results  with  Those  of  van  der  Meulen 

Since  the  present  experiments  are  closely  related  to  those  of 
van  der  Meulen  [59],  it  would  be  informative  to  compare  the  fully 
wetted  results  of  these  two  separate  investigations.  Van  der  Meulen 
studied  the  influence  of  dilute  polymer  solutions  upon  the  fully  wetted 
flow  and  upon  cavitation  inception  for  two  axisymmetric  bodies.  As  in 
the  present  case,  the  polymer  was  injected  into  the  boundary  layer  at 
the  stagnation  point.  The  two  bodies  studied  were  a 10mm  diameter 
hemisphere  nose  body  and  a blunt  nose  body  of  the  series  described  in 
Schiebe's  report  [88]  having  a final  diameter  of  9.88mm  and  a min- 
imum pressure  coefficient  of  -0.75.  The  test  Reynolds  number  range 

5 

was  1.2x10®  to  2.6x10  . 

To  observe  the  flow  vain  der  Meulen  used  pulsed  ruby  layer 
holography.  However,  to  make  the  flow  visible  he  added  salt  to  the 
polymer  solution.  The  injected  liquid  was  then  a 2 percent  - BOOwppm 
Polyox  (WSR  3011  solution.  However,  the  density  stratification  pro- 
duced in  the  boundary  by  the  injection  of  this  solution  is  known  to  be 
destabilizing  on  the  curved  sections  of  these  bodies  (see  for  example 
Ref  [89]).  The  magnitude  of  the  effect  relative  to  that  of  the  polymer 
is  unknown.  However,  the  good  agreement  between  the  laminar  sep- 
aration dimensions  on  the  hemisphere  nose  body  measured  by 
Arakeri  [40]  and  those  obtained  by  van  der  Meulen  when  he  injected  a 
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aalt  solution  suggests  that  the  effect  is  small  on  the  hemisphere  nose 
body.  The  magnitude  of  the  effect  on  the  "Schiebe”  body  is  unknown 
and  should  be  kept  in  mind  when  examining  the  ri-sulls  on  tiiat  bofly. 

On  the  hemisphere  nose  body  van  der  Meulen  oi)serve<l  lha'  liie 
injection  of  the  salt-polymer  solution  eliminated  the  laminar  separa- 
tion by  causing  an  early  transition  to  a turbulent  non-sejiarating  bound- 
ary layer.  On  the  "Schiebe"  body,  on  which  no  laminar  separation 
was  observed,  the  laminar  to  turbulent  transition  point  was  found  to 
move  to  a new  position  upstream  of  the  no-injection  position.  Com- 
parison of  quantitative  data  is,  however,  not  very  fruitful . Van  der  Meulen 
always  injected  at  a rate  such  that  the  injection  layer  thickness  at  the 
tangent  point  on  the  body  was  between  2 and  4 microns  (i  e.  values  of 
G between  40  x lO"^  and  70  y 10‘^).  At  these  injection  rates  the  lam- 
inar separation  on  the  hemisphere  nose  was  always  eliminated  and  thus 
no  comparison  of  laminar  separation  dimensions  with  polymer  present 
coul'"’  be  made.  Since  the  "Schiebe"  body  geometry  was  not  tested  in 
the  present  experiments,  no  quantit  ‘ive  comparison  can  be  made. 

The  results  of  the  prtse  tex,  iments  confirm  van  der  Meulen's 
observations  (on  the  hemisphere  nose  body  at  least)  that  the  presence 
of  small  quantities  of  Polyox  in  the  boundary  layer  causes  a destabiliz- 
ing effect.  Such  an  effect  is  rather  contrary  to  the  general  impression 
obtained  from  the  available  literature  on  the  effects  of  drag-reducing 
polymer  on  fluid  friction  (see,  for  example,  Hoyt  [90]).  For  this  rea- 
son we  summarize  here  from  the  literature  the  influence  of  these  poly- 
mers on  flow  past  bodies  having  either  a laminar  or  turbulent  separa- 


tion or  an  attached  laminar  to  turbulent  separation. 
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IV . 7 Other  Work 

The  influence  of  drag-reducing  polymers  upon  separation  has 
been  investigated  mainly  in  drag  - reduction  experiments  with  spheres 
and  cylinders.  Sarpkaya  et  al  [9l]  have  reviewed  the  available  litera- 
ture on  this  topic:  it  is  worthwhile  to  mention  a few  of  these  investiga 
tions  here  — especially  those  in  which  some  flow  visualization  was 
done,  Lang  and  Patrick  [92]  observed  the  flow  about  freely  falling 

# 

spheres,  cylinders  and  cones  in  solutions  of  Polyox  in  water.  They 
found  the  drag  on  a sphere  in  the  subcritical  Reynolds  number  range 
was  considerably  reduced  by  the  addition  of  polymer.  Flow  visualiza- 
tion by  dye  showed  that  the  laminar  boundary  layer  separation  was 
moved  further  downstream  and  that  wake  size  was  decreased.  White  [93] 
also  carried  out  free-fall  experiments  with  spheres  in  a polymer  solu- 
tion and,  like  Lang  and  Patrick,  noted  a large  drag  reduction  at  sub- 
critical  Reynolds  numbers.  He  also  concluded  that  a delayed  separa- 
tion and  a smaller  wake  size  were  responsible  for  the  reduction. 
Brennen  [66]  studied  the  influence  of  polymer  on  fully  developed  cavi- 
ty flow  past  two  spheres  and  a cylinder.  He  observed  that  the  polymer 
caused  a distortion  in  the  separation  line  and  also  on  the  cavity  surface 
for  all  three  bodies.  However,  whereas  the  separation  line  on  the 
spheres  moved  downstream,  the  separation  line  on  the  cylinder  did  not 
change  position.  Brennen  attributed  the  above  observations  to  a poly- 
mer induced  instability  in  the  wetted  surface  flow  around  the  headform 
which,  upon  convection,  dramatically  altered  the  cavity  surface  appeal^ 
ance.  Sarpkaya  et  al  [9l]  investigated  the  influence  of  polymer  solu- 
tions on  flow  past  cylinders.  Using  dye  injection  to  visualize  the  flow 


it  was  observed  that  the  separation  point  moved  downstream  slightly 
and  became  distorted.  Sarpkaya  also  claimed  that  the  dye  visualiza- 
tion tests  revealed  a laminar  to  turbulent  transition  just  ahead  of  sep 
aration. 

In  each  of  the  above  experiments  a delay  in  separation  was  ob- 
served. However,  it  is  not  clear  that  this  delay  is  due  to  an  early 
polymer-induced  transition  to  a more  separation- resistant  turbulent 
boundary  layer.  The  distortion  of  the  separation  line  noted  by  Brennen 
and  Sarpkaya  et  al  suggest  that  this  may  be  the  case.  To  obtain  addi- 
tional information  about  the  influence  of  polymer  on  separations,  we 
turn  to  the  effect  of  polymer  solutions  on  the  flow  past  hydrofoils. 

IV . 8 Flow  of  Polymer  Solutions  Past  Lifting  Surfaces 

The  flow  field  about  a lifting  surface  can  be  changed  substanti- 
ally by  changing  the  angle  of  attack  of  the  foil  and  thus  provides  a bet- 
ter opportunity  to  study  polymer  effects  on  separation  and  treinsition 
than  does  a cylinder  or  sphere.  Most  investigators  [94,  95,  96,  97, 
98]  have  found  that  the  presence  of  polymer  decreases  the  performance 
of  a lifting  surface.  Latto  and  Czaban  [99]  found,  however,  an  in- 
crease in  turbomachinery  performance  with  the  addition  of  polymer. 
Lehman  and  Suessmann  [ lOO],  who  did  both  polymer  injection  and 
polymer  "ocean"  tests  with  hydrofoils,  found  the  lift  could  be  either 
increased  or  decreased  depending  on  the  injection  process  and  the 
drag,  although  generally  reduced,  could  in  some  circumstances  actual- 
ly be  increased.  Fruman  et  al  [lOl,  102]  also  did  injection  tests  on  a 
hydrofoil  with  results  qualitatively  similar  to  Lehman  and  Suessmann. 
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In  the  above  experiments  flow  visualization  had  a low  priority  and  con- 
sequently very  little  was  learned  about  any  pol  yrne  r - induced  niodifiia 
tions  in  the  basic  flow  field. 

In  an  attempt  to  untangle  the  above  observations  and  gain  in 
sight  into  what  the  polymer  was  actually  doing,  Sarpkaya  [ 103]  tested 
six  different  hydrofoils  (NACA  0006,  0012,  0024,  4412,  66j-012  and 
65006)  in  a recirculating  water  tunnel  filled  with  a polymer  solution. 

He  found  that  the  NACA  0012,  0024,  4412,  and  66^-012  foils  showed  a 
decrease  in  lift  and  an  increase  in  drag  whereas  the  NACA  0006  and 
0065  foils  showed  eui  increase  in  lift  and  a decrease  in  drag.  Flow 
visualization  by  means  of  dye  injection  showed  that  the  position  of  lam- 
inar separation  on  the  upper  side  moved  forward  slightly  and  that  the 
separation  line  was  distorted  and  oscillated  back  and  forth.  Sarpkaya 
explains  these  results  as  follows:  the  NACA  0006  and  0065  foils  pos- 
sess a forward  laminar  separation.  Subsequent  transition  on  the  free 
shear  layer  and  reattachment  as  a turbulent  boundary  layer  give  rise 
to  a separation  bubble.  If  this  bubble  is  long,  as  is  the  case  for 
NACA  0006  and  0065,  it  will  increase  the  drag  and  reduce  the  lift. 

For  these  types  of  foils,  Sarpkaya  says  "the  instability  and  the  transi- 
tion to  turbulence  brought  about  by  the  polymer  help  to  replace  the  long 
bubble  by  a turbulent  boundary  layer  and  partly  restore  the  lift  other- 
wise lost  due  to  the  bubble.  Furthermore,  the  relative  decrease  in  the 
apparent  size  of  the  foil  as  well  as  the  friction  reduction  with  polymer, 
in  turbulent  flow  over  the  upper  and  lower  sides  of  the  foil,  results  in 
a Reynolds -number-dependent  drag  reduction.”  Presumably  by 
"apparent  size”  Sarpkaya  is  referring  to  the  foil  thickness. 
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As  for  the  remaining  foils,  Sarpkaya  explains  as  follows:  on 
the  NACA  0012,  0024  and  4412  foils,  which  have  a turbulent  separa- 
tion near  the  trailing  edge,  the  early  polymer-induced  transition  to  a 
turbulent  boundary  layer  causes  a thicker  boundary  layer  at  the  trail- 
ing edge  and  thus  a thicker  wake.  This  thicker  wake,  it  is  claimed, 
will  cause  a reduction  in  lift  and  the  thicker  boundary  layer  will  cause 
an  increase  in  effective  foil  size  and  thus  a larger  drag.  This  specu- 
lation by  Sarpkaya  is  supported  by  the  results  of  Tagori  et  al  [ 104]. 
They  studied  the  flow  past  NACA  00 12  and  66  ^-0 12  foils  in  a water  solution 
of  polyacrylamide  (separan  AP-300)  and  used  air  bubbles  injected  into 
the  freestream  to  visualize  the  flow.  These  results  clearly  show  a 
thickened  turbulent  boundary  layer  when  the  polymer  is  present. 

It  is  clear  from  the  present  experiments  and  the  ones  just  dis  - 
cussed that  the  presence  of  a drag- reducing  polymer  can  greatly  modi- 
fy the  laminar  boundary  layer  including  the  features  of  separation  amd 
trainsition.  Some  indications  of  what  may  happen  for  separation  have 
been  sketched  above.  We  now  review  what  little  information  on  the 
stability  of  a laminar  boundary  layer  is  available  with  a view  of  obtain- 
ing more  ideas  as  to  the  effects  to  be  expected  for  transition. 

IV . 9 Stability  of  Polymer  Flows  — Experimental 

As  would  be  expected*,  a large  fraction  of  the  experimental 
work  on  the  laminar  to  turbulent  transition  in  dilute  polymer  solutions 
has  been  done  in  pipe  flows.  The  majority  of  the  studies  indicate  a 
delay  in  transition  [90].  An  example  of  this  is  given  by  Castro  and 
Squire  [ 105]  who  measured  the  "intermittency"  (ratio  of  length  of  time 
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the  flow  is  turbulent  to  the  length  of  time  the  flow  is  laminar)  of  the 
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flow  through  a pipe.  In  teste  on  a number  of  different  molecular 
weight  polymers,  they  found  the  presence  of  polymer  decreased  both 
the  intermittency  and  the  intensity  of  the  turbulent  flashes.  Giles  and 
Pettit  [106]  used  dye  to  visualize  the  transition  process  in  glass  tubes. 
They  found  very  large  delays  in  the  occurrence  of  transition  and  the 
delay  increased  with  polymer  concentration.  White  and  McEligot  [ 1 07] 
made  studies  of  transition  in  small  tubes  for  mixtures  of  polymers 
with  different  molecular  weights.  They  also  found  delays  in  the 
occurrence  of  transition  with  the  delay  being  influenced  most  by  the 
polymer  with  the  highest  molecular  weight. 

However,  there  are  some  studies  which  found  no  delay , and  in 
a few  cases  an  even  earlier  transition  is  found  than  for  pure  water. 
Virk  et  al  [108]  tested  five  polyethylene  oxides  in  pipe  flows  and  found 
in  the  majority  of  the  tests  no  delay  in  transition  occurred.  Goldstein 
et  al  [109]  used  a laser  doppler  anemometer  to  investigate  turbulent 
and  transitional  pipe  flows.  They  found  that  although  the  presence  of 
the  polymer  did  not  alter  the  longitudinal  turbulent  intensity , an  un- 
steady flow  existed  at  Reynolds  numbers  for  which  the  water  flow  is 
normally  laminar.  More  recently  Patersonand  Abernathy  [ 1 10]  have 
studied  laminar  to  turbulent  transition  in  pipe  flows  of  dilute  solutions 
of  polyethylene  in  which  both  a smooth  and  an  abrupt  inlet  section  were 
used.  For  the  high  disturbance  or  abrupt  inlet  case  the  polymer  was 
found  to  have  no  effect  on  transition.  Whereas  for  the  low  disturbance 
or  smooth  inlet  case  the  polymer  was  found  to  cause  an  earlier  transi- 


tion to  occur . 

Besides  pipe  flows , investigations  into  the  stability  of  flow  between 
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two  concentric  cylinders  with  relative  rotation  have  been  carried  out. 
However,  since  these  flows  are  dominated  by  secondary  fluid  motions 
and  do  not  have  a laminar  to  turbulent  transition,  those  investigations 
will  not  be  discussed  here.  As  for  flows  over  flat  plates  and  discs, 

Hoyt  [90]  has  reviewed  the  available  information.  However  in  these 
investigations  the  main  effort  was  to  measure  drag  reduction  in  a tur 
bulent  boundary  layer  and  few  attempts  were  made  to  see  what  was  hap- 
pening to  the  pre-existing  laminar  boundary  layer. 

The  experimental  findings  reviewed  above  are  certainly  contra- 
dictory and  confusing.  This  confusion  is  mainly  due  to  a lack  of  under- 
standing of  the  basic  mechanism  by  which  the  polymer  interferes  with 
the  original  flow  and  thus  an  incomplete  understanding  of  the  variables 
affecting  the  polymer  behavior.  Another  factor  contributing  to  the  con- 
fusion is  the  precise  determination  of  the  laminar  to  turbulent  transi- 
tion. In  the  works  just  reviewed,  three  methods  were  used:  the 
"break  point"  in  the  pipe  friction  factor  curve,  measurement  of  the  "in- 
termittency"  and  flow  visualization.  Use  of  the  break  point  to  define 
transition  is  misleading  since  both  Paterson  and  Abernathy  [llO]  and 
Goldstein  et  al  [109]  found  unsteady  flow  to  exist  for  Reynolds  numbers 
in  the  laminar  portion  of  the  pipe  friction  factor  curve  for  dilute  poly- 
mer solutions.  Flow  visualization  is  a better  technique  but  it  has  the 
disadvantage  of  being  dependent  on  the  observer.  The  concept  of  inter- 
mittency  appears  to  be  best  because  an  arbitrary  value  can  be  given  to 
the  intermittency  factor  and  applied  to  all  experiments.  Further  it  is 
not  observer  dependent. 

Since  the  experimental  investigations  (into  the  laminar  to  tur- 
bulent transition  in  polymer  flows)  are  so  inconclusive,  it  was  hoped 
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that  the  theoretical  investigations  would  yield  more  consistent  results. 

A brief  review  of  some  of  these  theoretical  works  now  follows. 

IV.  10  Stability  of  Polymer  Flows  — Theoretical 

In  the  theoretical  works  reviewed,  the  approach  to  the  problem 
of  determining  laminar  to  turbulent  transition  is  the  same  as  for 
Newtonian  fluids,  i.e  solution  of  the  Orr-Sommerfeld  equation.  The 
only  modification  to  the  equation  to  account  for  the  presence  of  the 
polymer  is  to  change  the  constitutive  relation  to  include  viscoelastic 
effects . 

Fong  and  Walters  [ill]  used  a model  of  a slightly  viscoelastic 
fluid  to  investigate  the  stability  of  plane  Poiseuille  flow.  Their  solu- 
tions predicted  a destabilization  due  to  the  presence  of  elasticity. 

Since  they  felt  this  result  was  in  contradiction  to  experimental  observa- 
tions, they  repeated  the  analysis  using  another  constitutive  relation 
which  also  included  a variation  in  viscosity  with  the  rate  of  shear  [ll2]. 

‘ Again  their  solutions  predicted  a destabilizing  influence.  They  then  re- 

' turned  to  the  original  constitutive  relation  and  <malyzed  the  stability  of 

a laminar  boundary  layer.  This  time  they  found  the  presence  of  elas- 
ticity to  have  a stabilizing  effect.  Gupta  [ 113]  investigated  the  flow  of 
a viscoelastic  liquid  flowing  down  an  inclined  plane.  He  found  the  pres- 
ence of  elasticity  caused  the  flow  to  be  unstable  at  any  Reynolds  number  1 
■<  Craik  [ll4]  showed  that  this  unsatisfactory  result  was  due  to  an  inade- 

quate fluid  model.  More  recently  Mook  [ll5]  found  elasticity  to  be 
stabilizing  in  combined  Couette- Poiseuille  flow  and  also  in  a flat  plate 
boundary  layer  flow. 


i 


-b9- 


Ln  these  studies  it  has  been  assumed  that  the  staijility  of  two  (li- 
mensional  disturbances  will  control  transition  as  in  the  case  of 
Newtonian  flow  [ll6].  However,  Lockett  [ll7]  suggests  that  the  visco- 
elasticity may  make  it  necessary  to  study  the  full  three  diiTiensionalit> 
of  the  flow.  This  is  substantiated  somewhat  by  Brennen's  [66]  obser- 
vations of  three-dimensional  instabilities  on  fully  developed  cavity  sur- 
faces in  a drag- reducing  flow  past  spheres  and  cylinders.  Kundu  [ll8] 
using  an  energy  method  studied  the  stability  of  two  and  three  dimen- 
sional disturbeuices  in  plane  Couette  flow  of  a second-order  fluid.  He 
found  the  presence  of  elasticity  to  be  stabilizing  to  two  dimensional  dis- 
turbances but  three  dimensional  disturbeinces  could  be  either  stabilized 
or  destabilized  depending  on  the  consteints  in  the  constitutive  relation. 

From  all  this  work  it  seems  then  that  theoretical  calculations 
are  hampered  by  the  lack  of  the  "correct"  constitutive  relation  and  per- 
haps the  need  for  a full  three  dimensional  consideration  of  disturbances. 
A summary  of  both  the  experimental  and  the  theoretical  investigations 
can  at  best  say  that  the  influence  of  a drag- reducing  additive  upon  lam- 
inar to  turbulent  transition  may  be  stabilizing  or  destabilizing  depend- 
ing upon  the  flow  geometry,  type  of  polymer  and  its  concentration. 

IV , 1 1 A Correlation  Between  Disturbance  Amplification  and  G 

It  had  been  hoped  that  a review  of  the  literature  would  indicate 
a practical  modification  to  the  approximate  stability  calculation  method 
of  Jaffe  et  al  [45]  that  would  make  it  applicable  to  polymeric  flows. 
However,  from  the  preceding  review  no  such  simple  modification  is 
yet  possible  with  the  present  theoretical  means.  Yet,  there  is  some 
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advantage  to  calculating  the  added  disturbance  amplification  due  to  the 
presence  of  the  p>olymer  and  see  if  this  can  be  related  to  the  amount  of 
polymer  present.  For  a given  Reynolds  number  the  actual  amplifica- 
tion at  the  position  of  laminar  separation  was  estimated  by  the  method 
of  Jaffe  et  al  [45].  Then  using  the  e criterion  the  extra  amplifica- 
tion required  to  produce  transition  at  the  position  of  separation  was 
calculated.  This  added  amplification  plotted  against  the  value  of  G at 
which  the  laminar  separation  first  disappears  is  shown  in  Fig.  55. 
Comparison  of  these  results  with  those  of  van  der  Meulen  [59]  was  not 
possible  since  he  always  injected  polymer  solution  at  a rate  such  that 
the  laminar  separation  was  eliminated. 

As  CeUi  be  seen  in  Fig.  55,  there  is  a good  deal  of  scatter  but  a 
correlation  does  seem  to  exist.  But  perhaps  it  is  only  fortuitous  since 
the  calculations  based  on  two-dimensional  small  disturbance  theory  are 
almost  certainly  invalid  at  the  Reynolds  numbers  of  the  present  experi- 
ments. Moreover,  it  is  likely  [66]  that  the  disturbances  are  three- 
dimensional  although  this  remains  to  be  confirmed  in  fully  wetted  flow. 

IV . 1 2 Polymer  Mechanism 

The  previous  sections  have  discussed  the  modifications  to  the 
laminar  to  turbulent  transition  produced  by  the  presence  of  polymer, 
but  no  mention  of  a mechanism  to  explain  these  effects  has  been  made. 
Hoyt  [90]  has  reviewed  the  existing  mechanisms  proposed  to  explain 
friction  reduction,  but  these  theories  are  based  mainly  upon  observa- 
tion from  turbulent  pipe  flow  of  polymer  solutions.  First  it  is  not  clear 
that  the  polymer  acts  the  same  in  a laminar  boundary  layer  as  in  a 
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turbulent  one  and  second,  established  turbulent  pipe  flow  is  quite  dif- 
ferent from  the  flow  over  spheres,  cylinders,  etc.  There  are,  never- 
theless, a few  topics  gained  from  experience  in  pipe  flows  that  may  be 
discussed  with  respect  to  the  present  experiments; 

One  of  these  is  that  of  a critical  or  onset  shear  stress.  It  has 
been  noted  in  pipe  flows  of  drag  reducing  polymer  solutions  that  no  ef- 
fect is  observed  until  the  wall  shear  stress  reaches  a certain  critical 
value.  The  maximum  shear  stres ses  on  the  NSRDC  and  hemisphere  nose 
bodies  were  calculated  for  various  values  of  the  Reynolds  number  and 
compared  to  the  critical  value  of  2.  1 dynes /cm  for  Polyox  (WSR  301) 
[90].  In  all  cases  the  calculated  shear  stress  was  much  greater  than 
the  critical  value.  However,  in  the  present  experiments  no  tests  were 
done  for  which  the  maximum  shear  stress  was  less  than  the  critical 
value,  so  no  definitive  conclusion  about  the  significance  of  the  shear 
stress  effect  in  activating  the  polymer  can  be  made  for  these  experi- 
ments . 

A concept,  which  has  arisen  to  try  to  explain  these  various 
polymer  effects  in  laminar  and  turbulent  flow,  is  the  idea  of  molecular 
stretching.  Tulin  [H9]  and  Lumley  [ 120,  121  ] tried  to  explain  the 
turbulent  drag- reduction  by  suggesting  that  if  the  polymer  molecules 
were  stretched  from  their  naturally  coiled  condition  they  could  become 
large  enough  to  interfere  with  turbulent  eddies.  Brennen  [66]  has  used 
this  idea  in  a somewhat  different  way  to  forward  an  explanation  for  the 
observed  destabilization  of  the  boundary  layer. 

Brennen  proposed  that  the  highly  favorable  pressure  gradient  on 
the  forward  portions  of  spheres  and  cylinders  stretches  the  long  chain 
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molecules  out  of  their  naturally  coiled  state.  Once  past  the  minimum 
pressure  point  an  adverse  pressure  gradient  exists  and  the  polymer 
molecule  is  allowed  to  relax.  As  the  molecule  relaxes,  it  would  put 
energy  into  the  boundary  layer  and  if  the  boundary  layer  could  be  suf- 
ficiently energized,  transition  could  occur  upstream  of  the  position  at 
which  it  normally  occurs.  Essentially  then  the  molecule  is  considered 
as  a spring  which  is  first  stretched  and  then  released  as  it  passes  over 
the  body.  If  this  analogy  is  applicable,  then  it  would  be  expected  that 
the  molecule  would  oscillate  when  released  at  a frequency  equal  to  the 
reciprocal  molecular  relaxation  time.  For  Polyox  (WSR  301)  this  is 
1000  Hz  [so].  Now  it  is  known  from  solutions  of  the  Orr-Somme rfeld 
equation  that  a certain  rcuige  of  disturbance  frequencies  is  amplified  in 
the  boundary  layer  to  a much  greater  extent  than  other  frequencies.  In 
particular  there  is  one  frequency  (called  the  critical  frequency)  at 
which  maximum  amplification  occurs.  Using  the  approximate  calcula- 
tion technique  [4S]  the  critical  frequencies  for  the  present  bodies  were 
calculated  for  several  different  Reynolds  numbers  and  are  presented  in 
the  table  below. 


TA  BLE  IV.  1 Critical  Boundarv 


luencies 


(at  separation) 
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Comparison  of  these  critical  boundary  layer  frequencies  with  the  reso- 
nant frequency  of  the  molecule  indicates  that  the  boundary  layer  would 
be  particularly  responsive  to  disturbances  imposed  on  it  by  oscillations 
of  the  polymer  molecule.  The  NSRDC  body  would  seem  especially  sus- 
ceptible to  polymer  influence  since  its  critical  frequencies  are  very 
near  the  polymer  frequency  and  also  the  boundary  layer  is  highly  am- 
plified to  begin  with.  This  was  substantiated  by  experiment  in  that  for 
a given  Reynolds  number  and  concentration  a lower  injection  rate  was 
required  for  the  NSRDC  body  than  for  the  hemisphere  nose  body  to 
eliminate  the  separation,  eg.  see  Figs.  50,  51.  Thus  the  "stretching" 
hypothesis  seems  to  be  a plausible  mechanism  by  which  the  laminar 
boundary  layer  could  be  destabilized.  However,  in  view  of  the  lack  of 
direct  evidence  that  polymer  stretching  actually  does  occur  in  the  pres- 
ent flows,  this  mechanism  must  still  be  regarded  as  being  a hypotheti- 
cal one. 

In  summary  then  there  is  substantial  experimental  evidence 
that  the  presence  of  drag-reducing  polymers  causes  a destabilization 
of  the  laminar  boundary  layer.  Theoretical  evidence  is  less  clear  and 
predicts  either  a stabilization  or  destabilization  of  the  laminar  boundary 
layer  can  occur. 
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Chapter  V 


PRESENTATION  AND  DISCUSSION  OP'  CAVITATION  RESULTS 


V.  1 Freestreaim  Turbulence  Level  and  Cavitation  Inception  — Presen- 
tation of  Results 

The  cavitation  inception  characteristics  of  the  hemisphere  nose 
body  were  studied  in  the  HSWT  and  the  LTWT  at  its  lowest  turbulence 
level.  The  NSRDC  body  was  tested  in  the  LTWT  only  and  it  was  tested 
at  two  freestream  velocities  and  two  freestream  disturbance  levels. 

V.  1.  1 Visual  Observations  — Hemisphere  Nose  Body 

The  cavitation  inception  sequence  upon  the  hemisphere  nose 
body  was  the  same  in  both  the  HSWT  and  the  LTWT  and  has  already 
been  described  in  detail  in  section  1.7.  A ^'  otograph  illustrating  the 
appearance  of  developed  cavitation  on  this  body  can  be  seen  in  Fig.  58a. 

V.  1 . 2 Visual  Observations  — NSRDC  Body 

At  the  low  turbulence  level  (0.  05  percent),  band  type  cavitation 
inception  always  occurred  with  no  initial  appearance  of  steady  bubble 
type  cavitation.  At  the  high  turbulence  level  (0.65  percent)  cavitation 
was  of  the  same  type  and  occurred  at  the  same  value  of  the  cavitation 
number.  However,  at  the  low  disturbance  level,  the  subsequent  devel- 
oped cavitation  was  steady  whereas  the  developed  cavitation  at  the 
high  disturbance  level  was  unsteady.  At  the  0.65  percent  disturbance 
level  the  separation  line  became  very  irregular  and  erratic  and 
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occasionally  portions  of  the  body  would  become  rewetted.  This  differ- 
ence in  behavior  between  the  low  and  high  turbulence  levels  is  denion- 
strated  in  the  photographs  of  Fig.  56, 

V.  2 LTWT  Limitations 

Before  discussing  the  cavitation  experiments,  we  need  to  explain 
their  limited  test  range.  Only  a small  velocity  range  is  available  since 
the  LTWT  has  a maximum  test  section  velocity  of  only  25  feet  per  sec- 
ond and  the  lowest  velocity  at  which  bodies  of  the  present  shape  can  be 
made  to  cavitate  is  about  18  feet  per  second.  Hence,  tests  we  re  carried 
out  at  two  velocities  only.  There  is  another  limitation  . Although  five 
freestream  turbulence  levels  are  available,  only  two  could  be  used  — 
the  0.05  and  0.65  percent  levels.  This  is  so  because  turbulence  grids 
located  at  the  entrance  to  the  test  section  would  cavitate  before  the  test 
model  and  thus  three  of  the  turbulence  levels  were  eliminated  from  consid- 
eration. With  these  restrictions  then,  the  experiments  conducted  in  the 
LTWT  to  determine  the  influence  of  freestream  turbulence  level  upon 
cavitation  inception  were:  for  the  hemisphere  nose,  Pej^=  3.0  x 10 

and  3.8  X 10^  at  0.05  percent  turbulence  level:  and  for  the  NSRDC 

5 5 

body,  RCjj  = 3.3  x 10  and  3.8  x 10  at  0.05  percent  and 
5 

Rep  = 3.3  X 10  at  0.65  percent  turbulence  level. 


parameter.  Included  in  Fig.  57  for  comparison  are  cavitation  data  ob- 
tained on  the  same  model  in  the  NSRDC  facility.  The  cavitation  results 
for  the  hemisphere  nose  body  in  both  the  LTWT  and  the  HSWT  have 
been  included  in  Fig.  60.  As  mentioned  earlier,  the  high  disturliance 
level  caused  no  change  in  the  inception  value  on  the  NSRDC  body  but  it 
did  produce  a draunatic  change  in  the  appearance  of  the  developed  cavi- 
tation . Schlieren  visualization  of  the  fully  wetted  flow  at  the  0.65  per- 
cent disturbance  level  showed  only  a slight  reduction  in  size  of  the  lam- 
inar separation  region.  Comparison  of  the  scatter  in  the  measured 
laminar  separation  dimensions  between  the  low  and  the  high  disturbance 
levels  show  only  slightly  higher  scatter  at  the  higher  level  — thus  sug- 
gesting that  the  laminar  separation  is  only  slightly  more  unsteady  at  the 
high  turbulence  level  than  at  the  low  turbulence  level.  However,  the 
cavitation  indicates  that  the  freestream  disturbances  influence  the  vis- 
cous flow  to  a greater  extent  than  would  be  predicted  from  only  the 
above  schlieren  observations.  Or,  perhaps  the  explanation  is  that  the 
developed  cavity  is  more  sensitive  to  freestream  disturbances  than  the 
fuily  wetted  shear  layer.  If  this  were  so,  the  value  of  the  cavitation 
number  at  inception  would  not  change,  but  the  appearance  of  the  devel- 
oped cavity  would. 

V.  4 Comparison  with  Other  Investigations 

In  general  the  main  purpose  of  the  investigation  of  freestream 
turbulence  level  upon  cavitation  occurrence  and  type  was  to  determine 
if  it  could  be  a contributing  factor  in  the  differences  in  cavitation  re- 
sults on  identical  bodies  tested  in  different  facilities.  In  particular  we 


-77- 


wanted  to  know  if  the  observed  differences  in  cavitation  on  the  same 
NSRDC  body  when  tested  in  the  NSRDC  facility  versus  the  HSW'J  at 
CIT  could  be  explained  by  different  turbulence  levels  in  these  facilities 

It  had  been  observed  [l22]  that  in  the  NSRDC  facility  the  critical 
velocity  for  this  body  was  9.3  feet  per  second.  All  cavitation  tests 
were  carried  out  at  supercritical  velocities  and  cavitation  was  always 
of  the  traveling  bubble  type.  In  the  HSWT  the  same  body  was  found  to 
have  a critical  velocity  of  30  feet  per  second  and  the  type  of  cavitation 
depended  on  velocity  (see  Section  1.7).  It  was  recognized  that  the  dif- 
ferent viscous  flows  on  the  two  bodies  was  at  least  partially  responsible 
for  the  differences  in  cavitation.  However,  it  was  not  known  why  the 
critical  velocity  in  the  NSRDC  facility  was  so  low  compared  to  the 
HSWT.  It  was  suspected  [44]  that  the  freestream  turbulence  level 
might  be  responsible  and  thus  the  tests  in  the  LTWT  were  carried  out. 

At  the  lowest  turbulence  level  (0.05  percent)  in  the  LTWT  the 
leiminar  separation  was  found  to  be  still  present  at  a velocity  of  24  feet 
per  second  and  thus  supports  the  result  in  the  HSWT  (measured  turbu- 
lence level  is  0.25  percent'll).  And  also,  band  type  inception  always 
occurred.  To  reduce  the  critical  velocity  to  10  feet  per  second  in  the 
LTWT  required  a turbulence  level  of  3.6  percent.  As  can  be  seen  from 
Table  I.  1,  this  is  an  exceptionally  high  turbulence  level  in  a water  tun- 
nel test  section.  Initially  it  was  thought  unlikely  that  the  disturbance 
level  in  the  NSRDC  facility  is  this  high.  However,  after  inspecting  a 
drawing  of  the  NSRDC  facility  (Fig.  2.3,  I'Inapp  et  al  [26]),  this  result 
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does  not  seem  so  unlikely.  The  bulk  of  the  circuit  is  18-inch  diameter 
pipe  with  turning  vanes  in  only  one  corner;  it  has  a small  (6,25:1)  con- 
traction ratio  and  only  one  honeycomb^  located  immediately  after  the 
tunnel  pump.  This  type  of  water  tunnel  would  be  expected  to  have  a 
high  turbulence  level  and  in  view  of  the  changes  in  cavitation  made  by  a 
0.65  percent  disturbance  level,  it  seems  likely  that  the  freestream 
turbulence  is  responsible  for  the  low  critical  velocity  and  thus  partially 
responsible  for  the  differences  in  cavitation  results.  However,  it  is 
believed  differences  in  freestream  nuclei  populations  between  these 
facilities  must  also  be  included  in  a discussion  of  this  topic  and  will  be 
done  so  in  the  next  chapter. 

V . 5 Polymer  Additives  and  Cavitation  Inception 

V.  5.  1 Visual  Observations  - Hemisphere  Nose  Body 

As  in  the  fully  wetted  flow  visualization  experiments,  the  effect 
of  the  polymer  caused  a significant  chsmge  in  the  flow  and  here  in  par- 
ticular the  appearance  of  the  cavitation  on  the  hemisphere  nose  body 
was  strongly  dependent  on  the  amount  of  polymer  present  in  the  bound- 
ary layer.  For  a fixed  concentration  and  freestream  velocity  the  fol- 
lowing changes  were  observed  to  take  place  as  the  injection  rate  was 
increased:  At  zero  injection  rate,  incipient  band  type  cavitation  always 
occurred.  At  injection  rates  less  than  the  critical  value  at  which  the 
separation  would  disappear,  incipient  band  type  inception  would  still 
occur  but  as  can  be  seen  in  Fig.  58b  the  surface  of  the  developed  cav- 
ity has  a definite  wave  structure  and  the  separation  line  has  become 
very  irregular.  Inspection  of  schlieren  photographs  of  the  fully  wetted 
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flow  at  this  injection  rate  showed  a reduced  laminar  separation  re^>ion 
with  transition  on  the  free  shear  layer  having  moved  upstream  from 
the  no-injection  case.  With  a further  increase  in  the  injection  rate  to 
values  near  the  critical  rate,  different  types  of  cavitation  were  ob- 
served depending  upon  the  facility.  In  the  HSWT  band  type  inception 

I would  occur  in  patches  with  irregular  separation  lines  and  surfaces  as 

( 

j is  shown  in  Figs.  58c, d.  Furthermore  this  type  of  cavitation  would 

i 

i occur  in  an  intermittent  fashion.  At  even  higher  injection  rates  above 

the  critical  one,  a similar  behavior  was  observed  but  with  the  flow  al- 
tering between  fully  wetted  and  patchy  band  cavitation  more  rapidly.  A 
decrease  in  cavitation  number  at  this  injection  condition  would  make 
this  type  of  activity  more  violent,  but  no  steady  attached  cavitation 
could  be  obtained.  At  this  same  injection  rate,  the  fully  wetted  obser- 
vations showed  the  laminar  separation  to  be  eliminated  for  the  most 
part  with  only  an  occasional  reappearance.  That  is,  the  fully  wetted 
flow  is  almost  always  turbulent  in  this  region  of  interest.  If  then  at 
this  stage  the  injection  rate  was  suddenly  reduced  to  zero,  a lar^e 

t 

» steady  cavity  would  quickly  form  on  the  body.  In  the  LTWT  the  same 

sequence  of  cavitation  development  would  take  place.  However,  in  this 
facility  traveling  bubble  type  and  band  type  cavitation  would  occur  simul- 
taneously when  polymer  was  injected  into  the  boundary  layer  unlike  the 
HSWT  where  rm  bubble  type  cavitation  was  observed.  We  will  return  to 
this  importamt  difference  in  Chapter  VI. 

V.5.2  Visual  Observations  — NSE  DC  Body 


The  NSRDC  body  was  tested  in  the  LTWT  only  and  it  too  was 
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observed  to  go  through  a sequence  of  cavitation  development  similar 
to  that  described  above  for  the  hemisphere  nose  body  in  the  LTW'l  ; 
namely,  that  the  injection  of  polymer  at  subcritical  rates  changed  the 
original  bcUid  type  cavitation  into  alternating  band  and  bubble  type. 

And  at  supercritical  injection  rates  this  intermittency  became  more 
rapid  and  no  attached  cavitation  could  be  sustained.  Examples  of  these 
types  of  cavitation  on  the  NSRDC  body  are  shown  in  Fig.  59  wherein  it 
may  be  seen  that  the  polymer  injection  has  an  equally  dramatic  effect 
as  on  the  hemisphere  nose  body.  This  is  demonstrated  particularly 
well  in  Fig.  59d  where  only  one  cavitation  bubble  can  be  seen  at  the 
extremely  low  cavitation  number  of  0.34! 

V.  6 Discussion 

Results  of  the  cavitation  tests  with  polymer  injection  are  sum- 
marized in  Figs.  60  and  61  where  the  values  of  the  inception  cavitation 
number  are  plotted  versus  Reynolds  number  for  various  values  of  the 
parameter  G.  We  pause  here  to  discuss  the  parameter  G:  In  Fig.  62 
the  inception  cavitation  number  has  been  plotted  against  the  polymer 
injection  rate  for  two  concentrations,  but  in  Fig.  63  the  same  data  have 
been  replotted  against  the  factor  G instead  of  the  injection  rate.  It 
may  be  seen  there  that  both  sets  of  data  have  collapsed  onto  one  curve. 
A similar  result  was  also  found  for  the  laminar  separation  dimensions 
when  plotted  versus  G.  We  may  recall  from  the  definition  of  G 
(Chapter  IV)  that  for  a fixed  freestream  velocity  and  body  size,  G 
becomes  proportional  to  the  mass  flowrate  of  the  polymer  in  the  bound- 
ary layer.  So  the  correlation  of  the  cavitation  index  with  G implies 
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that  the  polymer  "effectiveness"  is  proportional  only  to  the*  amounl 
present  within  the  boundary  layer.  Evidently  from  this  finding;,  the 
amount  present  can  either  be  controlled  by  the  concentration  or  the*  in- 
jection rate,  and  hence  it  does  appear  that  only  G is  significant  for 
the  injection  effect. 

Returning  to  the  results  of  Figs.  60  and  61  then,  it  can  be  seen 
that  increasing  amounts  of  polymer  in  the  boundary  layer  produce  an 
increasing  suppression  in  cavitation  inception.  There  is  a limit,  how- 
ever, beyond  which  no  further  increase  in  cavitation  suppression  oc- 

. -6 

curs;  in  the  present  case  on  the  hemisphere  nose  body  G = 7 x 10 
and  by  observation  of  the  sclilieren  photographs,  this  limit  concurs  with 
the  removal  of  the  laminar  separation.  Van  der  Meulen  [59]  injected 

a 500  wppm  salt-polymer  solution  in  the  range  40x10"^  < G < 72x  10"^ 

\ 

and  for  all  of  these  conditions  and  test  velocities  the  laminar  separa- 
tion was  found  to  have  been  removed.  These  and  the  present  results 
are  compared  in  Fig.  64  where  the  maximum  percent  reduction  in  cavi- 
tation index  has  been  plotted  versus  the  Reynolds  number.  The  agree- 
ment is  believed  to  be  reasonably  good  for  experiments  of  this  type. 
Values  from  Baker  et  al  [52],  Holl  et  al  [53]  and  Ellis  et  al  [50j  who 
did  polymer  "ocean"  tests  on  hemisphere  nose  bodies  are  also  shown. 
However,  the  information  from  their  reports  is  limited  and  all  that  can 
be  said  is  that  they  give  values  approximately  the  same  as  those  noted 
in  the  present  study. 

During  their  cavitation  tests  Baker  et  al  ( 52]  and  Holl  et  al  [53] 
noted  a change  in  the  appearance  of  the  developed  cavitation.  F rom 
photographic  observations  of  these  changes  Holl,  Baker  et  al  speculated 
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that  the  cavitation  attenuation  was  tluc  to  a "flow  n-or icnlal  ion  in  Itic 

region  of  the  laminar  separation  bubble."  They  further  speeui  ,il  cf! 

[54]  that  the  amount  of  attenuation  would  depend  oi»  Itie  r-ilio  of  lin- 

molecular  relaxation  time(T)  to  a flow  time  scale  based  on  tin  fn  e- 

stream  velocity  and  the  boundary  layer  displacement  thickness  at  si'(j- 

aration  (6*)  i.e.  (T  • U lb  ' ).  This  idea  is  equivalent  to  tiiat  cjf  a 
s 00  s 

critical  shear  stress  discussed  in  Chapter  IV.  The  work  of  van  det 
Meulen  [59]  and  the  present  investigation  verify  the  first  speculation 
about  flow  reorientation.  However,  no  correlation  of  the  present  re- 
sults with  the  parameter  T • U / b^~-  could  be  obtained.  The  absence 

00 

of  a suppression  effect  on  an  eight  inch  diameter  nose  model  observed 
by  Arndt  et  al  [54]  is  believed  to  be  due  to  the  laminar  separation  be- 
ing very  much  reduced  in  size  or  eliminated  by  natural  transition. 

The  Reynolds  number  of  the  test  was  4 x 10^  [ 54]  and  the  freeslream 
turbulence  level  has  been  reported  to  be  about  0.8  percent  [49].  The 
critical  Reynolds  number  for  this  type  of  body  has  been  estimated  to 
be  5 X 10^  at  a zero  freestream  turbulence  level  [4l]. 

V.  7 Comparison  with  Boundary  Layer  Trip  Experiments 

The  suppression  of  cavitation  upon  these  bodies  is  a result  of 
the  elimination  of  the  laminar  separation  by  the  polymer.  It  is  worth- 
while to  compare  the  present  results  with  those  in  which  the  separa- 
tion is  eliminated  by  another  method.  Arakeri  and  Acosta  [41  ] did  a 
series  of  tests  with  a hemisphere  nose  body  and  the  ITT C body  using 


boundary  layer  trips  to  reduce  the  critical  Reynolds  number.  They 
found,  with  the  trip  present  and  at  velocities  above  the  new  critical 
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velocity,  that  the  occurrence  of  cavitation  was  si^;nifican(ly  siij)pt  t-ssed 
and  that  at  higher  velocities  the  tunnel  would  choke  before  the  body 
could  be  made  to  cavitate.  The  present  polymer  tests  show  a very 
large  effect  on  inhibiting  cavitation  but  not  as  dramatic  as  the  tripped 
tests,  particularly  at  the  higher  velocities.  Since  in  each  case 
the  schlieren  photographs  showed  the  laminar  separation  had  been  elim- 
inated, some  other  factor  must  be  the  controlling  influence  in  the  pres- 
ent experiments.  One  possibility  is  that  the  injected  solution  itself 
supplied  nuclei  to  initiate  cavitation.  Even  though  the  poly  me  r solution 
is  mixed  with  deaerated  tunnel  water,  it  is  under  an  air  pressure  of 
80  psi  in  the  holding  tank  and  may  therefore,  eventually  become  satu- 
rated at  a higher  pressure.  However,  if  this  were  the  case , some  t rav- 
eling bubble  type  cavitation  would  be  expected  to  occur.  The  photo- 
gra  'ig.  58  show  that  no  bubble  type  cavitation  occurs  and  in  fact 

th  ance  of  the  cavitation  is  much  the  same  as  that  observed  on 

the  tripped  bodies  [4l]. 

Another  explanation  is  that  the  body,  when  installed  in  the  HSWT , 
has  a slight  positive  angle  of  attack.  This  angle  of  attack  would  result 
in  a somewhat  uneven  distribution  of  the  polymer  tending  to  send  more 
over  the  top  portion  of  the  body.  From  the  photographs  of  the  physical 
appearance  of  the  cavitation  on  the  hemisphere  nose  body  presented  in 
Fig.  58  it  can  be  seen  that  the  cavitation  first  disappears  on  the  top  of 
the  model.  This  trend  was  repeatable  and  supports  the  above  specula- 
tion.  Further  support  is  provided  by  Arakeri  and  Acosta's  [43]  obser- 
vations of  the  viscous  flow  on  a 1.5  calibre  ogive  body  in  the  HSWT. 

Schlieren  flow  visualization  showed  a laminar  separation  to  be  present 
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on  top  of  the  body  at  a velocity  of  ()5  feet  per  second  wtu-reas  on  the 
bottom  side  no  separation  was  visibh'  — thus  suf'^c-sting  a slif.’h(  p<isi 
tive  angle  of  att»ck. 

Up  to  this  point  the  influence  of  viscous  effects  upon  the  type 
and  occurrence  of  cavitation  have  been  emphasized.  For  a complete 
understanding  of  the  inception  process,  the  nature,  number  and  size 
distributions  of  cavitation  nuclei  and  how  they  interact  with  the  viscous 
flow  past  the  body  must  be  known.  The  next  chapter  presents  some 
measurements  of  nuclei  populations  obtained  in  the  LTWT  and  discus- 
ses their  interaction  with  the  flow  past  the  present  bodies. 
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FREESTRPIAM  NUCLEI  AND  CAVITATION  INCEF'l  ION 


VI.  1 Presentation  of  Result s 

Some  typical  results  tjf  nuclei  population  measurements  made 
during  cavitation  tests  in  the  LTWT  are  presented  in  Fig.  65.  Included 
there  for  comparison  are  nuclei  distributions  measured  in  other  lacili- 
ties  by  different  investigators  The  rmclei  histograms  from  these  in- 
vestigations (including  tlu-  [present  one)  have  been  reduced  to  the  num- 
ber density  function,  N(R),  by  the  following  approximation: 


R, 


+ R 
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2\ 

= 


number  of  nuclei  with  radii  between  Rj  and  R^  per  unit  volume 


(R^  - R^: 


Estimates  of  the  slopes  of  the  lines  through  the  data  have  been 
made  and  are  given  in  Table  VI.  1 shown  on  the  following  page. 


VI.  2 Comparison  of  Nuclei  Distributions 

As  can  be  seen  in  Table  VI.  1,  most  of  the  distributions  have 
slopes  between  -3.  0 and  -4.  0 over  at  least  some  portion  of  the  distri- 
bution. That  the  slope  varies  so  little  is  rather  surprising  in  view  of 
the  many  different  sets  of  environmental  conditions  (ambient  pressure, 
air  content,  velocity)  existing  in  these  facilities  when  the  measure- 
ments were  made.  For  example,  Gavrilov  [1^3]  measured  the  nuclei 
distribution  in  quiescent  tap  water  whereas  in  the  present  study  in  the 
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LTWT  the  water  had  a velocity  of  approximately  20  feel  per  second,  (larl 
an  air  content  of  about  7 ppm  and  experienced  an  ambient  pressun  ol 
about  2 psia.  The  relati'  e insensitivity  of  the  sloiie  to  these:  elramatii 
changes  in  the  environment  implies  that,  although  the  nui  lei  [lopulation 
can  change  by  several  orders  of  magnitude,  ihc  distribution  of  the  nu- 
clei sizes  remains  constant!  Whetlier  or  not  this  eonstan'  y of  slope-  is 
always  true  or  for  what  range  of  nuclei  sizes  it  is  valid  will  reejuirc 
many  more  data  than  are  presently  available. 


VI.  3 Freestream  Nuclei  and  Cavitation  Inception  Tests  in  the  LTWT 

Aside  from  Keller's  work  [12]  no  systematic  study  of  the  effects 
of  these  different  nuclei  distributions  on  cavitation  inception  has  been 
done.  And,  although  the  main  purpose  of  the  present  nuclei  counting 
was  to  obtain  an  accurate  estimate  of  nuclei  distributions  in  the  L.TWT, 
a number  of  interesting  observations  made  during  cavitation  tests  on 
both  bodies  prompted  a few  experiments  in  which  it  was  deliberately 
attempted  to  modify  the  nuclei  distribution.  First,  these  observations 
will  be  described  and  then  the  results  of  some  experiments  designed  to 
explain  these  observations  will  be  discussed. 

As  will  be  recalled  from  Chapter  III,  the  LTWT  has  no  resorber. 
This  necessitates  collecting  cavitation  data  before  the  freestream  bub- 
bles created  at  the  pump  reach  the  test  section.  On  a number  of  occa- 
sions the  cavitation  on  the  body  was  deliberately  maintained  and  the 
freestream  gas  bubbles  from  the  pump  allowed  to  pass  through  the  test 
section.  As  the  number  of  free  gas  bubbles  increased,  the  initially 
occurring  band  t-ype  cavitation  was  gradually  destroyed  and  replaced  by 
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traveling  bubble  type  cavitation.  If  tlie  tunnel  pn  ssuri-  wa.s  lowi  rrd 

so  that  the  body  does  not  cavitate  but  the  pump  still  does  and 

again  the  freestream  bubble  population  allowed  to  increasi-  with 

time,  the  body  would  eventually  cavitate  - inception  then  always 

being  of  the  traveling  bubble  type.  Also,  the  models  could  be 

made  to  cavitate  at  almost  any  value  of  the  cavitation  number  (but 

not  above  O = -c  ) by  merely  allowing  the  number  of  freestream 
^min 

bubbles  to  increase. 

To  try  to  determine  the  reasons  for  the  above  behavior,  the 
effect  of  many  freestream  bubbles  upon  the  fully  wetted  flow  was  inves- 
tigated by  the  schlieren  technique.  Also,  nuclei  populations  were  meas- 
ured when  band  type  inception  occurred  and  when  bubble  type  inception 
occurred  (deliberately  promoted  by  the  technique  described  in  the  pre- 
vious paragraph)  and  compared.  The  influence  of  freestream  nuclei  on 
the  fully  wetted  flow  was  accomplished  by  circulating  water,  which  was 
saturated  with  air,  at  a slightly  reduced  tunnel  pressure  and  using  the 
schlieren  technique  to  visualize  the  flow.  The  results  of  the  schli  en 
study  are  presented  in  Fig.  66  and  the  measured  nuclei  distributions 
are  given  in  Fig.  67. 


VI.  4 Discussion  of  Nuclei  Tests  in  the  LTWT 

Steady  light  source  schlieren  observation  showed  that  as  the 
number  of  freestream  nuclei  increased,  the  laminar  separation  on  the 
hemisphere  nose  body  became  very  unsteady  and  was  finally  eliminated. 
The  schlieren  photographs  in  Fig.  66  were  taken  at  random  as  the  nu- 
clei population  increased  emd  clearly  show  the  effect  of  the  bubbles  on 


the  laminar  separation.  The  mechanism  by  which  these  gas  bubbles 
effect  this  change  is  unknown.  It  has  been  speculated  that  the  presence 
of  the  freestream  bubbles  could  produce  a high  freestream  disturl>ance 
level  and  thus  promote  an  early  laminar  to  turbulent  transition  on  ttie 
body.  However,  in  light  of  the  results  of  llie  turbulence  level  tests  on 
the  hemisphere  nose  body  pre.sented  in  Chapter  IV  this  seems  an  un- 
likely mechanism.  In  any  event,  whatever  the  mechanism  may  be,  the 
observed  elimination  of  the  laminar  separation  by  the  freestream  gas 
bubbles  helps  explain  the  observations  described  in  the  previous  section. 

In  the  case  where  band  type  inception  is  already  present,  the 
increasing  freestre.im  bubble  population  removes  the  laminar  separa- 
tion and  thus  causes  a change  from  band  type  cavitation  to  well  devel- 
oped traveling  bubble  type  cavitation.  For  the  case  where  the  tunnel 
pressure  was  held  at  an  initially  non-cavitating  level  for  the  body,  the 
increasing  number  of  freestream  nuclei  again  remove  the  laminar  sep- 
aration and  thus  traveling  bubble  type  inception  occurs  without  any  pre- 
vious appearance  of  band  type  cavitation.  In  the  latter  case  nuclei  dis- 
tributions were  obtained  at  the  instant  inception  occurred  and  two  of 
these  distributions  (a^  = 0.  58,  0.  73)  are  compared  in  Fig.  67  with  two 
distributions  (a^  - 0.  44)  obtained  when  band  type  inception  occurred, 
i.  e.  with  the  laminar  separation  present.  As  can  be  seen  in  this  figure, 
for  nuclei  w’th  radii  less  than  100  microns  all  the  distributions  are  es- 
sentially the  same  whereas  for  nuclei  greater  than  100  microns  in  ra- 
dius the  bubble  type  inception  distributions  have  many  more  nuclei  than 
the  band  type  inception  ones.  Thus  in  those  facilities  in  which  the  lam- 
inar separation  has  been  removed  either  by  many  freestream  bubbles 
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or  a high  freestream  turbulence  level,  traveling  bubble  type  inception 
will  occur  and  the  value  of  the  cavitation  number  at  which  it  occurs  will 
depend  on  the  nuclei  population.  To  see  if  these  ideas  are  consistent 
with  experiment  consider  the  results  of  the  cavitation  tests  on  the 
NSRDC  body  in  the  NSRDC  facility  and  the  LTWT.  (In  these  tests  the 
same  body  was  tested  in  facilities  for  which  nuclei  distributions  are 
known.  ) 

In  the  NSRDC  facility  at  all  test  velocities  the  laminar  separa- 
tion was  not  present  and  traveling  bubble  type  inception  occurred  at 
vadues  of  the  cavitation  index  between  0.  6 and  0.  8.  In  the  LTWT,  when 
bubble  type  inception  was  deliberately  promoted,  the  laminar  separa- 
tion was  also  absent  and  inception  took  place  at  values  of  the  cavitation 
index  between  0.  58  and  0.  73.  Using  the  static  stability  analysis  for  a 
spherical  gas  bubble,  the  critical  bubble  size  was  calculated  for  the 
given  values  of  the  cavitation  number  and  freestream  velocity.  Then, 
using  the  measured  nuclei  distributions,  the  number  of  "cavitatable" 
nuclei  per  cubic  centimeter  was  estimated.  The  results  of  these  calcu- 
lations are  presented  in  the  Table  VI.  Z shown  on  the  following  page. 

The  nuclei  distributions  for  the  NSRDC  tests  were  also  obtained  by 
holography  [12]  and  an  example  of  one  of  the  distributions  is  given  in 
Fig.  67. 

It  can  be  seen  that  with  the  exception  of  the  one  test  in  the 
LTWT  (a^  = 0.  58)  that  good  agreement  (for  these  types  of  experiments) 
is  obtained  between  the  two  facilities.  If  we  now  follow  Johnson  and 
Hsieh's  example  [39]  and  assume  that  only  nuclei  within  ± 0 01  body 
radii  of  the  body  centerline  are  possible  candidates  for  cavitation  sites. 
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the  number  of  cavitation  events  per  second  can  be  calculated  and  the  re- 
sults of  these  calculations  are  also  given  in  Table  VI.  2,  Pete r son  rneas - 
ured  the  actual  number  of  events  for  the  NSRDC  tests  by  an  acoustic 
method  and  chose  the  value  of  one  event  per  second  to  be  inception  .since 
at  this  rate  good  agreement  was  obtained  with  the  "visual  call"  of  in<  op- 
tion. In  the  LiTWT  inception  was  determined  visually  and  no  measure- 
ments of  event  rates  were  made. 

Comparing  the  calculated  event  rate  between  the  two  facilities, 
it  can  be  seen  that  again  good  agreement  is  obtained.  This  implies  for 
one  thing  that  the  different  experimenters  give  the  same  "visual  call"  of 
inception  and  this  is  good  to  know.  Comparing  the  calculated  event 
rates  to  measured  rates  it  can  be  seen  that  the  calculated  rate  is  in  gen- 
eral higher.  This  is  attributed  to  the  fact  that  not  all  the  calculated 
"cavitatable"  nuclei  actually  cavitate.  This  relates  back  to  the  question 
of  what  the  nuclei  actually  are  and  the  relative  importance  of  the  differ- 
ent types  as  cavitation  nuclei  (e.  g.  particulates  versus  gas  bubbles). 

The  results  of  these  experiments  suggest,  in  those  cases  in 
which  viscous  effects  do  not  control  Inception,  traveling  bubble  type 
cavitation  will  occur  and  the  value  of  the  cavitation  index  at  inception 
will  depend  on  the  nuclei  distribution  and  the  critical  nuclei  size. 

It  has  been  proposed  earlier  by  Silberman  and  Schiebe  [88,  129. 
130,  131]  that  by  measuring  the  number  of  cavitation  events  per  unit 
time  on  a standard  body  as  a function  of  the  cavitation  number  the  cavi- 
tation "susceptibility"  of  the  facility  water  could  be  determined.  This 
proposal  is  the  reverse  of  the  method  used  in  the  present  investigation 
in  which  the  nuclei  population  was  first  determined  and  then  the 
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cavLtation  event  rate  inferred.  By  using  Silberman's  te<  hniquc  a dis- 
tribution is  obtained  of  only  those  nuclei  which  actually  "cavitate"  in- 
stead of  the  entire  nuclei  spectrum  and  thus  presurtiably  a i>etter  corre- 
lation between  calculated  and  measured  event  rates  would  be  aciuevi-d 
The  results  of  the  present  experiments  and  those  of  Peterson  ( 1 j ti-nd 
to  support  Silberman  et  al's  proposal  provided  ol  course  that  the  stan- 
dard body  have  no  laminar  separation  and  that  inception  is  not  cont  rolled 
by  the  laminar  to  turbulent  transition. 

VI.  5 Cavitation  Nuclei  and  Boundary  Layer  Interaction 

In  the  preceding  section  the  nuclei  distributions  were  discussed 
with  reference  to  traveling  bubble  type  cavitation  inception  on  bodies 
which  do  not  have  a laminar  separation.  As  has  been  mentioned  before, 
in  cases  where  the  laminar  separation  is  present  a good  correlation 
can  be  made  between  the  incipient  cavitation  number  and  the  pressure 

coefficient  at  the  position  of  the  separation  [41  ].  For  the  case  of  an 
attached  laminar  to  turbulent  transition,  Arakeri  [42]  found  good 
agreement  between  the  pressure  coefficient  at  the  calculated  position 
of  transition  and  the  desinent  cavitation  data  of  Parkin  and  Holl  [46]  for 
a 1.5  calibre  ogive  body.  However,  how  the  nuclei  interact  with  these 
viscous  effects  and  cause  the  cavitation  to  occur  is  not  fully  under- 
stood [3].  For  example,  on  the  hemisphere  nose  body  it  was  found 

that  a.  ^ -c  whereas  on  the  ITTC  body  o.  ^ -c  [40].  No  explana- 
^ Ps  ^ ^s 

tion  for  this  behavior  has  been  proposed.  Another  unexplained  obser- 
vation concerning  the  interaction  between  the  freestream  nuclei  and 
the  viscous  flow  was  made  in  the  present  experiments,  i.e.  in  both 
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the  HSWT  and  the  LTWT,  cavitation  tests  on  the  NSRDC  body  pro- 
duced band  type  cavitation  inception  at  approximately  the  same  vahn- 
of  the  cavitation  number  (0.40  in  the  HSWT  and  0.44  in  the  LTWT). 
Based  on  the  cavitation  tests  in  the  HSWT  on  a tripped  lu  tmsphen  no.st 
body  it  is  estimated  that  the  nuclei  distribution  in  the  HSWT  is  1 to  4 
orders  of  magnitude  less  than  that  in  the  LTWT.  Yet  both  facilities 
exhibited  the  normal  band  type  of  inception  at  subcritical  Reynolds 
numbers  at  about  the  same  inception  index'.  Since  band  type  inception 
is  controlled  by  the  presence  of  a laminar  separation,  it  is  speculated 
that  the  nuclei  population  in  the  separated  bubble  is  substantially  differ- 
ent from  that  in  the  freestream.  However,  at  the  present  time  no 
measurements  of  nuclei  distributions  in  the  separated  bubble  are  avail- 
able and  the  above  speculation  must  await  experimental  verification. 

In  summary  then  the  following  observations  were  made  about 
freestream  nuclei  and  cavitation  inception: 

1.  Nuclei  populations  appear  to  have  a "universal"  size  distribution 
although  the  absolute  concentrations  may  vary  by  several  orders 
of  magnitude. 

2.  If  there  are  many  large  nuclei  in  the  freestream  they  may  substan- 
tially alter  the  fully  wetted  flow  and  thus  change  the  value  of  the 
cavitation  number  and  the  type  of  inception. 

3.  Band  type  inception  appears  to  be  insensitive  to  freestream  nuclei 
populations  (provided  of  course  that  the  population  is  not  so  large 
as  to  remove  the  separation). 
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Chapter  VII 


SUMMARY  AND  CONCLUSIONS 


The  influence  upon  the  basic  viscous  flow  on  two  axisyniniel  ric 
bodies  of  freestream  turbulence  intensity  and  sinall  amounts  of  a dra^- 
reducing  polymer  in  the  boundary  layer  was  investigated.  J'bc  <dl<  i ts 
of  the  resulting  modifications  to  the  viscous  flow  upon  cavitation  incep- 
tion were  then  studied.  At  the  same  time  freestream  nuclei  popula- 
tions were  measured  and  preliminary  tests  carried  out  to  determine  the 
influence  of  different  distributions  on  inception. 

This  investigation  produced  the  following  results:  first,  the 
viscous  flow  on  the  hemisphere  nose  body  was  insensitive  to  the  present 
freestream  disturbances.  The  NSRDC  body,  by  comparison,  was  v'ery 
susceptible  to  the  freestream  turbulence  intensity  having  its  critical 
Reynolds  number  reduced  from  5 y 10^  at  0.04  percent  turbulence  level 
to  1.6  X 10^  at  3.75  percent.  In  the  cavitation  tests  the  NSRDC  body 
was  studied  at  two  turbulence  levels  (0.04  and  0.65  percent).  It  was 
found  that,  although  the  inception  index  did  not  change,  the  developed 
cavitation  became  very  unsteady  at  the  higher  turbulence  level. 

Second,  the  injection  of  small  quantities  of  a drag- reducing 
polymer  into  the  boundary  layer  was  found  to  also  lower  the  critical 
Reynolds  number  but,  unlike  the  freestream  turbulence  level,  the  poly- 
mer was  effective  on  both  bodies.  The  "effectiveness"  of  the  polymer 
in  lowering  the  critical  Reynolds  number  was  found  to  be  approximately 
proportional  to  the  amount  of  polymer  present.  The  subsequent 
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cavitation  tests  showed  a larj>e  suppression  in  the  inception  index  an>l 
a change  from  bcuid  type  inception  to  simultaneously  occurring  travel- 
ing bubble  and  "rough"  band  type. 

Third,  the  measured  nuclei  distributions  were  loiind  to  be  in 
general  accordance  with  other  populations  in  similar  facilities.  It  wa.s 
found  that,  if  there  were  many  freestream  gas  bubbles,  the  noriTally 
present  laminar  separation  was  eliminated  and  traveling  bubble  type 
cavitation  occurred  — the  value  of  the  inception  index  then  depended 
upon  the  nuclei  population.  In  cases  where  the  separation  was  present 
it  was  found  that  the  value  of  the  inception  index  was  insensitive  to  the 
freestreaun  nuclei  populations. 

The  conclusions  are  then:  both  the  increased  freestream  turbu- 

lence level  and  the  injection  of  Polyox  remove  the  laminar  separation 
by  destabilizing  the  laminar  boundary  layer  and  promoting  an  early 
transition  to  a turbulent  non-separating  one.  However,  the  mechan- 
isms by  which  either  the  turbulence  level  or  the  polymer  interact  with 
the  laminar  boundary  layer  are  not  understood  and  thus  no  explanation 
for  the  lack  of  effect  of  the  turbulence  intensity  on  the  hemisphere  nose 
body  is  forwarded. 

On  the  present  bodies  in  the  LTWT  and  the  HSWT  the  type  of 
cavitation  and  the  value  of  the  inception  index  are  determined  by  the 
presence  of  a normally  occurring  laminar  separation.  Elimination  of 
the  separation  by  freestream  turbulence,  by  the  injection  of  Polyox  or 
by  the  presence  of  many  free  gas  bubbles  in  the  inconiing  flow  causes  a 
change  in  the  type  of  inception  from  band  to  travelling  bubble  type. 

The  value  of  the  cavitation  number  at  inception  is  then  determined  by 
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I 

the  freestream  nuclei  population. 

It  was  found  that  when  traveling  bubble  type  ineeplion  oi.  ctit  ,s 
that  it  is  very  desirable  to  have  first  a universally  accepted  definiti(ni 
of  inception  based  on  the  number  of  cavitation  events  per  unit  ti..  f per 
unit  of  surface  area  and  second  a standard  procedure  for  counting  the 
number  of  events. 

In  cases  where  the  viscous  effects  control  inception  it  is  the  | 

interaction  between  these  viscous  effects  (o.g.  laminar  separation, 
attached  transition)  and  the  nuclei  which  determine  inception.  How 
these  interactions  take  place  is  not  fully  understood;  however,  in  the 
particular  case  when  a laminar  separation  occurs,  it  is  believed  that 
the  insensitivity  to  freestream  nuclei  populations  is  due  to  a substan-  j 

tially  different  population  in  the  separated  region  where  inception  j 


begins.  This  speculation  awaits  experimental  verification. 
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Fif>.  <i8.  Schlierfti  photoyraph  of  an  example  of  stable  injeilion 
of  water  at  staj^nafion  point  of  hemii  sphere  n(jse  tjoby. 
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Fig.  29.  Schlicren  photograph  of  an  e.xamplo  i>f  unstable  injio  lion 
of  water  at  stagnation  point  of  henus))heri-  nose  body. 
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Plot  of  the  regions  of  stable  and  unstable  injection  with  water 
on  the  hemisphere  nose  body  for  two  injection  tube  diameters. 
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Fig.  32.  Schematic  drawing  of  the 
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Fig.  34.  bchematic  drawings  showing  the  definitions 
of  the  laminar  separation  dimensions. 
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Fig.  43.  Comparison  of  calculated  and  measured  values  of  the  critical 
Reynolds  number  for  the  hemisphere  nose  and  NSRDC  bodies. 
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Fig.  45  Effect  of  injection  of  water  on  the  laminar  separation 
on  the  NSRDC  body.  Flow  is  from  right  to  left. 

Reo  = 3. 2 X 10^. 

(a)  Q (injection  rate)  =0.0  ml/sec 

(b)  Q = 1 . 80  ml  /sec 

(c)  Q = 3 . bO  ml  / sec 

(d)  Q = 6 . 60  ml  /sec 

(e)  Q = 9 . 80  ml  /sec 
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Fif>.  47  Effect  of  injection  of  Polyox  (WSR  301)  on  th-,-  laminar 
separation  on  the  hemisphe rt‘  nose  body.  Flow  is 
frorti  rif^ht  to  left.  Ri'o-  3.9  X 1 0'  . 

Concentration  100  wppm. 

(a)  C (inji-ction  rale)  0.0  ml/sec,  Ci  0.0 
(h)  Q - d.  5 i ml /sec,  C;  0. 9 X 1 0 

(c)  Q 1.10  ml/sec,  C.  I . I X 10''' 

(d)  6 - 1 . 79  ml/.-ec,  C I . 7 X lO''" 

(e)  Q C 00  ml/sei  , (■  3.9x10" 


Fig.  48  Effect  of  injection  of  Polyox  (WSR  301)  on  the 
laminar  separation  on  the  NSRDC  body.  Flow 
is  from  right  to  left.  Roq  - l.b  x 10"’. 
Concentration  = 500  wppm. 

(a)  Q (injection  rale)  =0.0  ml /sec,  G = 0. 0 

(b)  Q = 0.  1 ml/sec,  G = 0. 5 x lO'" 

(c)  (!l  = 0.3ml/sec,  G=  1.5x10"'’ 

(d)  6 = 0.  5 ml/sec,  G = 2.  5 X lO"'"’ 


Effects  of  injection  of  Folyox  (WSR  30  1) 
laminar  separation  on  the  NS13DC  body, 
is  from  right  to  left.  Re^  - 3.4  x 10  . 
Conc^entration  = 500  wppm. 

(a)  Q (injection  rate)  = 0.0  tnl/st'c,  (3  - 

(b)  9 ^ 0.  1 ml/sec,  G = 0.  U.  X 1 0 

(c)  Q = 0. 3 ml/sec,  G=  1.08x10  ' 
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BODY  REYNOLDS  NUMBER -UD/z. 

Effect  of  injection  of  Polyox  'WSR  301)  on  the  length 
of  the  laminar  separation  on  the  NSRDC  body. 


HEMISPHERE  BODY 
O 50  WPPM 
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Fig.  57.  The  effect  of  freestream  turbulence  level  upon  the 
cavitation  inception  number  on  the  NSRDC  body. 


polymer  injection.  Flov.'  is  from  right  to  left.  Re;.  = 3. 

(a)  Q (injection  rate)  = 0.0  ml/sec.  , C = 0.  0,  C7  = 0.  44 

(b)  (i  = 2.40  ml/sec.  , c (concentration)  = 20  v^^^pm,  G 

(c)  Q = 2.40  ml/sec.  , c = 20  wppm,  G = 0.  34  X 10  , a 

(d)  Q = 2 00  ml/ sec.,  c = 500  wppm,  G=7.11xl0®, 
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BODY  REYNOLDS  NUMBER  - UD/Z/ 


POLYMER  INJECTION  RATE-  tnl/sec 

Fig.  62.  The  effect  of  injection  of  Polyox  fWSR  3011  on  the  cavitation 

inception  number  on  the  hemisphere  nose  body  as  a function  of 
injection  rate  for  two  concentrations.  Reg  = 7,  S y 10^, 


Fig.  63.  The  same  data  as  presented  in  Fig.  62  but  now  plotted  v^ersus  the 
parameter  G to  show  the  correlation  of  the  cavitation  suppress- 
ion with  the  amount  of  polymer  injected  into  the  boundary  layer. 


BODY  REYNOLDS  NUMBER  - UD/l/ 

ig.  64  Comparison  of  the  cavitation  suppression  results  of  the 
present  study  with  those  of  several  other  investigations 


NUCLEI  NUMBER  DENSITY  DISTRIBUTION  FUNCTION  (m: 
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Fig.  65 . 


Comparison  of  nuclei  number  density  distri- 
bution functions  derived  from  the  experimen 
tal  results  of  several  investigations. 


Schlierf-n  photOL'raphs  of  the  effeiM  of  increasin;^  numbers  of  freestreair 
separation  on  the  hemisphere  nose  boriy.  I'low  is  from  left  lo  ri^ht.  H< 


NUCLEI  NUMBER  DENSITY  DISTRIBUTION  FUNCTION  (m"'') 
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NUCLEI  RADIUS  (m) 

Fig.  67  Some  nuclei  distributions  obtained  in  the  LTWT  during 
* the  present  investigations  and  one  distribution  obtained 
in  the  NSRDC  facility  by  Peterson  [l5].  The  holo- 
graphic technique  was  used  in  both  investigations. 
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Appendix  I 


MEASUREMENT  OF  FREESTREAM  TURBULENCE 
LEVEL  IN  THE  LTWT 


I.  1 Flow  Measuring  Equipment 

The  test  section  longitudinal  velocity  component  was  measured 
using  a DISA  55D05  constant  temperature  hot  film  anemometer.  The 
anemometer  consists  basically  of  a bridge  circuit  with  a feedback  am- 
plifier. The  probe,  which  forms  one  arm  of  the  bridge,  was  a DISA 
wedge-shaped  film  type.  Two  types  of  probes  were  used:  a 90  degree 
hot  film  probe  (55R36)  and  a straight  hot  film  probe  (55R32).  Initially 
it  had  been  desired  to  measure  the  velocity  and  turbulence  profiles 
across  the  test  section.  For  this  purpose,  the  90  degree  probe  was 
purchased.  This  probe  was  mounted  in  the  top  test  section  window  and 
was  "swept”  from  the  tunnel  center-line  up  to  the  window.  However,  it 
was  found  that  the  resultant  signal  was  so  dominated  by  the  vibration  of 
the  probe  that  no  useful  information  could  be  extracted.  Subsequently 
the  straignt  probe  was  purchased  and  firmly  held  on  the  tunnel  center - 
line  by  a sting  support.  Photographs  of  the  180  degree  probe  itself  and 
the  probe  mounted  in  the  tunnel  are  presented  in  Figs.  13  and  14  of  the 
main  text. 

Normally  this  DISA  anemometer  is  a battery  operated  unit. 
However,  because  of  the  high  probe  currents  required  for  measure- 
ments in  water,  it  was  necessary  to  use  an  auxiliary  power  supply  in- 
stead of  the  batteries.  A Hewlett-Packard  constant  voltage  power 
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supply  'HP  6216A)  was  used.  The  output  signal  from  the  anemometer 
was  first  low-pass  filtered  by  a simple  RC  circuit  to  remove  high  fre- 
quency noise  from  the  anemometer  feedback  amplifier.  A circuit  dia- 
gram and  the  measured  frequency  response  of  this  filter  are  given  in 
Figs.  I.  1 and  1.2  respectively.  The  filtered  signal  was  then  divided 
into  its  D.  C.  and  A.C.  components.  A Keithley  digital  multimeter 
(model  171)  was  used  to  extract  the  D.  C.  component  while  a Hewlett- 
Packard  true  rms  voltmeter  (model  3400A)  measured  the  A.  C.  compo- 
nent. Since  it  was  very  difficult  to  follow  the  fluctuating  needle  readout 
on  the  true  rms  meter,  use  was  made  of  the  output  provided  by  this  mete  r . 
The  true  rms  meter  provides  a D.  C.  voltage  from  zero  to  one  volt 
which  is  proportional  to  the  needle  deflection.  This  output  was  connec- 
ted to  a Hewlett-Packard  integrating  D.  C.  voltmeter  (Counter-DMV 
model  5326B)  which  averaged  the  signal  over  a ten  second  time  period. 

Spectral  analysis  of  the  signal  was  performed  by  a Hewlett- 
Packard  spectrum  analyzer  (model  3580A).  The  results  were  shown 
on  a cathode  ray  tube  display  which  was  photographed  for  permanent 
record.  Also  photographed  at  the  same  time  was  the  anemometer  sig- 
nal's displayed  on  an  oscilloscope.  Examples  of  both  the  signal  and 
its  spectrum  are  presented  in  Figs.  1.3  and  1.4  respectively.  Figure 
I.  5 is  a block  diagram  of  the  turbulence  level  measuring  system  and 
Fig.  I.  6 is  a photograph  of  the  system. 

I.  2 Data  Recording 

Turbulence  levels  were  measured  at  a fixed  location  on  the 
tunnel  center-line  as  a function  of  velocity  for  a given  configuration  of 
turbulence  manipulators.  At  the  beginning  of  each  test,  the  probe 
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cold  resistance  auid  zero  flow  D.  C.  voltage  were  recorded.  If  the 
probe  was  new,  a frequency  response  adjustment  was  also  made.  The 
velocity  was  then  gradually  increased  in  increments  of  two  to  three 
feet  per  second.  At  each  velocity,  the  D.  C.  component  voltage  and  ten 
ten- second  averages  of  the  true  rms  value  of  the  fluctuating  component 
voltage  were  recorded.  Also,  the  pressure  drop  across  the  nozzle 
was  recorded.  At  the  conclusion  of  each  test  (which  usually  consisted 
of  eight  velocities),  the  probe  cold  resistance  and  zero  flow  D.  C.  volt- 
age were  checked  to  see  that  no  change  had  occurred  during  the  test. 

I.  3 Data  Reduction 

First,  the  freestream  velocity  was  calculated  from  pressure 
drop  data  across  the  nozzle.  Then,  on  log-log  paper,  the  freestream 
velocity  was  plotted  versus  the  difference  between  the  anemometer  D.  C. 
voltage  and  the  anemometer  zero  flow  D.  C.  voltage.  As  expected  from 
theory,  this  gave  a straight  line  the  slope  of  which  was  measured  direc- 
tly from  the  graph.  Next,  the  ten  recorded  true  rms  values  of  the  fluc- 
tuating component  were  averaged.  Then,  since  the  fluctuating  values 
were  small  in  comparison  with  the  D.  C.  or  mean  value,  the  turbulent 
intensity  was  calculated  by  the  following  relation; 


■/  V 
— - rnis 

U ~ AV 


X SLOPE 


where 

u^  = the  rms  value  of  the  fluctuating  component  of  the  longi- 
tudinal velocity. 

U = mean  value  of  the  longitudinal  velocity. 
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\ms  ~ average  of  the  true  rms  values  of  the  fluctuating  compo- 
nent of  the  anemometer  signal. 

AV  = the  anemometer  D.  C.  voltage  minus  the  anemometer 
zero  flow  D.  C.  voltage. 

SLOPE  = slope  of  the  log-log  plot  of  velocity  versus  /SV . 

This  value  was  then  converted  to  percent  and  plotted  versus  freestream 
velocity.  And  thus  a graph  was  obtained  for  each  configuration  ofturbu- 
lence  meinipulators . The  summary  of  the  results  are  presented  in 
Fig.  15  of  the  main  text. 

I.  4 Comparison  of  Results 

The  only  other  measurements  of  freestream  turbulence  level 
in  the  LTWT  were  carried  out  by  Vauioni  and  Brooks  [ 124].  They 
photographed  the  trajectories  of  neutrally  buoyant  particles  and  esti- 
mated the  turbulence  level  from  the  position  and  time  coordinates  of 
the  particles.  It  is  very  difficult  to  accurately  determine  the  turbu- 
lence intensity  by  this  method  and  they  estimated  that  their  measure- 
ment procedure  would  produce  a turbulence  level  of  0.  4 percent  even  if 
the  actual  level  were  zerol  However,  at  the  higher  disturbance  levels 
where  this  error  is  minimized,  good  agreement  between  the  present 
values  and  those  of  Vanoni  and  Brooks  is  obtained  if  the  injector  of  the 
neutrally  buoyant  particles  was  located  in  the  grid  and  not  in  a strut 
(see  Fig.  1.7). 

L 5.  Estimate  of  the  Error  in  the  Turbulence  Level  Measurements 

At  the  low  turbulence  level  (0.  05  percent)  it  is  believed  that 


I 


-173- 


"noise"  from  the  anemometer  feedback  amplifier  dominates,  so  no  esti- 
mate of  the  error  is  possible.  However,  at  the  higher  turbulence  levels 
where  other  errors  dominate,  an  estimate  can  be  made. 


Error  in  AV  : 

The  specifications  for  the  Keithley  multimeter  give  the  accuracy 
as  ±0.02  percent  of  reading  + 0.01  percent  of  range.  However,  this  was 
not  the  major  concern  in  measuring  the  D.  C.  voltage.  At  the  higher 
turbulence  levels  the  D.  C.  voltage  could  vary  as  much  as  0.03  volts  be- 
cause of  the  turbulent  fluctuations,  but  only  one  value  was  recorded. 
Also,  the  D.  C.  level  could  change  this  much  because  of  changing  probe 
resistance.  So  AV  could  be  in  error  by  as  much  as  ±0.  03  volts.  At 
lower  velocities  where  the  D.  C.  voltage  was  approximately  3 volts  this 
implies  an  error  of: 

QQ  X 100  = 2 percent 

while  at  high  velocities  where  the  D.  C.  voltage  was  approximately  6 
volts  the  error  would  be: 

QQ  X 100  = 1 percent. 


Error  in  V : 

rms 

Specifications  of  the  true  rms  meter  give  the  accuracy  of  the 
A.C.  to  D.  C.  conversion  as  ±5  percent  for  the  frequency  range  of  inter 
est.  The  integrating  voltmeter  is  accurate  to  ±0.04  percent.  So  the 
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limitation  here  is  the  A.C.  to  D.  C.  conversion  taking  place  inside  the 
true  rms  meter. 

Error  in  Slope: 

The  estimation  of  the  slope  was  found  to  be  not  better  than 
±5  percent.  This  figure  was  arrived  at  by  (i)  drawing  several  best  fit 
lines  and  (ii)comparison  with  the  slope  of  a least  squares  fit  line. 

So  the  error  in  turbulent  intensity  is: 


(1±0.  05)V 

X (1±0.  05)SLOPE  =»  ±12  percent 

(1±0.  02)AV 


at  low  velocities,  and: 


(1±0. 05)V 


rms 


(1±0.  01)AV 


X ( 1±0.  05)SLOPE  =9  ±11  percent 


at  high  velocities. 

The  majority  of  the  error  is  introduced  by  the  error  in  the  estimate  of 
the  slope  of  the  velocity  versus  AV  line. 

Repeatability  of  measured  turbulence  levels  was  found  to  be  bet- 
ter than  5 percent. 


V-.- 
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LIST  OF  FIGURE  CAPTIONS  — APPENDIX  I 

Fig.  I.  1.  Circuit  diagram  of  the  low-pass  filter. 

Fig.  1.2.  Frequency  response  of  the  RC  filter  as  determined  with  a 
signal  generator.  Maximum  frequency  available  from  the 
signal  generator  is  2 KUz. 

Fig.  1.3.  Example  of  the  oscilloscope  trace  of  the  aneniometer  signal. 

Freestream  velocity  is  9.33  feet  per  second  and  the  turbu- 
lence level  is  0.65  percent.  Vertical  scale  is  lOmv/cm. 
Horizontal  scale  is  10  msec/cm. 

Fig.  1.4.  Exaunple  of  the  spectrum  of  the  anemometer  signal.  Free- 
strceim  velocity  is  9.81  feet  per  second  and  the  turbulence 
level  is  0.65  percent.  Horizontal  scale  is  200  Hz/cm. 

Fig.  I.  5.  Block  diagraun  of  the  turbulence  measuring  system  and  the 
signal  reduction  equipment. 

Fig.  1.6.  Photograph  of  the  turbulence  level  measuring  equipment. 

The  anemometer  is  on  the  table  beside  the  spectrum  analyzer . 
The  equipment  in  the  rack  is  used  to  determine  the  A.  C.  and 
D.C.  components  of  the  anemometer  signal. 

Fig.  I.  7.  Comparison  of  present  turbulence  level  measurements  with 
those  of  Vanoni  and  Brooks  [ 124] . 
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Fig.  1.3  Example  of  the  oscilloscope  trace  of  the  signal 

from  the  hot  film  anemometer.  !■' reest  ream  veloc 
ity  is  9.33  feet  per  second  and  the  turl)ulenc:e  l(!vel 
is  0.65  percent.  Vertical  scale  is  10  mv/i:m 
Horizontal  scale  is  10  msec/cm. 


Fig.  1.4  Example  of  the  spectrum  of  the  aneinometer 
signal.  Free.siream  velocity  is  'ESI  )«-el  per 
8»TOtHl  and  the  (nrludence  levtd  is  O.ii*'  pel' 
cent.  Ikiri'/.onlal  scale  is  ^iOO  II/,  ^cm. 


Fig.  1.6  Photograph  of  the  turbulence  level  measuring 
equipment.  The  anemometer  is  on  the  table 
beside  the  spectrum  analyzer.  The  equipment 
in  the  rack  is  used  to  determine  the  A C.  and 
D.  C.  components  of  the  anemometer  signal. 
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Appendix  II 


TUNNEL  CALIBRATION  PROCEDURE 


1 


As  there  is  a boundary  layer  growth  along  the  tunnel  walls,  it  is 
necessary  to  calibrate  the  tunnel  both  bare  and  with  a model  installed  in 
order  to  determine  the  static  pressure  and  the  velocity  at  the  location  of 
the  model.  Figure  II.  1 is  a schematic  diagram  showing  the  location  of 
the  pressure  taps  used  to  determine  the  static  pressure  at  the  model 
location.  Pressure  p^  is  the  pressure  at  the  entrance  to  the  nozzle. 
Pressure  P2  is  the  pressure  at  the  entrance  to  the  test  section  and  is 
the  average  of  four  pressure  taps  in  a piezometer  ring.  Pressure  p^ 
is  the  pressure  at  the  model  location;  it  consists  of  a single  orifice 
0.030  inches  in  diameter  drilled  in  a brass  window  plug. 


The  calibration  consists  of  measuring  the  pressure  differences  j 

D -D.,.  P,-P  and  p - Pt  as  a function  of  tunnel  velocity.  The  differ- 
^1  ^1  ^00  00 

ence  p - Po  was  measured  by  a mercury-water  manometer;  the  differ- 

ence  p^- P by  a differential  pressure  transducer  and  the  difference 

o - p by  a water  manometer.  This  was  done  for  the  bare  tunnel  at  two 
^oo  *^2 

turbulence  levels  (0.05  percent  and  0.65  percent)  and  with  the  NSRDC 
body  installed  at  one  turbulence  level  (0.65  percent). 

These  pressure  differences  were  then  converted  into  pressure 
coefficients  as  follows: 


p - p 

Pi-P? 


, n = 1,  2 . 
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These  coefficients  were  then  plotted  versus  Pj-  p^.  the  pressure  drop 
across  the  nozzle,  and  are  presented  in  Figs.  II.  2 and  II.  3. 

We  now  list  the  working  formulas  for  the  tunnel  speed  and  cavita- 
tion number . 

The  working  formula  for  the  tunnel  speed  U is: 


which  is  obtained  by  applying  the  Bernoulli  equation  between  the  model 
location  and  the  entrance  to  the  nozzle.  By  applyim  the  Bernoulli  equa- 
tion and  continuity  across  the  nozzle,  we  get: 


where  A^/A^  is  the  nozzle  area  contraction  ratio.  Substituting  this 
into  the  equation  for  U: 


”7 

1 P1-P2  Pj-Poo  . 

1 

ip 

For  the  LTWT  (Aj/A^l 

2 

1 = 14. 535,  90 

U = 


1 

0.005) 


. - 
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The  working  formula  for  the  cavitation  number  o is: 


^oo 

° = T-rr  • 

2 pU 


where  P is  the  vapor  pressure.  Or 


Poo-Pz  . Pz-Pv 
= 1 z + T~^ 

2PU  2PU 


which  upon  substitution  for  U gives: 


a = 


-C  + 


Pz-Pv 


P,-P 


1 *^ZJ 


fc  + 0.005 

L Pi  J 


In  carrying  out  the  actual  calculations,  those  values  of 


1 


and  C corresponding  to  the  bare  tunnel  were  used.  Thus  no  account 

Pz 

of  tunnel  blockage  is  made  in  determining  either  the  cavitation  param- 
eter or  the  tunnel  velocity. 
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LIST  OF  FIGURE  CAPTIONS  — APPENDIX  II 


Fig.  II.  1 Schematic  drawing  showing  the  location  of  the  pressure  taps 
used  to  calibrate  the  LTWT. 


Fig.  II.  2 


Plot  of  the  calibration  factor 
across  the  nozzle 


versus  the  pressure  drop 


Fig.  II.  3 


Plot  of  the  calibration  factor  -C 

Pz 

drop  across  the  nozzle. 


versus  the  pressure 


Fig.  II.  1 Schematic  drawing  showing  the  location  of  the 
pressure  taps  used  to  calibrate  the  LTWT. 
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Appendix  III 


BOUNDARY  l.AYER  CALCULATIONS 


III.  1 Laminar  Boundary  Layer  Growth 

The  growth  of  the  laminar  boundary  'ayer  was  estimated  by 
Thwaites  method  which  is  detailed  in  Ref.  [ 125]  for  two  dimensional 
bodies  and  in  Ref.  [ 126]  for  axisymmetric  bodies.  Application  of  the 
Thwaites  method  requires  the  evaluation  of  the  parameter  X which  is 
defined  as: 


e^u' 

e 


where  0 is  the  boundary  layer  momentum  thickness,  U'^  is  the  local 
velocity  gradient  at  the  edge  of  the  boundary  layer  and  v is  the  kine- 
matic viscosity.  It  is  determined  by  evaluating  the  following  expres- 
sion: 


X = 


0. 45  U' 
e 

2"6’“ 


r®  2 ‘ 


ds 


(III.  1) 


Where  s is  the  coordinate  measured  along  the  surface  of  the  body,  r^ 
is  the  local  radius  of  revolution  and  U^  is  the  local  velocity  at  the 
edge  of  the  boundary  layer  — s and  r^  are  shown  defined  in  F ig.  III.  1 
For  the  purposes  of  evaluation,  Eq.  UL  1 was  made  non-dimensional: 
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(-dC  /d(s/D)) rJ[l-C  ]^'^d(s/D)  ail-i) 

■I  p ayrr^  (1-c  )^‘*0  0 P 

0 P P 


where  the  velocity  has  been  non-dimensionalized  by  which  is  the 

velocity  far  from  the  body  and  the  distances  have  been  non-dimension- 

alized  by  the  body  diameter  D.  Also,  the  relationship 

where  C is  the  local  pressure  coefficient  and  is  defined  by 
P 

Cp=  (P  - Pqq)/(1/2  pU^)  has  been  used/  is  the  static  pressure  far 

from  the  body  cuid  p is  the  liquid  density. 

The  necessary  pressure  distribution  was  computed  numerically 
with  a Douglas -Neumann  potential  flow  program.  Pressure  distribu- 
tions calculated  with  this  program  for  the  hemisphere  nose  body  and  the 
NSRDC  body  are  presented  in  Fig.  III.Z.  Equation  (III.  2)  was  then  evaluat- 
ed on  a high  speed  computer  using  Simpson's  rule  to  evaluate  the  integral 
and  a spline  interpolating  routine  to  determine  derivatives  of  the  pres- 
sure distribution.  Then,  knowing  X and  0 as  a function  of  s/D,  ref- 
erence to  the  table  on  page  306  of  [ 126]  provides  values  of  the  shape 
factor  Hg  as  a function  of  s/D  which  thus  provides  the  displace  thick- 
ness  6 . 

The  position  of  laminar  separation  was  determined  by  finding 
the  value  of  s/D  at  which  X = -0.09.  Comparison  of  the  calculated 
position  of  separation  and  the  observed  position  for  the  hemisphere  nose 
and  NSRDC  bodies  can  be  made  in  Fig.  42  of  the  main  text.  The  maxi- 
mum error  is  about  3 percent. 
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III.  2 Approximate  Calculation  of  Spatial  Amplification  of  Boundary 
Disturbances 

Calculation  of  the  position  of  transition  on  the  present  axi sym- 
metric bodies  was  carried  out  by  using  the  semi  - empirical  method  pro- 
posed by  Smith  and  presented  in  Ref.  [45]  . Essentially  Smith's  method 
states  that  transition  from  a laminar  to  a turbulent  boundary  layer  occurs 
when  the  amplitude  of  aui  initial  boundary  layer  disturbance  has  been  in- 

7 9 

creased  by  a factor  of  between  e and  e . In  the  present  calculations 
7 

the  factor  e was  used. 

The  spatial  amplification  of  an  initial  boundary  layer  disturbance 
was  calculated  from  the  spatial  stability  charts  computed  by  Wazzan 

et  al  [127]  for  Falkner-Skan  velocity  profiles.  The  results  of  their  cal- 

* 

culations  are  presented  in  the  form  of  curves  of  constant  a^6 
a plot  of  dimensionless  frequency  (u)  ) versus  Reynolds  number 
where  is  the  imaginary  part  of  the  dimensional  wave  number  asso 
elated  with  a Tollmien-Schlichting  disturbance.  There  is  a plot  for 
each  Falkner-Skan  profile  (i.e.  for  each  value  of  the  Hartree  P)  used 
in  the  calculation. 

Then,  if  s^^  denotes  the  position  of  neutral  stability  (the  posi- 
tion from  which  disturbances  begin  to  grow)  the  amplification  (A)  at 
any  point  downstream  of  Sj^  is  given  by: 

s/D 

InAl  = - f (a,6*/Re^)(U/U  )(Re  )d(s/D)  (III.  3) 

L/n  ‘ 6*  e 00  D 

®N/D 

for  a fixed  value  of  the  reduced  frequency  w^- 


111,2.1  Present  Computational  Method 


Fix  the  Reynolds  number  Re^^.  From  the  Thwaites  calculation, 

the  local  Reynolds  number  based  on  displacement  thickness  Re  ...  and 

6 

the  shape  pareuneter  are  known  as  a function  of  s/D.  Correspond- 
ing to  each  a value  of  P is  obtained  from  Fig.  III.  3.  Fix  a value 

of  the  reduced  frequency  lu^.  Then  from  the  stability  chart  in  Ref.  [127] 
corresponding  to  the  calculated  P and  fixed  U)  . ^ value  of  cr.  6 /Re 

is  obtained  at.  each  value  pf  s/D.  And  thus  the  value  of  the  integrand  in 
Eq.  (III.  3)  is  known  as  a function  of  s/D  and  the  integration  may  be 
carried  out  numerically.  This  has  to  be  done  for  different  values  of 
(j)  until  the  value  producing  the  maximum  amplification  at  the  desired 
position  is  found. 

Some  of  the  calculations  made  by  this  method  on  the  hemisphere 
and  NSRDC  bodies  are  presented  in  Figs.  III.  4 and  III.  5. 
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LIST  OF  FIGURE  CAPTIONS  — APPENDIX  III 


Fig.  III.  1 Schematic  diagram  showing  the  definitions  of  the  body  coor- 
dinates s and  Tq. 

Fig.  III.  2 Pressure  distributions  for  the  hemisphere  nose  and  NSRDC 

bodies  calculated  by  the  Douglas -Neumann  potential  flow  pro- 

5 

gram  for  a body  Reynolds  number  of  4.  17x10  . 


Fig.  III.  3 Relationship  between  the  Hartree  P and  the  shape  param- 
eter H . 
s 

Fig.  III.  4 Estimated  maximum  disturbance  amplification  as  a function 
of  the  streamwise  coordinate  for  several  Reynolds  numbers 
on  the  hemisphere  nose  body. 


Fig.  III.  5 


Estimated  maximum  disturbance  amplification  as  a function 
of  the  streamwise  coordinate  for  several  Reynolds  numbers 


on  the  NSRDC  body. 
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